Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,7050GT25.4% 102.4 / 10.6 63intergenic (+38/‑126)ymiC/acnAprotein YmiC/aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (15/32);  new base T (16/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 8.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCT  >  NC_000913/1335462‑1335943
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
tcatccTTCATATTAACCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAATAGCTTGAGTGATCTCAACACAAAATAATAACCGTTAAGGGTGCCGCCTATGATCCACACAAATATGAGATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGTACGCCTCTGTACCATTACTCACGAATCCTTGCCACCGCGGCGGCGGCG‑TTGAGATTTGCATTGGcacctg                                                                                                                                                                                                                                             >  2:542802/1‑248 (MQ=255)
       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                                                                                                                                                                               <  2:330613/251‑1 (MQ=255)
                                                            aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                               <  1:401177/251‑1 (MQ=255)
                                                            aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                               <  1:534523/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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGG‑TTGTCATTTGCTTTGCCACAAGG‑TTTCTCCTCTTTTATCAATTTGGGTTGTTAT‑CAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAA‑TCAGGGCTT‑CGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCT  >  NC_000913/1335462‑1335943

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