Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,6240TG29.0% 114.9 / 10.3 69G252G (GGT→GGG) ptsGglucose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base G (17/3);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 6.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGAC  >  NC_000913/1158378‑1158867
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
cATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTCCGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGc                                                                                                                                                                                                                                                         <  1:232926/250‑1 (MQ=255)
     ggCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTg                                                                                                                                                                                                                                                   <  2:113697/251‑1 (MQ=255)
        cgccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGT‑GGGCTTCCTGTTc                                                                                                                                                                                                                                                >  1:404884/1‑251 (MQ=255)
             attgtttcttcaagccatgccTTCTCTCATTGTGCTG‑TTATCCAGAACCCGTTAGTGGCTGTTGTCATTGTGGGTTGCATCGAACGTGGCCTGGTACCGTTTGGGCTGAACCACATATGTAACGGACCTTTCCAGATGGATATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCTACATTCCGTGTTATATTGCTGGAGACCCGCCTGCGGGGAAACTGTCTGG‑TGTCTTACTGTTCaaaa                                                                                                                                                                                                                                            <  2:224039/231‑1 (MQ=255)
             aTTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTTCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAAt                                                                                                                                                                                                                                           <  1:467045/251‑1 (MQ=255)
                   tCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACg                                                                                                                                                                                                                                      <  1:465201/250‑1 (MQ=255)
                     tGCAATACAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGAACCGACTGCGGGGAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGGACGGTc                                                                                                                                                                                                                                   >  1:431140/1‑251 (MQ=255)
                       cAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTg                                                                                                                                                                                                                                 <  1:231667/251‑1 (MQ=255)
                        aaTCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGc                                                                                                                                                                                                                                <  1:176161/251‑1 (MQ=255)
                             cgACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCt                                                                                                                                                                                                                           <  1:50491/250‑1 (MQ=255)
                             aGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCt                                                                                                                                                                                                                           <  1:44967/251‑1 (MQ=255)
                                ccTTCTCGCATTGGGCTGCTTG‑CCAGATCCCGGTAGTTGCGTTTGGCATTTCCGGTTTCATCGAACGTTGCCTGGTCCCGTGTGGGCGGCCCCACATCTGGAACGTTGCTTTCCAGATGCAGATTGGTGAATACTCCAACGCTGCAGTTCCGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGGCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGcc                                                                                                                                                                                                                        <  2:201730/251‑1 (MQ=255)
                                       tcAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGCGTCCGGGGAAACTGTCTGT‑GGGCTTCCTGTTCGAAATAGACGGTCTACCAGCAGCCCCAATTg                                                                                                                                                                                                                 >  2:143146/1‑251 (MQ=255)
                                        cAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGGTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAAATGTCTGG‑GGGGTTTCCGTTCAAAATGTAAGGTCTGCCAGCCGGCGCCATTGc                                                                                                                                                                                                                >  2:179556/1‑251 (MQ=255)
                                         aGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTAAGGTCTGCCAGCTGCCGCAATTGCt                                                                                                                                                                                                               >  1:524094/1‑251 (MQ=255)
                                                 gctTT‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGGTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGc                                                                                                                                                                                                        <  1:266761/250‑1 (MQ=255)
                                                 gctTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCCCATGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCCATTGCTTTCTGGc                                                                                                                                                                                                        >  2:213120/1‑250 (MQ=255)
                                                         aGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCCGGGAAAATGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCt                                                                                                                                                                                               >  1:101653/1‑251 (MQ=255)
                                                             cccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtgg‑t                                                                                                                                                                                                                                                            <  2:229877/187‑1 (MQ=255)
                                                                                aTTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGTAAACTGTCTGG‑GGGCGTCCCGTTCAAAATATACGGTCTGCCAGCTGCCGCCATTGCTATCCGGC‑ACCCTGCTAAACCAGAA‑AACCGC‑GCGGAAgg                                                                                                                                                                      >  2:13524/1‑250 (MQ=255)
                                                                                         ttCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑ta                                                                                                                                                             <  1:104794/250‑1 (MQ=255)
                                                                                           cATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑tatt                                                                                                                                                           <  1:3333/250‑1 (MQ=255)
                                                                                            aTCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCAAATGCGGGTAAAATGTCTGG‑GGGCTTCCTGTTCAAAAAGTACGGGCTGGCAGCTACCCGAATTGGTAACTGGC‑ACTCTGCGAAACCCAAA‑AACCGCCGCGaagggggcggggttat                                                                                                                                                           >  2:346432/1‑237 (MQ=255)
                                                                                               gAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGAt                                                                                                                                                      <  1:172979/251‑1 (MQ=255)
                                                                                                 aCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGGACCGTCTGCCAGCTGCCGCAATTGCTATTTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGAtct                                                                                                                                                    >  1:131745/1‑251 (MQ=255)
                                                                                                  cGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCAACTGCGGGTCAACTGGCTGG‑GGGCTTCCTGTTCAAAATGGACGGTCCGCCCGCTGCCGCACTTGCTATCTGGC‑ACTCTGCTAAAACACAA‑AACACC‑GCGGAAGTgggggggttttttatttc                                                                                                                                                    >  2:52788/1‑236 (MQ=255)
                                                                                                        cTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCCATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCCG‑GAATATGATCTCCgaggc                                                                                                                                             >  1:104568/1‑251 (MQ=255)
                                                                                                         tGGTACCGTTTGGCCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGc                                                                                                                                                                                                                            <  2:105181/175‑1 (MQ=255)
                                                                                                         tGGTACCGTTTGGCCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGc                                                                                                                                                                                                                            >  1:105181/1‑175 (MQ=255)
                                                                                                                     gTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAAATGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGGAAGGGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTc                                                                                                                              >  1:228178/1‑251 (MQ=255)
                                                                                                                       cTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCCGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GGGAAAGGGGGCGG‑GATTATGATTTCCGCGGCGC‑TGAcccctttccc                                                                                                                             >  2:148626/1‑243 (MQ=255)
                                                                                                                         gCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTg                                                                                                                           <  1:116976/250‑1 (MQ=255)
                                                                                                                                   tGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGAAATTCCGCGTTATATGGCGGGAGACACGCCTGCGGGGAAACTGTCTGG‑TGGGTTCATGTTCACAATGTACGGTCTGCCAGCTGACGTAAATGCTATCTCGA‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAACAAGGGCGG‑GATTATGATCTCCGCGGCTC‑TGACCC‑CGGTCCTGACCGGTATCAc                                                                                                                >  1:541596/1‑251 (MQ=255)
                                                                                                                                     gcaCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCCACTGCGGGTAAACCGTCCGG‑TGGCTTCCTGTTCAAAATGTACGGTCTTCCCGCTGCCCCCATTGCTATCTGGC‑CATCTGCTCAACCAGAA‑CACCGC‑GCGCAAGTGGGCGG‑TATTCTGATTTCCGCGGCGC‑TGACCT‑CGGTCCTggcctgtgtccccc                                                                                                              >  2:192542/3‑238 (MQ=255)
                                                                                                                                      aaCGTACATTTCCAGATGCACATTGGTGAATACACCAACGCAGAAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCAGATGACCCGACTGCGGGTAAACTCGGTAG‑TGGCTTCCTGTTCGAAATATACGGTCTGCCAGCTGCCGAAATTTCTATCTTGC‑ACTCTTCTAAACCCGAA‑AACCAC‑CCGAAAGTGGGCGG‑AGTTATTATCTCCGCGGCGC‑GGCACT‑CGGTCCTGACCGGGATAACCgg                                                                                                             >  1:391908/1‑250 (MQ=255)
                                                                                                                                         gTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑GTGCCTCCCGGTCAAAATGAACCGTCTGCCAGCAGCCGCACTTGCTATCCGGCCAC‑CTGCTAAACCAGAACAACC‑C‑GCGACAGAGGGCGGGCTTTATG‑TATCCGCGGCGC‑TGACCTTCGTTC‑GGACCGGTACCACCgcggc                                                                                                          >  2:6815/1‑247 (MQ=255)
                                                                                                                                          tACCTTTCCAGATGCAGATTGGTGACTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGCCCCGACTGCGGGGAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTCCGGTCTGCCAGCTGCCGCAGTTGCTATCTGGC‑CCTCCGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TCTTATGCGCTCCGCGGCGC‑NGACCT‑CGTTCCTGACCGGTATCAcctcgccg                                                                                                         >  2:239797/1‑245 (MQ=255)
                                                                                                                                           aCCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAAATGGCTTG‑GGGCTTCATGTTCAAAATGCACGGTCTGCCAGCTGCCGCCATTGCTATCTTGG‑ACTCCGCTAAAACAGAA‑AACCGC‑GCGGAAGCGGGCGG‑GATTATGAACTCCGCGGCGC‑CNACCT‑CGTTCATGACCGGTATCACCGAGcagg                                                                                                        >  2:264704/1‑248 (MQ=255)
                                                                                                                                                      tGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑TGGGTTCCTGTTCAAAATGTAAGGTCTGCCAGATTCCGCAATTGCTATCTGGC‑ACTCTGCTCAACACGAA‑AACCGC‑GCGAAAGGGGGCGG‑TATTCTTATCTCCGCGGCGC‑TGACCT‑CGTTCCTCCACGGAATCCACGAGCCGATCGTG‑TTCTcc                                                                                            >  2:365289/1‑251 (MQ=255)
                                                                                                                                                          gATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑TGGCCTCCTGTTCAAAATGTACGGTCTTCCCGCTGCCGCAATTGCTATCTGGC‑CCTCTGCTAAACCAGAA‑AACCCC‑GCGAAAGTGGGCGG‑GATTATAATCTCCCCGGCGC‑TGACCC‑CGCTCCTGCCCGGTATCACCGAGCCGATCGAG‑TTCTCTTTc                                                                                         >  1:23851/1‑250 (MQ=255)
                                                                                                                                                          gATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGAACCGACTGCGGGTAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCCGCTGCCGCCATTGCTATCTGGC‑ACTCTGCCAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑GATTATGATCTCCGCGGGGC‑TTAACT‑CGTTCCTGGCCGGGATccacgcggcgaacggggtttcctctt                                                                                         >  1:540334/1‑226 (MQ=255)
                                                                                                                                                              ggTGAATACCCCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGt                                                                                     <  1:279070/250‑1 (MQ=255)
                                                                                                                                                               gTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAGCTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTc                                                                                   <  1:52788/251‑1 (MQ=255)
                                                                                                                                                                         cAACACAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAACTGTCTGGGTGG‑TTCCTGTTCCAAATGGCCGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑CCACTGCTTAACCCGAA‑AACCGC‑GCGACAGGGGGCGG‑TATTATGTTCTCCGCGGCGC‑TGGACT‑CGTTCCTGATCGCGATCACCCAGCCGAACGCG‑TTCTCCTTCCTGTTCTGTGCGccct                                                                         >  2:319538/1‑249 (MQ=255)
                                                                                                                                                                           aCGCCGCAGGGCAGGTTTCCCACGGCGACATTCCGCGTTATAGGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAACCCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATc                                                                       <  2:224470/251‑1 (MQ=255)
                                                                                                                                                                            cGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGCACCGACTGCGGGGAAAACGTCTGG‑GGGCCTCCTGTTCAAAATGTAAGGTCTGCCAGCTGCCGCCATTGCTATCTGGC‑ACTCTTCTAAACCAGAA‑CAACGC‑GCGGAAGGGGGCGG‑GATTATGATCTTCGCGGCGC‑TGACCT‑CGTTCCTTACCGGTATTCCCGCGCCGCTCGGG‑TTCTTCTTCATTCTCCGTGCGCCGATcc                                                                      >  2:115594/1‑251 (MQ=255)
                                                                                                                                                                                 caggtcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTAc                                                                 <  1:536863/251‑1 (MQ=255)
                                                                                                                                                                                   ggtcaggtTTTCCACGGCGACTTTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAATCTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTTCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCCCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTAcat                                                               <  2:137473/251‑1 (MQ=255)
                                                                                                                                                                                    gtcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTCCGGGGAACCTGTCTGG‑GGGCTTCCTGTTCAAAATGTACCGTCTGCCAGCAGCCGCCATTGCTATCTTGC‑ACTCTGCTAAACCCAAA‑AACCGC‑GCGGAAGGGGGCGG‑GATTATGATCTCCGCGGAGCCTG‑CCC‑CGTTACTGACCGGGATCAACCTGCCGACCGGGGTTC‑CCTTCCCGTTCGTTGCGCCGGTCCCGTAcaa                                                               >  2:434898/1‑249 (MQ=255)
                                                                                                                                                                                               cACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGAt                                                   <  1:148626/251‑1 (MQ=255)
                                                                                                                                                                                               cACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑GATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑GTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGAt                                                   >  1:108949/1‑251 (MQ=255)
                                                                                                                                                                                                aCGGCGACATTCCGCGTTGTATGGCGGGTGACCCGAATGCGGGGAAACTGTCTGG‑TGGCTTCCTGTTCCAAATGTACGGTCTTCCAGCTGCCGCAATTGCTCTCTGGC‑AATCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑GATTATTATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCgggat                                                                                                                   >  2:442284/1‑184 (MQ=255)
                                                                                                                                                                                                aCGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTAt                                                                                                                   <  1:442284/187‑1 (MQ=255)
                                                                                                                                                                                                 cGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑CCTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑GATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGGAGCACCGAGCCGATCGAG‑TTCTCCCTCCTGTTCGTTTCGCCGATCCTGTACAtcctcccccccgttc                                                 >  2:65747/1‑242 (MQ=255)
                                                                                                                                                                                                   gCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAAATGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTTACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACCTCATCCACGCGATTCTg                                               >  1:488154/1‑251 (MQ=255)
                                                                                                                                                                                                     gACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGc                                             >  1:342856/1‑251 (MQ=255)
                                                                                                                                                                                                      aCATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGGGGGCGG‑TATTATGATCTCCGCGGCGC‑TGGCCT‑CGTTCCTGACCGGTATCACCGAGCCCATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGcg                                            >  1:187579/1‑250 (MQ=255)
                                                                                                                                                                                                           ccGCGTTATATGGCTGGTGACCCGACTGCGGGTAAACTGTCTGG‑GGTCTTCCTGTTCAAAATGTACGTTCTGCCAGCTGCCGCAATTGCTATCTGGC‑TCTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGTCGC‑TGACCT‑CGTTCCTGACCGGTATCCCCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGcc                                        <  1:216801/250‑1 (MQ=255)
                                                                                                                                                                                                            cgcgTTATATGGCGGGTGACCCGGCTGCGGGGAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAAACAGAA‑AACCGC‑GGGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGAACT‑CGTTCCTGACCGGGATCAACGAGCCGATAGAG‑GTCTCCTTCATGTTCGTTGCGCCGGTCCTGTACATCATCCACGCGGTTCTGGCAGGCCTg                                      >  1:409012/1‑251 (MQ=255)
                                                                                                                                                                                                              cgTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑GGGCTTCCTGTTCAAAATGTACGGGCTGCCCGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑GATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGCG‑TTCTCCTTCCTGTTCGTTGCGCCGCTCCTGTACATCATCCACCCGATTCTGGCAGGCCTGGc                                    >  2:159384/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGCGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCCCGCGATTCTGGCAGGCCTGGCATTcccc                             >  2:51677/1‑250 (MQ=255)
                                                                                                                                                                                                                          ggTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTg                        <  1:65747/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGACTGCGGTTAAACTGTCTTG‑TGGCTTCCTGTTCAAAATGTCGGGTCTGCCAGCTTCCGCAATTGCGATCTGGC‑ACTCTGCTAAACCAGAA‑AGCCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCGGCGGCGC‑TGACCT‑CCTTCCTGACCTGTATCACTGAGCCGATCGAG‑TTCTCCTTCAGGTTCGTTGCGCCTATCCTGTACATCATCCTCGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttctt                 <  1:247418/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttctt                 <  1:130120/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAACGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTCttcttc                >  2:528136/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAC‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCNGGCCTGGCATTCCCAATCTGTAttactctgg             >  2:260899/1‑244 (MQ=255)
                                                                                                                                                                                                                                      gCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTg              <  1:213120/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACGCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTgggg           <  1:57745/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aaCTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGt     >  1:175333/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tGTCTGG‑TAGCTTCCTGTTCAAAATGTACGGTCTGACAGCTGCCGCACTTGCTATCTGGC‑ACTCTGATAAACCAGAA‑AACCGC‑GCGAAAGCGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACAT‑CGTTCCTGACCGGTATCACCGAGCCGATCGCG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTTTATTCTTCTGGGGTTGTGTGAc  >  2:398263/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
CATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTA‑CCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGG‑TGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGC‑ACTCTGCTAAACCAGAA‑AACCGC‑GCGAAAGTGGGCGG‑TATTATGATCTCCGCGGCGC‑TGACCT‑CGTTCCTGACCGGTATCACCGAGCCGATCGAG‑TTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGAC  >  NC_000913/1158378‑1158867

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑