Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,416 T→C 5.4% intergenic (+55/‑153) 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,4160TC 5.4% 792.6 / 21.8 264intergenic (+55/‑153)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (115/122);  new base C (9/5);  total (132/132)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 268 total aligned reads displayed.

AGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGC  >  NC_002947/4901270‑4901565
                                                                                                                                                  |                                                                                                                                                           
agcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCgcgca                                                                                                                                                          <  2:2985206/150‑2 (MQ=255)
 gcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                                                         <  1:175137/150‑1 (MQ=255)
   aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCgg                                                                                                                                                       <  2:4009886/150‑1 (MQ=255)
   aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTAAAACCGCGCCGACAGGCGCGGACCGGGCGCgg                                                                                                                                                       >  1:424158/1‑150 (MQ=255)
      aTTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGCCCGTGCGCGGTcc                                                                                                                                                    >  2:4753859/1‑150 (MQ=255)
         ttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                             >  1:2221632‑M1/94‑150 (MQ=255)
         gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                             <  1:1889185/150‑1 (MQ=255)
            cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTa                                                                                                                                        <  1:6451870/150‑1 (MQ=255)
            cGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTa                                                                                                                                        <  2:6410498/150‑1 (MQ=255)
             agggtgttggggcagattgccagcgacgataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAg                                                                                                                                       >  1:1403676‑M1/90‑150 (MQ=255)
               ggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGa                                                                                                                                     >  1:3404162‑M1/88‑150 (MQ=255)
                gtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAg                                                                                                                                    >  1:5240021‑M1/87‑150 (MQ=255)
                 tgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGc                                                                                                                                   >  2:182309‑M1/86‑150 (MQ=255)
                 tgttggggcagagggccagcgactatcacgggtgggtgtagacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGc                                                                                                                                   >  2:4228274‑M1/86‑150 (MQ=255)
                  gttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCg                                                                                                                                  >  2:2086741‑M1/85‑150 (MQ=255)
                   ttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCgg                                                                                                                                 >  1:6017856‑M1/84‑150 (MQ=255)
                   ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCgg                                                                                                                                 <  1:700696/150‑1 (MQ=255)
                    tggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCggg                                                                                                                                >  1:6321769‑M1/83‑150 (MQ=255)
                    cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCggg                                                                                                                                <  2:1459117/150‑1 (MQ=255)
                     ggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGt                                                                                                                               >  1:3934811‑M1/82‑150 (MQ=255)
                      caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGGGGACCGCGCACGGTccgcgc                                                                                                                                                <  1:1940715/138‑5 (MQ=255)
                      caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtt                                                                                                                              >  1:4047014‑M2/1‑148 (MQ=255)
                       ggcagagtgccagcgactataccgggtgggtttatacgggttaaagaacaaaaaatgcttgggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGGCC‑‑‑CCC‑‑GTAGGAGCGGGttt                                                                                                                             >  2:6627537‑M1/80‑150 (MQ=255)
                       aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGttt                                                                                                                             >  2:1559366‑M2/1‑147 (MQ=255)
                         gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGAtttac                                                                                                                           >  2:502363‑M2/1‑145 (MQ=255)
                         cagattgccagcgactataccgggtgggtttatacgggttacagaacaaaaaatgctttggttatgcttctgcctgtGCCGACCTCTTCGCGGGTACATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTAc                                                                                                                           >  2:2629450‑M1/78‑150 (MQ=255)
                          aagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGAtttacc                                                                                                                          >  1:946898‑M2/1‑144 (MQ=255)
                               gccagcgactataacgggtgggtttatacgggttaaagaacaaaaactgctttggttatgcttgtgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCATACAGGCGCGGACGGTGCGAGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGa                                                                                                                     >  1:5761330‑M1/72‑150 (MQ=255)
                               ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcga                                                                                                                     >  1:2733610‑M2/1‑139 (MQ=255)
                               ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCCCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑CTAGGAGCGGGtttacccgcga                                                                                                                     <  2:953578‑M2/150‑12 (MQ=255)
                                ccagcgactataacgggtgtgtttattcgggttaaagaacaaacaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGaa                                                                                                                    >  2:1802568‑M1/71‑150 (MQ=255)
                                  cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaaga                                                                                                                  >  2:2914004‑M2/1‑136 (MQ=255)
                                  cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaaga                                                                                                                  >  2:3566668‑M2/1‑136 (MQ=255)
                                   gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGAtttacccgcgaagag                                                                                                                 >  1:5207203‑M2/1‑135 (MQ=255)
                                    aGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATACGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagagg                                                                                                                <  2:6698582‑M2/150‑17 (MQ=255)
                                                              tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                                                      >  2:2218101/1‑150 (MQ=255)
                                                                   gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATACGCTCCAACAGGCGCGGACCCTGCGCGGTCC‑‑‑CCC‑‑GGAGGAGCGGGGTTACCCGCGAAGAGGCCGGCACACGCAGAAGCATGACCAAAGCAtt                                                                                 >  2:2875483/1‑150 (MQ=39)
                                                                          ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACACGCAGAAGCATAACCCAAGCATtttttttt                                                                          >  2:5252421/1‑150 (MQ=255)
                                                                             tttCCCAACGCGTGATCGCGCGTGCGCGGGCCTGTTCGCGGGTAAATCCGCACCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCG‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                                                       <  1:5551522/150‑1 (MQ=38)
                                                                             tttCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGGAAATACGCTCCAACCGGCGCGGACCCTGCGCGGTCC‑‑‑CCC‑‑CTAGGAGCGGGTTTACCCCCGCACAGGCCGGCCCAGGCAGAAGCATAACCCAAGCATTTTTTGTTAtt                                                                       >  2:1935089/1‑150 (MQ=255)
                                                                                   aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCCCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGCAGAGGCCGGCCCAGGCAGAAGCCTAACCAAAGCATTTTTTGTTATTTTAccc                                                                 >  2:5770560/1‑150 (MQ=38)
                                                                                          gATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGTAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                                          <  1:4464565/150‑1 (MQ=255)
                                                                                           aTAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAACGCATTTTTTGTTCTTTAACCCGTATAAac                                                         <  2:6135767/150‑1 (MQ=255)
                                                                                                 ccTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGCCCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccaccc                                                   <  1:6728924/150‑1 (MQ=255)
                                                                                                 ccTGCGCCGGCCGCTTCGCGGGTAAATACCCTCCGACAGGCGCGGCCCGGGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCAGGCAGAATCCTAACCAAAGCAATTTTTGTTCTTTAACCCGGATAAacccaccc                                                   >  2:3755155/1‑150 (MQ=37)
                                                                                                   tGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGt                                                 <  2:5553733/150‑1 (MQ=255)
                                                                                                     cgcgGGCCTCTGCGGGGTAAAAGCCGCACCCCCCGGCGCGGACGGGGCGGGGCCC‑‑‑CCC‑‑GTCGGAGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                                               <  1:1994034/150‑1 (MQ=255)
                                                                                                        CGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAg                                            <  2:5627076‑M1/150‑1 (MQ=255)
                                                                                                          GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                                          <  1:1104050‑M1/150‑1 (MQ=255)
                                                                                                          GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                                          <  1:4746502‑M1/150‑1 (MQ=255)
                                                                                                            CTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc                                        <  1:2265289‑M1/150‑1 (MQ=255)
                                                                                                             tctTCGGGTGTAACTCCGCTCCCCCTGGGGGGCCCCGTGGGGGGCCC‑‑‑CCG‑‑TTCGGTGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt                                       <  1:3019206/150‑1 (MQ=255)
                                                                                                             TCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTGTTCTTTAACCCGTATAAACCCACCCGGTATAGTCGCTg                                      <  2:6383067‑M1/150‑1 (MQ=255)
                                                                                                               TTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg                                     <  1:3795732‑M1/150‑1 (MQ=255)
                                                                                                                TCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGc                                    <  2:12143‑M1/150‑1 (MQ=255)
                                                                                                                  GCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTa                                                                                                                                        <  2:5876654‑M1/48‑1 (MQ=255)
                                                                                                                      GTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg                              <  1:4804731‑M1/150‑1 (MQ=255)
                                                                                                                      GTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCT‑‑GTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTg                              <  2:3313628‑M1/150‑1 (MQ=255)
                                                                                                                       TAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc                             >  2:5058973‑M1/1‑150 (MQ=255)
                                                                                                                       TAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc                             <  2:6865378‑M1/150‑1 (MQ=255)
                                                                                                                        AAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGcc                            <  1:2706965‑M1/150‑1 (MQ=255)
                                                                                                                          aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcacaggcactatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgggggtcttcttgt                          <  2:4182500‑M2/150‑103 (MQ=255)
                                                                                                                           tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                         <  2:4115953‑M2/150‑104 (MQ=255)
                                                                                                                           tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:1265143/1‑150 (MQ=255)
                                                                                                                           tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  1:1017357/1‑150 (MQ=255)
                                                                                                                           tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                         >  2:822864/1‑150 (MQ=255)
                                                                                                                           tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                         <  2:712623‑M2/150‑104 (MQ=255)
                                                                                                                           tCCGCTCATACAGGCGCTGACCGTGCGCGGTCC‑‑‑CCC‑‑GAAGGAGCGCGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGATATAGTCTCTGGCAATATGCCCCa                         >  2:2322260/1‑150 (MQ=255)
                                                                                                                            ccGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa                        >  1:5605391/1‑150 (MQ=255)
                                                                                                                              gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                      <  2:3142060‑M2/150‑107 (MQ=255)
                                                                                                                              gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCAt                                                                                             <  1:2151497/79‑1 (MQ=255)
                                                                                                                              gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca                      >  2:3621460/1‑150 (MQ=255)
                                                                                                                              gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                      <  1:178421‑M2/150‑107 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  2:2511667‑M2/150‑109 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  1:4684810‑M2/150‑109 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  1:359621‑M2/150‑109 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  2:3466489‑M2/150‑109 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                    >  1:3117639/1‑150 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                    >  2:5672553/1‑150 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  2:4599358‑M2/150‑109 (MQ=255)
                                                                                                                                tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  2:5794270‑M2/150‑109 (MQ=255)
                                                                                                                                 ccTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggt                   <  2:4940375‑M2/150‑110 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                  <  2:6419533‑M2/150‑111 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                  <  2:5429833‑M2/150‑111 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCt                  >  1:5399210/1‑150 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtc                  <  1:6254944‑M2/150‑111 (MQ=255)
                                                                                                                                  cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgcc                                                   <  2:904762‑M2/117‑78 (MQ=255)
                                                                                                                                   tACGGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCtt                 >  2:3242106/1‑150 (MQ=255)
                                                                                                                                   tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCGCCTGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                    <  2:2097896‑M2/150‑109 (MQ=255)
                                                                                                                                   tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg                                    >  1:2430284‑M2/1‑39 (MQ=255)
                                                                                                                                    actggCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTc                <  1:5639780/147‑1 (MQ=255)
                                                                                                                                    aCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTc                <  2:6019068/150‑1 (MQ=255)
                                                                                                                                     cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg               <  2:6471695‑M2/150‑114 (MQ=255)
                                                                                                                                      aGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGGTGACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCCAAGCCTGTTTTGGTCTTTAACCCGGCTAACCCCACCCGGGATAGTCGCTGGCCATCTGCCCCACCACCCTTCac              >  2:4398506/1‑149 (MQ=255)
                                                                                                                                        gcgcGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaa            <  1:5707468‑M2/150‑117 (MQ=255)
                                                                                                                                         gggggACCGCGCACGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGaa           >  2:3054096/4‑150 (MQ=255)
                                                                                                                                         cgcgGACCGTGCGCGGTCC‑‑‑CCG‑‑TTTGGCGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGaa           <  1:2543246/150‑1 (MQ=255)
                                                                                                                                            ggcccGGGCGCGGTCC‑‑‑CCC‑‑TTTGGTGCGGGTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACaa        <  1:3482183/147‑1 (MQ=255)
                                                                                                                                            ggACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACaa        >  1:67772/1‑150 (MQ=255)
                                                                                                                                                  tGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGc  >  2:4825803/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
AGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCC‑‑GTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGC  >  NC_002947/4901270‑4901565

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: