Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,439 G→A 9.0% intergenic (+78/‑130) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,901,4390GA 9.0% 620.3 / 38.1 217intergenic (+78/‑130)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (122/70);  new base A (7/12);  total (130/87)
Fisher's exact test for biased strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCGC  >  NC_002947/4901290‑4901588
                                                                                                                                                        |                                                                                                                                                     
tggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                                       <  2:4350947‑M1/68‑1 (MQ=255)
cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                                       <  1:4999977/150‑1 (MQ=255)
 ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                                      >  2:3576001‑M1/82‑150 (MQ=255)
 ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                                      >  1:21667‑M1/82‑150 (MQ=255)
 acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt                                                                                                                                                      >  1:716055‑M2/1‑149 (MQ=255)
 acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCCGGCGCGGGCC‑‑‑CCCGGAGGAGCGGGt                                                                                                                                                      <  1:7420815‑M2/150‑2 (MQ=255)
  gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                                                                     >  1:5130600‑M1/81‑150 (MQ=255)
  gggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtt                                                                                                                                                     >  1:6688621‑M1/81‑150 (MQ=255)
   ggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                                                                    >  1:1185078‑M1/80‑150 (MQ=255)
   ggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                                                                    >  2:2856017‑M1/80‑150 (MQ=255)
    gcagattgccagcgactataacgggtgggtctatccgggttacagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTa                                                                                                                                                   >  2:2072159‑M1/79‑150 (MQ=255)
     gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                                                  <  2:6429822‑M2/150‑6 (MQ=255)
     gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttac                                                                                                                                                  >  2:1423823‑M2/1‑145 (MQ=255)
      aagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacc                                                                                                                                                 >  1:4778808‑M2/1‑144 (MQ=255)
       gattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAccc                                                                                                                                                >  2:6500921‑M1/76‑150 (MQ=255)
         ttgccagctactataacgggtggttttataagggttaaagaacaataaatgctttggttatgcttctgcctgtGCCGGCCTCTTCTCGTTTAAATCCGCTCCTACAGGCGCGGACCGTTCGCTGTCC‑‑‑CCCGTAGGAGCGGGTTTACCcgc                                                                                                                                              >  1:560330‑M1/74‑150 (MQ=255)
         aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgc                                                                                                                                              >  1:1623208‑M2/1‑141 (MQ=255)
            ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGTGGGTTTACCCGCGaa                                                                                                                                           >  1:4627301‑M1/71‑150 (MQ=255)
            ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  2:605143‑M1/71‑150 (MQ=255)
            ccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  2:605657‑M1/71‑150 (MQ=255)
            cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAACTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                                                                           >  2:1105866‑M2/1‑138 (MQ=255)
            cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                                                                           >  2:2756453‑M2/1‑138 (MQ=255)
            cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaa                                                                                                                                           >  2:2800850‑M2/1‑138 (MQ=255)
              cgcgactataccgggtggggttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTTGGAGCGGGTTTACCCGCGAaga                                                                                                                                         >  1:7176389‑M1/69‑150 (MQ=255)
                  cAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGGTAAACCGCTCCGACCGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCGAGGGGCGGGtgtacccccgaagagggc                                                                                                                                     <  1:2309338‑M2/150‑19 (MQ=255)
                     ggggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTAGCGGGTAAATCCGCTCCGACAGGCGCGGACCGGGCGCGGGCC‑‑‑CCCCGAGGAGCGGGgttacccccgaaggggccggc                                                                                                                                  <  1:2388399‑M2/150‑22 (MQ=255)
                                gtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           <  2:5319907‑M1/80‑1 (MQ=255)
                                ctgttcgCGCGGTGGGAAGGGGTTGTTTCCCAACGCGTGATCGCGCCGGCGCCGGCCTCTGCGCGGGTAACTCCGCGCCTCCCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCactat                                                                                                                       <  2:3965269/150‑5 (MQ=37)
                                                   caacccattCCCACAGCGTGATAGCGCCTGCGCCGGCCTCATCGCGGGTAAAAACCCACCAACAGGCGCGGACCCTGCGCGGTCC‑‑‑CCCCTAGGAGCCGGTTTAACCCCCAAGAGGCCCGCACACGCAGAAGCAAAAACAAAACAattttt                                                                                                    >  1:7266087/8‑150 (MQ=255)
                                                    ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTg                                                                                                   <  2:4827378/150‑1 (MQ=255)
                                                          ttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCTCCGACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGGAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCGTAACCAAAGCATTTTTTGTTCttt                                                                                             >  1:6739018/1‑150 (MQ=39)
                                                                ccgCGAGATAGCGCCTGCGCCGGCCTCTTCGCGGTTAAATCCGCTCCTCCAGGGGCGGACCGGGCGCGGACC‑‑‑CCGGTTGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCg                                                                                       <  1:3691563/148‑1 (MQ=38)
                                                                aaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCCGGTAAATCCCCTCCTACCGGCCCGGCCCGGTCGCGGTCC‑‑‑CCCGGAGGCGCCGGTTTACCCGCGAAGAGGCCGGCACACGCCGAAGCAAACACAAAGCCTTTTTTGTTGTGTGACCCg                                                                                       >  2:1965548/1‑150 (MQ=255)
                                                                  gCGTGATAGCGCCTGCGCCGGCCTCGTCGCGGGTAAATCCGCACCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAACCGtg                                                                                     >  2:5680763/1‑149 (MQ=255)
                                                                                      GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                                                                 >  2:3465116‑M1/1‑150 (MQ=255)
                                                                                           TTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg                                                            >  2:5420674‑M1/1‑150 (MQ=255)
                                                                                                gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtggggct                                                       <  2:4158487‑M2/150‑97 (MQ=255)
                                                                                                GGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTCTAAACCCACCCGTTATAGTCGCTGGCAATc                                                       <  2:2124533‑M1/150‑1 (MQ=255)
                                                                                                 GGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGACCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                                      <  1:4313275‑M1/150‑1 (MQ=255)
                                                                                                    AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCACTCTGcc                                                   <  1:1247165‑M1/150‑1 (MQ=255)
                                                                                                    AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                                                   <  1:820177‑M1/150‑1 (MQ=255)
                                                                                                    AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                                                   <  1:1247372‑M1/150‑1 (MQ=255)
                                                                                                      aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt                                                 <  1:5370651‑M2/150‑103 (MQ=255)
                                                                                                      aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt                                                 <  2:4029381‑M2/150‑103 (MQ=255)
                                                                                                      aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt                                                 <  1:3002604‑M2/150‑103 (MQ=255)
                                                                                                      ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc                                                 <  2:7105224‑M1/150‑1 (MQ=255)
                                                                                                      ATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcccc                                                 <  2:4682320‑M1/150‑1 (MQ=255)
                                                                                                       tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                                                <  2:6011389‑M2/150‑104 (MQ=255)
                                                                                                       tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:7238308/1‑150 (MQ=255)
                                                                                                       tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  2:2233286/1‑150 (MQ=255)
                                                                                                       tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                                                >  1:1585959/1‑150 (MQ=255)
                                                                                                        ccGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa                                               <  2:6721887/150‑1 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                                             <  2:1595891‑M2/150‑107 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                                             <  1:3401979‑M2/150‑107 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGtt                                                                                                 <  2:6909942/98‑1 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccaccc                                                                          <  1:146911/121‑1 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccaccc                                                                          >  2:146251/1‑121 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                                                                      <  1:4672746/125‑1 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                                             >  1:6397836/1‑150 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                                             >  1:3592284/1‑150 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                                             >  2:5489428/1‑150 (MQ=255)
                                                                                                          gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                                             <  1:305097‑M2/150‑107 (MQ=255)
                                                                                                           cTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAacac                                            >  1:5943086/1‑150 (MQ=255)
                                                                                                            tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                                           >  1:4280952/1‑150 (MQ=255)
                                                                                                            tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                                           >  1:4235018/1‑150 (MQ=255)
                                                                                                            tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattc                                                                                               >  2:2432674‑M2/1‑42 (MQ=255)
                                                                                                              cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                                         <  2:1193255‑M2/150‑111 (MQ=255)
                                                                                                              cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgt                                                             <  2:5016054‑M2/130‑91 (MQ=255)
                                                                                                               tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                                        <  2:556221‑M2/150‑112 (MQ=255)
                                                                                                               tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                                                                      <  1:945228/120‑1 (MQ=255)
                                                                                                               tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                                        <  2:7301400/150‑1 (MQ=255)
                                                                                                               tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                                        >  2:4938361/1‑150 (MQ=255)
                                                                                                               tACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                                        <  2:2810571/150‑1 (MQ=255)
                                                                                                                aCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTc                                       >  1:551133/1‑150 (MQ=255)
                                                                                                                  aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggc                                     <  1:4982589‑M2/150‑115 (MQ=255)
                                                                                                                  aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggc                                     <  1:6034487‑M2/150‑115 (MQ=255)
                                                                                                                  aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggc                                     <  2:593157‑M2/150‑115 (MQ=255)
                                                                                                                    gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa                                   >  1:1945942/1‑150 (MQ=255)
                                                                                                                    gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa                                   >  1:1103581/1‑150 (MQ=255)
                                                                                                                     gggggACCGCGCACGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGaa                                  >  2:3773142/4‑150 (MQ=255)
                                                                                                                     gggggACCGCGCACGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGaa                                  >  2:5225737/4‑150 (MQ=255)
                                                                                                                         gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatg                              <  2:2924950‑M2/150‑122 (MQ=255)
                                                                                                                         gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACATGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAg                              <  2:4014688/150‑1 (MQ=255)
                                                                                                                         gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAg                              >  1:695868/1‑150 (MQ=255)
                                                                                                                            cgggcgcGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTCACCCGTATAAACCCACCCGGGAGAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc                           >  1:5232534/4‑150 (MQ=38)
                                                                                                                                       c‑‑‑cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCCACACCCGTCAAGAACAAGGACGCACCGTGaaa                >  1:3883083/1‑150 (MQ=255)
                                                                                                                                       c‑‑‑cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaa                >  2:288839/1‑150 (MQ=255)
                                                                                                                                       c‑‑‑cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaa                >  1:6874798/1‑150 (MQ=255)
                                                                                                                                           cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATTACCAAAGCATTTTTTGTTCTTTATCCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaa               >  2:4238500/1‑150 (MQ=255)
                                                                                                                                           cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaa               >  2:4926568/1‑150 (MQ=255)
                                                                                                                                            ccGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctcccca              >  1:4123625‑M2/1‑13 (MQ=255)
                                                                                                                                            ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                                                                 >  2:5150410/1‑99 (MQ=255)
                                                                                                                                               tAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAAct           >  1:5646236/1‑150 (MQ=255)
                                                                                                                                               tAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAAct           <  2:2690787/150‑1 (MQ=255)
                                                                                                                                                  gAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt        >  1:6295822/1‑150 (MQ=255)
                                                                                                                                                   aGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACACGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTc       >  1:1441542/1‑150 (MQ=255)
                                                                                                                                                    gCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctccccacgcatggc      <  1:2813178‑M2/150‑146 (MQ=255)
                                                                                                                                                    gCGGGATTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGCACGCAACGTGAAAAAACTCTTCa      <  2:746982/150‑1 (MQ=255)
                                                                                                                                                       ggTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCg   <  2:2309338/150‑1 (MQ=255)
                                                                                                                                                        gtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctccccacgcatggcattg  <  1:2609789‑M2/150‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCGC  >  NC_002947/4901290‑4901588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: