Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,901,419 G→A 6.2% intergenic (+58/‑150) 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,4190GA 6.2% 953.2 / 33.9 308intergenic (+58/‑150)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (142/147);  new base A (11/8);  total (153/155)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 317 total aligned reads displayed.

AGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAC  >  NC_002947/4901270‑4901568
                                                                                                                                                       |                                                                                                                                                         
   aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCgg                                                                                                                                                       >  2:6409404/1‑150 (MQ=255)
     gttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                     >  2:5242212‑M1/98‑150 (MQ=255)
     gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                     <  2:4081032/150‑1 (MQ=255)
     gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                     <  2:6992937/150‑1 (MQ=255)
     gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTc                                                                                                                                                     <  1:3231689/150‑1 (MQ=255)
      cttGCCCGACCTGCCACAAGACGACCCACGAGCAGGGGGCGGCGTTCTGTGCTCGCGGTGGGACGGGGGTGTTTCCCAACGCGGGAGAGCGCCGGCGCCGGCCTCTTCGCGGGTAAATCCGCGCCTACCGGCGCGGACCGTGCGCGGTcc                                                                                                                                                    <  2:6319343/149‑1 (MQ=255)
      agtGCCCGACCTGCCACAAGAAGGCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc                                                                                                                                                    <  1:6522328/148‑1 (MQ=255)
       ttGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                                                <  2:3676032/150‑1 (MQ=255)
         gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                              <  2:20972/150‑1 (MQ=255)
           gaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGt                                                                                                                                            >  2:4927627‑M1/92‑150 (MQ=255)
            aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTa                                                                                                                                           <  2:2630191‑M1/60‑1 (MQ=255)
             gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                                          <  2:81161/150‑1 (MQ=255)
             gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                                          <  2:8386979/150‑1 (MQ=255)
             gACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                                          <  1:3790592/150‑1 (MQ=255)
             agggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAg                                                                                                                                          >  2:3793495‑M1/90‑150 (MQ=255)
              aCCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAgg                                                                                                                                         <  2:1804102/150‑1 (MQ=255)
               ggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGa                                                                                                                                        >  2:5182335‑M1/88‑150 (MQ=255)
                   ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGTGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                    <  2:7312676/150‑1 (MQ=255)
                   ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                                                    <  2:7349431/150‑1 (MQ=255)
                    tggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                   >  2:1224093‑M1/83‑150 (MQ=255)
                    cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                                                   <  2:537064/150‑1 (MQ=255)
                    cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTGAATACGCTCCAACACGCGCGGCCCCGGCGCGGGGC‑‑‑CCCCGAGGGGCggg                                                                                                                                   >  1:1084287/1‑150 (MQ=255)
                     acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTAAAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGGAGGAGCGGGt                                                                                                                                  <  2:1703812‑M2/150‑2 (MQ=255)
                       aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                                                >  1:5841762‑M2/1‑147 (MQ=255)
                       aagaagACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCAACAGGCGCGGACCCGGCGCGGTCC‑‑‑CCCGGAGGAGCGGGttt                                                                                                                                <  2:2221498‑M2/150‑4 (MQ=255)
                         gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttac                                                                                                                              >  1:137058‑M2/1‑145 (MQ=255)
                         gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttac                                                                                                                              >  2:5490137‑M2/1‑145 (MQ=255)
                         gaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAAAACGCTCCGACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGgttac                                                                                                                              <  1:2779351‑M2/150‑6 (MQ=255)
                         cagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAc                                                                                                                              >  2:5851404‑M1/78‑150 (MQ=255)
                           agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttaccc                                                                                                                            >  1:958377‑M2/1‑143 (MQ=255)
                             aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgc                                                                                                                          >  1:2301224‑M2/1‑141 (MQ=255)
                             aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgc                                                                                                                          >  1:5502825‑M2/1‑141 (MQ=255)
                              cccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcg                                                                                                                         >  2:478997‑M2/1‑140 (MQ=255)
                               ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccccga                                                                                                                        <  2:7630081‑M2/150‑12 (MQ=255)
                               ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCCCGTCGCCGGGGAAAAACCGCCAACAAGCGCCGACCGGGGGCCGGGCC‑‑‑CCCGAGGAGCCGGGTTAccccccc                                                                                                                        >  1:1544290/1‑146 (MQ=255)
                                cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaa                                                                                                                       >  2:3137830‑M2/1‑138 (MQ=255)
                                 cagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGc                                                                                                                                      >  1:7613165‑M1/70‑134 (MQ=255)
                                 aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaag                                                                                                                      >  1:5240181‑M2/1‑137 (MQ=255)
                                   gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagag                                                                                                                    >  1:3542728‑M2/1‑135 (MQ=255)
                                      actataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCGGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGGCC‑‑‑CCCTTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                                 >  1:8537058‑M1/65‑150 (MQ=255)
                                            gcggcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTAGCGGGTAAATCCCCACCAACACGCGCGGACCCTGCGCGGGCC‑‑‑CCCCGAGGAGCGGGTTTACCCGCGCACAGGCCCGcccc                                                                                                           >  2:1564775/1‑147 (MQ=255)
                                                          tgggtgtgGTAATGGGTTGTTTCCCAATGCGTTTTCGCGCCTGCGCGGGCCTCTTCGGGGGTAAACCGGCTCCTACCGGCGGGGACGGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGTTTTACCGGCGAAGTGGCCGGCACAGGCAGAAGCATAAc                                                                                             <  2:821105/146‑1 (MQ=255)
                                                          tcGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCGCGTCGCGGGTTAATCCGCGCCTACAGGCGCGGACCGGGCGCGGTCC‑‑‑CCCGGAGGGGCGGGTTTAACCCCGACGAGGCCGGCACAGGCAGAAGCATAAc                                                                                             >  1:679696/1‑150 (MQ=255)
                                                              tgtgGGATTGGGTTGTTTCCCAACGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTCCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCaaa                                                                                         <  2:6264432/150‑1 (MQ=255)
                                                              tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑ccc                                                                                                                                              <  2:7523444/97‑1 (MQ=255)
                                                                        ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCCAAGCATTTTTTg                                                                               >  2:3553702/1‑150 (MQ=255)
                                                                             tttCCCAAAGCGGGATAGCGCCTGCGCCGGCCGCTTCGCGGGTAAATACGCTCCGACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTcgt                                                                          >  2:2502400/1‑148 (MQ=39)
                                                                                cccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTaa                                                                       <  2:7191857/150‑1 (MQ=255)
                                                                                 ccaccgCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAc                                                                      <  2:8813781/145‑1 (MQ=255)
                                                                                  cAATGCGTTATCGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACCGGCGCGGCCCGTGCGCGGTCC‑‑‑CCCGTAGTAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                                                     <  1:7115496/150‑1 (MQ=38)
                                                                                  cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                                                     <  2:610954/150‑1 (MQ=255)
                                                                                      gCGTGATCGCGCCTGCGCCGGCCTCTTCGCGGTTAACTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGta                                                                 <  1:7283936/150‑1 (MQ=255)
                                                                                          gATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa                                                             <  1:5098692/150‑1 (MQ=255)
                                                                                                     cgcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta                                                  <  2:8767819/150‑1 (MQ=255)
                                                                                                       cgggCCTCTTGGCGTTTAAATCCGCTCCTACCGGCGGGGCCGGTGGGGGGTCC‑‑‑CCGTTCGTTGGGGTTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                                                <  1:4075253/148‑1 (MQ=255)
                                                                                                       CCGGCCTCTTCGCGGTTAAATCCGCTCCTACCGGCGCGGACGGGGCGCGGTCC‑‑‑CCCGTTGGAGCGGTTTAACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                                                >  1:1930115‑M1/1‑150 (MQ=255)
                                                                                                          gCCTGTTCGGGGTTAATCCGGCTCGTACTGGGGCGGCCGGGGGGGGGGCC‑‑‑CCGGTGGGAGGGGGTTTCCCCGCGACGTGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc                                             <  1:4038362/150‑1 (MQ=255)
                                                                                                              CTTCGCGGTTAAATCCGCTCCTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTg                                         >  2:3202245‑M1/1‑150 (MQ=255)
                                                                                                                    cgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtc                                    <  1:2875762‑M2/150‑96 (MQ=255)
                                                                                                                    cgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtc                                    <  1:7317939‑M2/150‑96 (MQ=255)
                                                                                                                     ggttAAATCCGCTCCACCCGGCGCGGACGGGGCGGGGCCC‑‑‑CCGGTCGGTGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc                                   <  1:3078866/147‑1 (MQ=255)
                                                                                                                     ggttAAACCCGCTCCTACTGGCGCGGACCGGGGGCGGTCC‑‑‑CCGGTTGGTGGGGTTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc                                   <  1:3267321/147‑1 (MQ=255)
                                                                                                                     GGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGgcaaac                                   <  2:1849292‑M1/150‑3 (MQ=255)
                                                                                                                     GGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc                                   <  2:4040386‑M1/150‑1 (MQ=255)
                                                                                                                      GGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt                                  <  2:7726304‑M1/150‑1 (MQ=255)
                                                                                                                        tAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgggggtcttct                                <  1:5632962‑M2/150‑100 (MQ=255)
                                                                                                                            tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                            <  1:7474293‑M2/150‑104 (MQ=255)
                                                                                                                            tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtg                            <  2:1670324‑M2/150‑104 (MQ=255)
                                                                                                                            tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                            >  2:6709620/1‑150 (MQ=255)
                                                                                                                            tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                            >  2:4466239/1‑150 (MQ=255)
                                                                                                                            tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCa                            >  1:2615299/1‑150 (MQ=255)
                                                                                                                                gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                         <  2:5699117‑M2/150‑107 (MQ=255)
                                                                                                                                gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggca                         <  2:1279141‑M2/150‑107 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                       <  1:5147381‑M2/150‑109 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                       >  1:4592062/1‑150 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                       >  1:1809737/1‑150 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc                       >  1:6371855/1‑150 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                       <  2:506588‑M2/150‑109 (MQ=255)
                                                                                                                                  tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcagg                       <  2:6559363‑M2/150‑109 (MQ=255)
                                                                                                                                    cTACCGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCt                     <  2:2312005/150‑1 (MQ=255)
                                                                                                                                    cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgcccccggctcgtgggtcttcttgtggcaggtc                     <  1:4470574‑M2/150‑111 (MQ=255)
                                                                                                                                    cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTGTAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCt                     >  1:8039628/1‑150 (MQ=255)
                                                                                                                                    cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCt                     >  2:7480520/1‑150 (MQ=255)
                                                                                                                                     tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                    <  1:6118768‑M2/150‑112 (MQ=255)
                                                                                                                                     tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg                    <  1:6417481‑M2/150‑112 (MQ=255)
                                                                                                                                        aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTtata                                                <  1:556/119‑1 (MQ=255)
                                                                                                                                        aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCCGCCCCAACACCCTTCaa                 >  2:7647624/1‑150 (MQ=255)
                                                                                                                                         ggCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAg                >  2:2780627/1‑150 (MQ=255)
                                                                                                                                          gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa               >  1:8347128/1‑150 (MQ=255)
                                                                                                                                          gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTCTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGa               >  1:4330519/1‑150 (MQ=255)
                                                                                                                                          gcgcGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaa               <  1:1286818‑M2/150‑117 (MQ=255)
                                                                                                                                            gcgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc             <  1:8684131/150‑1 (MQ=255)
                                                                                                                                              ggACCGGGCGCGGTCC‑‑‑CCGGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACaa           <  1:3456259/150‑1 (MQ=255)
                                                                                                                                              ggACCGCGCACGGTCC‑‑‑GCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACaa           >  2:3402537/1‑150 (MQ=255)
                                                                                                                                                 ccGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCTCAACACCCTTCAAGAACAAGGa        <  1:4866206/150‑1 (MQ=255)
                                                                                                                                                 ccGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGa        <  2:685912/150‑1 (MQ=255)
                                                                                                                                                  cGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc       <  2:4224353/150‑1 (MQ=255)
                                                                                                                                                  cGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGAc       >  1:8756833/1‑150 (MQ=255)
                                                                                                                                                       gcgGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAc  <  1:2387455/150‑1 (MQ=255)
                                                                                                                                                       gcgGTCC‑‑‑CCCGTAGAAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAc  >  2:434883/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                         
AGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAC  >  NC_002947/4901270‑4901568

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: