Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,7170CT22.6% 117.1 / 10.4 62intergenic (+50/‑114)ymiC/acnAprotein YmiC/aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base C (17/31);  new base T (14/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.82e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGT  >  NC_000913/1335469‑1335961
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
tCTTAGTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTGAAGGGTGTAGTCTATGATCAACACAAATATGAAATATGGGTCCTGGATGGGCGCTTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTCCTTTGCCACAAGT‑TTtct                                                                                                                                                                                                                                                        <  2:448206/250‑1 (MQ=255)
         ccATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTTCTTTGCCACAAGG‑TTTTTCCTCTTTTAt                                                                                                                                                                                                                                              >  1:312165/1‑251 (MQ=255)
              tttAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          <  1:197424/250‑1 (MQ=255)
               ttAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTTG‑CATTTGCTTTGCCACAAGG‑GTTCTCCTCTTTTATCCAttt                                                                                                                                                                                                                                        >  1:479879/1‑251 (MQ=255)
                tAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGGCATTTCGTTTTCCTCAAGG‑TTTCTCCTCTTTTATCAATTTg                                                                                                                                                                                                                                       >  1:504263/1‑251 (MQ=255)
                  attaatttaatggCTTTTAATACACCGCAGCAAAAACAGCTTGAGTTATATCAACACAAAAGAATAAACGTTAAGGTTGTAGCCTATGATCAACACAAATAGGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGCCGATGCTCTTCTGGTACGATCTCCGCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTTTCATTTGCTTGGCCACAATG‑TTTCTCATCTTTTATCAATTTggg                                                                                                                                                                                                                                     <  2:314189/240‑1 (MQ=255)
                    gAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTTGATCATTACTCACTGATCCTTGGCCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGgtt                                                                                                                                                                                                                                   <  1:504562/251‑1 (MQ=255)
                          aaGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑c                                                                                                                                                                                                                            <  1:67590/251‑1 (MQ=255)
                                  ttAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTTCCACAAGG‑GTTCTCCTCTTTTATCAATTTGGGGTGTTAT‑CAAATCGtt                                                                                                                                                                                                                    >  1:52939/1‑251 (MQ=255)
                                                 aaTAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCtttt                                                                                                                                                                                                                                                <  2:344943/209‑1 (MQ=255)
                                                 aaTAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAAG‑GTTCTCCTCtttt                                                                                                                                                                                                                                                >  1:344943/1‑209 (MQ=255)
                                                  aTAACAGCTTGAGTTATCTTAAAACAAAATAATAAACGTTAACGGTGTAGCCTATGATCAACACAAATAGGAAATATTTGTCCTGGATGGGCGCGTTTTCTCTGTCGATTCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TGTCTCCTCTTTTATCAATTTGGGTTGGTAT‑CAAATCGTTACGCTATGTTTGTGtt                                                                                                                                                                                                    <  2:525308/251‑1 (MQ=255)
                                                   tAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTTCCACAAGG‑TTTTTCCTCTTTTATCAATTTTGGTTGTTTT‑CAAATCGTTTCGCGATGTTTGTGttt                                                                                                                                                                                                   >  1:330613/1‑250 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  1:534523/251‑1 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  2:330613/251‑1 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  1:401177/251‑1 (MQ=255)
                                                     aCAGCGTGAGTTATCACAACACAAAATAATAACCGTGAAGGGTGTAGCCTATGATCAACACAAATATTAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGTGCCGAACTCCTCTGGATCATTACTCACTGAGCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGTTTGTTAT‑CAAAGCGTTACGCGATGGTTGTGTTATc                                                                                                                                                                                                 <  2:145225/251‑1 (MQ=255)
                                                      cAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GCGTCATTTGTTTTTCCCCAAGG‑GTTTTCCTCTTTTATCACTTTGGGTTGTTTT‑CAAATCGTTACGCGATTTTTTTGTTATCt                                                                                                                                                                                                >  2:424114/1‑251 (MQ=255)
                                                      aatcTTGAGTTATATCAAAACAAAATAATAACCGTTTAGGGTGTAGCCTATGATCAACACAAATATAAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTGCTGGGCCGAACTCCTCGGGATCATTACTCACTGATCCTTGTCCCCGCTGCGGCGGTG‑TTGTCATTTGCTTTGCCACAAGG‑CTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCt                                                                                                                                                                                                <  2:336320/248‑1 (MQ=255)
                                                                       acacAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTTATTTTCTTTTTCCCAAGG‑GTTTTCCTCTTTTATCAATTTGTGTTGTTTT‑CACAACGTTACGCGNTGTTTGTGTTTTTTTTAattttttacctgt                                                                                                                                                                                >  2:257064/1‑239 (MQ=255)
                                                                              ataataACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTTTCCTCTTTTTTCAATTTGGGTTGTTAT‑CAAATCGTTTTGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑Tc                                                                                                                                                                       >  1:240790/1‑251 (MQ=255)
                                                                                 ataaCCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCTCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTGGAATATTCACCCTGAAGAGAA‑TCAgg                                                                                                                                                                    <  1:167207/251‑1 (MQ=255)
                                                                                 ataaCCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAgg                                                                                                                                                                    <  2:479879/251‑1 (MQ=255)
                                                                                       gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTTCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGGTGTTTT‑CAAAACGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑Cg                                                                                                                                                             >  1:344388/1‑251 (MQ=255)
                                                                                               tgtAGCCTATGATCAACACAAATATGAAGTATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGACCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCCTTTTCTTTTCCACAACG‑GTTCTTCTCTTTTTTCAATTTTGGCTGTTTT‑TAAAACGCGACTCAGTGTTTGCTTTTTCTTTAATATTTACCCTGAAGAGAAACCAGGG‑TT‑CGCCACCCTg                                                                                                                                                     >  2:312379/1‑251 (MQ=255)
                                                                                                  aGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCCTTTGCTTTGCCCCAAGG‑TTTCTCCTCTTTTTTCAATTTTGGGTGTTAT‑CAAATCGGTACGCGATGTTTTTGTTATCTTTAATATTTACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTcc                                                                                                                                                  >  1:161163/1‑250 (MQ=255)
                                                                                                       aTGATCAACACAACTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGGAATTTGTGTTGCCACAAAG‑TTTTTCCTCTTTTATCAATTTGGGTTGTTTTTCAAA‑CGTTACGCGGTGTTTTTTTTTTCTTTTATATTTACCCTGAAGAGAA‑TCAGGGATT‑TGCCACCCTGTCAttta                                                                                                                                              >  2:78194/1‑248 (MQ=255)
                                                                                                       aTGATCAACACAAATATGATATATTTGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTTTTCTGGGCCGAACTCCTCTGGATCATTACTCACTTATCCTTGACCCCTCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAATG‑TTTCTCCTCTTTTATCAATTTTGGTTGGTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCAttt                                                                                                                                               <  2:68412/249‑2 (MQ=255)
                                                                                                          aTCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTGCCACACAG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑TAAAACGTTTCGCGATGTTTGTGTTTTCTTTAATATTCACCCTGAAGAGAA‑TTAGGGCTT‑CGCAACCCTGTCATTAagg                                                                                                                                            >  1:499739/1‑249 (MQ=255)
                                                                                                           tCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTGCCACCAGG‑TTTCTCCTCCTTTTTCAATTTTGGTTGTTTT‑CAAATCGTTACGCCGTGTTTGTGTTATTTTTTATATTTACCCTGAAGAGAAATC‑GGGCTT‑CGCAACCCTGTCATTattgagg                                                                                                                                         >  2:501234/1‑245 (MQ=255)
                                                                                                            cAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACGCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAaggagg                                                                                                                                         <  1:47246/250‑1 (MQ=255)
                                                                                                               cacaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCATGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCt                                                                                                                                     <  1:199021/251‑1 (MQ=255)
                                                                                                               cacaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTTATTTGCTTTTCCACAAGG‑TTTTTCCTCTTTTATCAATTTTGGTTGTTTT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAAGGAA‑TCAGGGCTTTCCCAACC‑TGTCATTAAGGAGGAGCt                                                                                                                                     >  1:264057/1‑251 (MQ=255)
                                                                                                                   aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTTCTTTTCCACAAGG‑TTTCTCCTCTTTTTTCAATTTTGGGTGTTAT‑TAAATCTGTACGCGGTGTTTGTGTTATCTTTAATATTTACCCTGTAGAGAA‑TCAGGGCTT‑CGCAACCCTTTCATTAAGGAGGAGCTATgt                                                                                                                                 >  1:398930/1‑251 (MQ=255)
                                                                                                                      atGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATgtcgt                                                                                                                              <  1:4304/251‑1 (MQ=255)
                                                                                                                      atGAAATATTGGTACTGGATTGGCGCGTTTTCTCTGTCGATGCTCTTCTTGGCCGAACTCCTCTGTATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTTCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATgtcgt                                                                                                                              <  1:270244/251‑1 (MQ=255)
                                                                                                                           atatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTTCCACAAGGGTTTTTCC‑CTTTTATCAATTTGGGGTGTTAT‑CCAATCGTTACGCGATGTTTGTGTTTTCTTTTATATTTACCCTGAAGGGAA‑TCAGGGCGT‑CGCAACCCCGTCCTTACAGAGGCGCAATGTCGTCaaac                                                                                                                         >  1:387728/1‑249 (MQ=255)
                                                                                                                           atatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTGCCCCCAGG‑GTTTTCCTCTTTTTTCAATTTGGGGTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTTACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAAcc                                                                                                                         >  1:3898/1‑251 (MQ=255)
                                                                                                                            tatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCCTTTGCTTTTCCACAAGG‑TTTTTCCTCTTTTTTCAATTTGGGTTGGTAT‑CAAAACGTTACGCGATGTTTGTGTTTTTTTTAATATTCACCCCGTAGAGAA‑TCAGGGGTT‑TGCAACCCTTTCATTTAGGAGGAGCTATGTCGTCAAccc                                                                                                                        >  2:526164/1‑251 (MQ=255)
                                                                                                                                ggTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTTCCACAAAG‑TTTTTCCTCTTTTATCAATTTGGGGTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTTTTTTTAATATTCAACCTGAAAGGAA‑TCAGGGCTT‑CCCCACCCTGTCATTTAGGGGGGGGTTTGTCGTCaaacccacg                                                                                                                    >  1:308670/1‑244 (MQ=255)
                                                                                                                                      ggATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAg                                                                                                               <  1:526164/250‑1 (MQ=255)
                                                                                                                                          gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          <  2:277217/126‑1 (MQ=255)
                                                                                                                                          gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCCTTTTCTTTGCCACAAGG‑GTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          >  1:277217/1‑126 (MQ=255)
                                                                                                                                                gTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTTTCATTTGCCTTTCCACCAAG‑TTTCTCCTCTTTTATCAATTTGGGTGTTTAT‑CAAATCGTTACGCGCTGTTTGTGTTATCTTTAATATTTACCCTGACGAGAA‑TCAGGGCTT‑CGCACCCCCGGCATTAAGGAGGAGCTATGTCGGCAAccccccgaagagccagcaaggcc                                                                                                    >  2:27988/1‑231 (MQ=255)
                                                                                                                                                gTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTTCCACCAGG‑TTTCTCCTCTTTTATCAAATTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACAAGAAGCCAGTCAGGa                                                                                                     >  1:315001/1‑250 (MQ=255)
                                                                                                                                                         gTCGATGCTCTTCTGGGCCGAACTCCTCTGGAACATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTGGGGTTGTTAT‑CAAATCGTTACGCGATGGTTGTGTTATCTGTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCTTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAg                                                                                           <  1:193400/251‑1 (MQ=255)
                                                                                                                                                            gATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAACTTGGGGTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCCGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAAGACAACTTGCAGGcc                                                                                        >  1:76434/1‑251 (MQ=255)
                                                                                                                                                               gCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCGTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTAGCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGGGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCaaa                                                                                     <  1:78194/251‑1 (MQ=255)
                                                                                                                                                                 tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTa                                                                                                                                                                                                   <  1:485543/141‑1 (MQ=255)
                                                                                                                                                                 tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCCTTTTCTTTTCCACCAGG‑TTTTTCCTCTTTTTTTAATTTTGGGTGTTAT‑TAAATCGTTACGCGGTGTTTTTGttt                                                                                                                                                                                                   >  2:485543/1‑140 (MQ=255)
                                                                                                                                                                  cttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAg                                                                                    <  1:225923/249‑1 (MQ=255)
                                                                                                                                                                    tctGGGCCGAACTCCTCTGGATCATTACCCACTTATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTAACAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATa                                                                                 <  1:363974/250‑1 (MQ=255)
                                                                                                                                                                            gAACTCCTCTGGATCATTACTCACTGATCCTNGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTAcc                                                                        <  2:264057/251‑1 (MQ=255)
                                                                                                                                                                                      ggATCATTACTCACTGGTCCTTGACCCCNCTGCGGCTGGG‑TTGTCAGTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGGTTGTGTTATCGTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCActacta                                                                 <  2:266829/248‑1 (MQ=255)
                                                                                                                                                                                       gATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCCTTTGCTTTTCCACAAGG‑TTTTTCCTCTTTTTTAAATTTTGGATGTTTT‑CAAATCGTTACGCCATGATTTTGTTATTTTTTATATTCTCCCCGAAGGGAA‑TCAGGGCAT‑CGCAACCCTGTCATTAAGGAGGAGCTATTTCGTCAAACCTAAGTGAAGCCAGGAAGGGCACGATGCAGGACAACGATAAAACTTACCACTACTAAAGc                                                             >  1:405173/1‑251 (MQ=255)
                                                                                                                                                                                        aTCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑GTGTCATTTTCTTTTCCACAAGG‑GTTTTCCTCTTTTATCAATTTGGGGTGTTTT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTTCCCCTGAAGAGAA‑TCAGGGCTT‑TGCAACCCTGCCATTTAGGAGGAGCTATGTAGTCAAACCTACGAGAAGCCAGTAAGGACAAGTTGCAGGCCCAAGATAAAACTTAACACCACTACAGcc                                                            >  1:101903/1‑251 (MQ=255)
                                                                                                                                                                                                 cACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGGTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCGT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACGTACCACTACTACAGCCTGCCGCTtg                                                   <  1:173857/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGCCCCCGCTGCGGCGGGG‑TTGTCATTTGTTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGATATGTCGTCAACCCTACGAGAAGCCAGTAAGGATACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgct                                                 <  1:484085/249‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCGGCCGCTtgctgc                                               <  1:110560/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGACCCCGCTGCGGCGGGG‑GTTTCCTTTTCTTTGCCCCAAGG‑TTTCTCCTCTTTTTTCAATTTTGGTTGTTAT‑CAAAATGTTACGCGCTGTTTGTGTTTTCTTTAATATTCTCCCCGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACACTACGAGAAACCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAACCTGCCGCTtgctgc                                               >  1:383213/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCt                    >  1:450917/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tCATTTGCTTTGCCACAAGG‑TTTCTCCTCTNTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCt                    <  2:240790/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGCCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGt  >  2:497131/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
TCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGT  >  NC_000913/1335469‑1335961

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