Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	691,883	0	C	T	28.6%	82.2 / 12.6	49	E153E (GAG→GAA)	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base C (11/24);  new base T (14/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 9.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAAC  >  NC_000913/691637‑692129                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                              tAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGTCAGTCCAA‑AGCCCGGC‑AAATTATTCTTTTTCgg                                                                                                                                                                                                                                                           >  2:212447/1‑251 (MQ=255)        gTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCACTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGTC‑GCCCGTCAAA‑AGCCTTGC‑AAAAAATTTTTTTTCGGAAAT‑Gt                                                                                                                                                                                                                                                    >  1:390665/1‑250 (MQ=255)              cccGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATc                                                                                                                                                                                                                                             <  2:465/251‑1 (MQ=255)                cGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGCCAGTCAAA‑TTCCCTGC‑AAAATATTCTTTTTCGGCAAT‑GTACGGATccc                                                                                                                                                                                                                                           >  1:134635/1‑250 (MQ=255)                 ggAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTc                                                                                                                                                                                                                                          <  1:317800/251‑1 (MQ=255)                         ggATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TGCGCCAGCCCAA‑AGGCTGGC‑AAATTATTTTTTTTTGGCAAT‑GTACGGAGCCACC‑TAGcccc                                                                                                                                                                                                                                 >  1:422894/1‑250 (MQ=255)                                           cGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑CAATTATTTTTTTTCGGCAAT‑GTACGGAACCTCA‑TAGCCCCGAGCAAACAAAATctctt                                                                                                                                                                                                               >  1:74687/1‑250 (MQ=255)                                              cACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGttt                                                                                                                                                                                                            <  2:301689/251‑1 (MQ=255)                                                     tttCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAggg                                                                                                                                                                                                     <  1:328707/250‑1 (MQ=255)                                                       tCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑Ct                                                                                                                                                                                                  <  2:65462/250‑1 (MQ=255)                                                          gTTCACCGTGGAAAAGTTGGATGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGAGG‑TCCTCCAGATAAA‑CGCCGGGC‑ACATTATTCTTTTTCGGCAAC‑ATACGGCTCATCT‑AATCACAGAGCAAGCCAAATATCTGTTT‑Cgtgcgccctcc                                                                                                                                                                                               >  2:146835/1‑241 (MQ=255)                                                                   ggAAAAGTTGGCTGGGTAACAGGGAGATAATTCCCTTCTTGTTTCTTATCGTGTCACTACTGTGTGAATTTTCGTCTCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATTAAAACGCTGTCGTTAATCTCCATCTGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGTCTAT‑GTACGGATCCTCA‑TATCCCAGAGCAAGCATAATCTCTGTTC‑CAAGGG‑CTTCC‑ATttcttc                                                                                                                                                                                      <  2:185405/249‑1 (MQ=255)                                                                        aGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGCCCGGCCAA‑AGCCTGGC‑CAATTATTTTTTTTCGGCAAT‑GTACGGATCCTCA‑TAACCCAACGCAAGCATAATTTCTTTTT‑CGAGGG‑CTTTC‑ATttcttttgtctcc                                                                                                                                                                               >  2:125300/1‑243 (MQ=255)                                                                            ggCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCGCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑GTTCC‑ATTTCTTCTGCTTCGTCa                                                                                                                                                                            <  1:64248/250‑1 (MQ=255)                                                                             gCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATc                                                                                                                                                                          <  1:298421/251‑1 (MQ=255)                                                                             gCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATc                                                                                                                                                                          <  1:1674/251‑1 (MQ=255)                                                                                  gTAACAGGGAGAACAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGACGCTAATGGCCGCGTATGGGTTCTCATGTTCGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCACTCAAC‑CACCTGGC‑AAAATACTCTTTCTCGGCAATCGAACG‑ATTCTCC‑AACCCCAGAGCACGCATAATTTCTGTTT‑TGAGGG‑CTTTC‑ATTTTTTTTGCTTccgccatctcg                                                                                                                                                                      >  1:384670/1‑240 (MQ=255)                                                                                    aaCAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAAATCCCTG‑C‑AAAATATTCTTTTTCCGCAAT‑GTACCGAACCTCA‑TAGCCCAGAGCAAACATAATCTCTGTTT‑CGCGGG‑CTTCCCATTT‑TTCTTCTTCGTCATCTTCCAtg                                                                                                                                                                   >  2:416589/1‑251 (MQ=255)                                                                                     aCAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑CAATTATTCTTTTTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTTC‑ATTTCTTCTGCTTCGTCATCTTCGAttt                                                                                                                                                                  >  1:111122/1‑249 (MQ=255)                                                                                            gaAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGACGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTTGC‑AAATTATTATTTCTCGGCAAT‑GTAAAGGTCCTCA‑TAGCCCAGAGCAAGCATAATCTCTTTTT‑CTCGGG‑GTTCC‑AATTCTTCTTCTTCGTCCTCTTCGCTGTGATCGt                                                                                                                                                            >  1:336900/1‑250 (MQ=255)                                                                                               aaaTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACC‑GCCGCCAGGCAAA‑TACCGGGC‑AAAATATTTTTTTCCGGCAAT‑GTACGGAGCCCCA‑TAGACCAGAGTAAGCATAAACTCTGGTT‑CGAGGGGCTCCC‑TTTTCT‑CTGCTTTCTTATCTTCGCTGTGATCGTAAc                                                                                                                                                         >  1:425448/1‑250 (MQ=255)                                                                                                aaTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCt                                                                                                                                                       <  1:423113/251‑1 (MQ=255)                                                                                                          tgttgCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TATGCCAAGCAAA‑TGCCTGGC‑AAATTATTCTTTTTCGGCAAT‑GTACGGGTTCTCA‑TAACCCAGGGCAAGCATAAACTCTGTTTTCG‑GGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGTGCGGAACCTCACa‑aaaga                                                                                                                                             >  1:154717/1‑247 (MQ=255)                                                                                                              gCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGAc                                                                                                                                         <  2:349233/250‑1 (MQ=255)                                                                                                                    cGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGt                                                                                                                                   <  2:4901/250‑1 (MQ=255)                                                                                                                          actactGTTTGAATTGTCGTCGCGCATGGCGTCGTATGGGTTCTCATGTTAGGTAATCAAACCGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGTTCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGcaccac                                                                                                                            <  2:260464/251‑1 (MQ=255)                                                                                                                             actGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCaccacc                                                                                                                           <  1:67824/249‑1 (MQ=255)                                                                                                                                    gAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑CTTGCCCGGCAAA‑AGCCCGGCCAAATT‑TTCTTTTTCGGCAAA‑GTACGGGTCCTCA‑TAACCCCGAGGAAGCAAAATCTCTTTTT‑CGGGGG‑CTTCC‑CTTTCTTCTGCTTCTTCCTCTTCGTTGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCCCCCTATgcccc                                                                                                                  >  2:54928/1‑248 (MQ=255)                                                                                                                                            gtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑CTTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑CGTgcgc                                                                                                         >  1:221979/1‑251 (MQ=255)                                                                                                                                                  cATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAg                                                                                                   <  1:36245/251‑1 (MQ=255)                                                                                                                                                       cgtcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGCCCGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGGTCCTCA‑TAGCCCAGAGCAAACATAATCTCTGTTT‑TGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAAACTAACA‑AAAGCAGACTGCCCTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGtt                                                                                              >  1:134011/1‑251 (MQ=255)                                                                                                                                                          cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCGTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGGTCGCCTCCAGTGGTTTg                                                                                            <  2:147966/250‑1 (MQ=255)                                                                                                                                                           gtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑CCTGCCAGTACAA‑TGCCCTGC‑AAAATATTTTTTTTCCGAAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑GTTCC‑ATTTCTTCTGCTTTGTCATCTTCGATGTGATCGTAAACTAACA‑AATACAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCACCTCCAGCGGGTTGcc                                                                                          >  1:338673/1‑251 (MQ=255)                                                                                                                                                                       aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTTCCAGTCAAA‑ATCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CCTCC‑CTTTCTTCTGCTTCGTCATCTTCCATGTGATCGTAACCTAACA‑AATGCAGACTTCCGTGCACCACCATATGCGCCCCAG‑GCGCCTCCCATGGTTTGCCCTGCTCCTgggc                                                                              >  1:266865/1‑248 (MQ=255)                                                                                                                                                                          ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGAAG‑TTTGTCAGACAAA‑TGACTTGC‑AAAATATTCTTTTTCGGGAAT‑GTTCGGATACTCAATA‑CCCACAGCAAACAAAATTTCTTTTT‑CGTGGG‑CTTCC‑CGTTTTTCTTCTTCGGCATCTTCGCTGTGATCTGAACCCAACACAAT‑CAGACTACCGTGCGCCCCAACATGAGCCC‑CGCGAGCCTCCAGTGGTTTTCCCTGCTCCTGAGcttc                                                                           >  2:58100/1‑251 (MQ=255)                                                                                                                                                                            agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGcttcct                                                                         <  1:437493/251‑1 (MQ=255)                                                                                                                                                                             gttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCG‑TCTGCCAATCAAA‑TGCCTGGC‑AAAATATTCTTTCTCGGCAAT‑GTAAGGAGCCTCA‑TAACCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCCTCGATGTGATCGGAACCTAACA‑AATGCAGACTGCCGTGCACCACCCTATGCGCCC‑AGTGCGCCTCCAGTGGGTTGCCTTGCTCCTGAGcttcctt                                                                        >  1:185087/1‑251 (MQ=255)                                                                                                                                                                                aaTCAAAGCGCAGTCGATAATCACCAACGGCGGGGACG‑CCTGCCAGACAAA‑TTTCTGGC‑AAATTATGCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCGCTGTTT‑CGAGGG‑CTTCC‑ATTGCTTCAGCTTCGTCATCTTCGATGTGATCGTATCCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGACTTGCTCCTGAGcttcctgct                                                                      <  2:424266/250‑4 (MQ=255)                                                                                                                                                                                       aCGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑TTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGCGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCCCCAGTGGGTTGCCTTGCTCCTGAGCTTCCTTCTCAaccacc                                                              >  1:238662/1‑251 (MQ=255)                                                                                                                                                                                                     cACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATaa                                                 <  1:314368/250‑1 (MQ=255)                                                                                                                                                                                                                 ggACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGATCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGTTCGCCCAg                                    <  1:230788/251‑1 (MQ=255)                                                                                                                                                                                                                   aCG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGt                                   >  2:54726/1‑250 (MQ=255)                                                                                                                                                                                                                     g‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGc                                <  1:416097/251‑1 (MQ=255)                                                                                                                                                                                                                           ccAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACAt                           <  1:339519/251‑1 (MQ=255)                                                                                                                                                                                                                              gTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTc                        <  2:221979/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCAAggcg             >  2:153808/1‑251 (MQ=255)                                                                                                                                                                                                                                           gC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAggcgg            >  1:146792/1‑251 (MQ=255)                                                                                                                                                                                                                                            c‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAggcggc           >  1:221358/1‑251 (MQ=255)                                                                                                                                                                                                                                                   attCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTc      <  1:398182/250‑1 (MQ=255)                                                                                                                                                                                                                                                      ctttctCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAAc  >  1:245002/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                              TAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACG‑TCTGCCAGTCAAA‑TGCCTGGC‑AAATTATTCTTTCTCGGCAAT‑GTACGGATCCTCA‑TAGCCCAGAGCAAGCATAATCTCTGTTT‑CGAGGG‑CTTCC‑ATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACA‑AATGCAGACTGCCGTGCACCACCATATGCGCCC‑AGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAAC  >  NC_000913/691637‑692129

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