Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,898,904 G→C 7.3% R662G (CGG→GGG)  PP_5173 ← RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,898,9040GC 7.3% 245.4 / 5.2 78R662G (CGG→GGG) PP_5173RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (3/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

ATTTCCACGCGCAGCACCTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAGGG  >  NC_002947/5898756‑5899049
                                                                                                                                                      |                                                                                                                                                 
aTTTCCACGCGCAGCACCTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATcggccg                                                                                                                                                   >  5:313305/1‑149 (MQ=255)
            agcaCCTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACcggcgg                                                                                                                                       <  5:52661/149‑1 (MQ=255)
              caCCTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGc                                                                                                                                     <  1:197592/149‑1 (MQ=255)
                ccTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGcc                                                                                                                                    >  1:278363/1‑148 (MQ=255)
                     cccGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCAtct                                                                                                                               <  1:361341/148‑1 (MQ=255)
                      ccGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCAtctc                                                                                                                              >  5:4754/1‑148 (MQ=255)
                       cGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGCTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTcc                                                                                                                             <  1:109155/149‑1 (MQ=255)
                         gCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGc                                                                                                                          <  1:275407/149‑1 (MQ=255)
                             gTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGAtt                                                                                                                      >  8:182815/1‑149 (MQ=255)
                                 cAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCAc                                                                                                                  >  6:217952/1‑149 (MQ=255)
                                   gTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGt                                                                                                                >  4:135014/1‑149 (MQ=255)
                                         gATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTgc                                                                                                          <  6:313305/149‑1 (MQ=255)
                                             atcTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTgccgac                                                                                                      >  1:193861/1‑149 (MQ=255)
                                                tcGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTgccgacgcc                                                                                                   >  3:282536/1‑149 (MQ=255)
                                                 cGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTgccgacgcc                                                                                                   >  6:350008/1‑148 (MQ=255)
                                                  gCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTgccgacgc                                                                                                    <  5:217952/146‑1 (MQ=255)
                                                        gtgCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAAtt                                                                                           >  1:71517/1‑149 (MQ=255)
                                                          gCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAAttct                                                                                         >  8:186575/1‑149 (MQ=255)
                                                          gCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAAttct                                                                                         >  7:130137/1‑149 (MQ=255)
                                                           cTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAAttct                                                                                         >  6:2079/1‑148 (MQ=255)
                                                             cGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCa                                                                                      <  3:86586/149‑1 (MQ=255)
                                                              gTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCa                                                                                      <  7:266719/148‑1 (MQ=255)
                                                              gTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCAt                                                                                     >  3:200571/1‑149 (MQ=255)
                                                               ttCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGAtt                                                                                                                      <  8:176315/115‑1 (MQ=255)
                                                               ttCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGAtt                                                                                                                      >  7:176315/1‑115 (MQ=255)
                                                                    gTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCa                                                                                  <  5:253728/146‑1 (MQ=255)
                                                                    gTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCa                                                                                  <  4:282536/146‑1 (MQ=255)
                                                                    gTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGa                                                                               >  5:286706/1‑149 (MQ=255)
                                                                    gTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGa                                                                               >  6:133307/1‑149 (MQ=255)
                                                                         gCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGAtcttc                                                                          <  5:143026/149‑1 (MQ=255)
                                                                           aaGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTg                                                                        <  8:29567/149‑1 (MQ=255)
                                                                            aGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGGCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGc                                                                       <  8:130137/149‑1 (MQ=255)
                                                                               tGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAgg                                                                    <  5:350008/149‑1 (MQ=255)
                                                                                 gCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATGGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCTGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCg                                                                  <  4:30427/149‑1 (MQ=255)
                                                                                 gCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCg                                                                  >  6:179459/1‑149 (MQ=255)
                                                                                 gCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCg                                                                  <  2:243365/149‑1 (MQ=255)
                                                                                   cAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGc                                                                                                                                     >  5:293835/1‑80 (MQ=255)
                                                                                   cAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGc                                                                                                                                     <  6:293835/80‑1 (MQ=255)
                                                                                   cAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCTGGATCTTCTGCAGGCGGt                                                                >  4:112274/1‑149 (MQ=255)
                                                                                   cAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAAATCTTCATGCAGGATCTTCTGCAGGAGGt                                                                >  3:134641/1‑149 (MQ=255)
                                                                                   cAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCCCTGACCCGGTACTGAAACGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGCGGCCTAGGAATTCTTCATGCAGGATCTTATGCAGGCgg                                                                 >  8:292392/1‑148 (MQ=255)
                                                                                    aGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTc                                                               <  7:241513/149‑1 (MQ=255)
                                                                                     gCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTc                                                               >  3:9715/1‑148 (MQ=255)
                                                                                      cTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCgg                                                                 <  6:373594/145‑1 (MQ=255)
                                                                                      cTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCgg                                                                 >  5:373594/1‑145 (MQ=255)
                                                                                      cTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCca                                                             <  1:182814/149‑1 (MQ=255)
                                                                                       tCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCGCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCca                                                             >  1:78030/1‑149 (MQ=255)
                                                                                       tCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCca                                                             >  1:132808/1‑148 (MQ=255)
                                                                                            ggCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAg                                                                     >  4:293233/1‑135 (MQ=255)
                                                                                            ggCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAg                                                                     <  3:293233/135‑1 (MQ=255)
                                                                                                  tttGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCACGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGa                                                  <  2:193861/148‑1 (MQ=255)
                                                                                                  tttGCGCACCTGGTCGATGTCGGCACCGCTGAACCGGTACTGAATCGGCCGGCCCACCGGCGGGCGCCTCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTGTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGAc                                                 >  2:351709/1‑149 (MQ=255)
                                                                                                    tGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACg                                                <  5:2079/148‑1 (MQ=255)
                                                                                                      cgcACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGctc                                             <  6:4754/149‑1 (MQ=255)
                                                                                                          ccTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCg                                         <  7:186575/149‑1 (MQ=255)
                                                                                                           cTGGTCGATGTCGGCAGCGCTGACCCGGTACTGAATCGTCCGGCCCACCGGCTGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTGTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTc                                          <  3:42015/147‑1 (MQ=255)
                                                                                                                tCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCCGACGCTCTTCGAAGcc                                    <  5:133307/149‑1 (MQ=255)
                                                                                                                    tGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTg                                <  5:179459/149‑1 (MQ=255)
                                                                                                                      tCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCt                              <  2:278363/149‑1 (MQ=255)
                                                                                                                        ggCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCGCAAAATTTATTCATAAAGGATCTTATGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTg                             >  1:123756/1‑148 (MQ=255)
                                                                                                                        ggCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGa                            >  2:51113/1‑149 (MQ=255)
                                                                                                                        ggCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGa                            <  7:182815/149‑1 (MQ=255)
                                                                                                                              gCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGACGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCa                      >  5:68810/1‑149 (MQ=255)
                                                                                                                               cTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCAc                     <  2:71517/149‑1 (MQ=255)
                                                                                                                               cTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCAc                     <  4:9715/149‑1 (MQ=255)
                                                                                                                                    ccGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTg                             >  6:331652/1‑136 (MQ=255)
                                                                                                                                    ccGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTg                             <  5:331652/136‑1 (MQ=255)
                                                                                                                                       cgctggagatgGGCCCGCCGGTGGGCCGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTgc              <  4:81733‑M1/138‑1 (MQ=255)
                                                                                                                                         tACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGTATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTgcgc            <  2:132808/149‑1 (MQ=255)
                                                                                                                                         ctggagatgGGCCCGCCGGTGGGCCGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCa                                                                      <  7:21368‑M1/82‑1 (MQ=255)
                                                                                                                                         ctggagatgGGCCCGCCGGTGGGCCGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCa                                                                      >  8:21368‑M1/10‑91 (MQ=255)
                                                                                                                                         ctggagatgGGCCCGCCGGTGGGCCGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCca                                                             >  5:195398‑M1/10‑100 (MQ=255)
                                                                                                                                         ctggagatgGGCCCGCCGGTGGGCCGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCca                                                             <  6:195398‑M1/91‑1 (MQ=255)
                                                                                                                                            tGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGcgcg           >  8:46980/1‑147 (MQ=255)
                                                                                                                                              aaTCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCgta         <  3:135014/147‑1 (MQ=255)
                                                                                                                                                  GGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAgg   >  2:78030‑M1/1‑149 (MQ=255)
                                                                                                                                                  GGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAgg   >  8:98042‑M1/1‑149 (MQ=255)
                                                                                                                                                  GGCCCGCCCACCGGCGGGCCCATCTCCCGCGATTGCACGGTGGTGCCGACGCCGACGAAGTCTTCATTCCCGATCTTCTGCAGGCGCTCCATCATGGCCTGACGCTCATCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAgg   >  6:34628‑M1/1‑149 (MQ=255)
                                                                                                                                                   GCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGCCGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAggg  <  2:188226‑M1/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
ATTTCCACGCGCAGCACCTTACCCGGCTCGTTCCAGTCGTAGATCATCTCGCCGATGTGCTCGTTCTGGTCGAGCAAGGTGGCCAGCTCGATGGCGTGTTTGCGCACCTGGTCGATGTCGGCACCGCTGACCCGGTACTGAATCGGCCGGCCCACCGGCGGGCCCATCTCCAGCGATTGCACGTTGGTGCCGACGCCGACGAATTCTTCATGCAGGATCTTCTGCAGGCGGTCCATCATGGCCTGACGCTCTTCGAAGCCCTTGCTGACAATCACCAGCTGCGCGTAGTAAGGG  >  NC_002947/5898756‑5899049

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: