Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	862,559	0	G	T	39.4%	31.2 / 13.2	33	H16Q (CAC→CAA)	ybiW	putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (5/15);  new base T (13/0);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 2.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTTCGATCCACGAGACAGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAAT  >  NC_000913/862316‑862809                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                tttttCGATCCACGAGACAGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAA‑A‑TAGCACC‑AGCGCA‑‑TTTTTTTGGGCtt                                                                                                                                                                                                                                                          >  1:297407/1‑251 (MQ=255)  tttgtgATCCACTAGACAGTATATTCCGTGAAGAGCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTGGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCTCGACGTACCGGGATCGGCTTATCGAGATTTTTTTGATACATCTCGTTATAGTGCTGCGCGCGCTCGGGACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑ACCGCA‑‑TTTTTGTgcgatt                                                                                                                                                                                                                                                          <  2:243031/245‑4 (MQ=255)          ccACGAGACAGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCgg                                                                                                                                                                                                                                                  <  1:217488/250‑1 (MQ=255)            aCGAGACAGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGAGTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑ATTTTGTGCGCTTTAATGCGGTCg                                                                                                                                                                                                                                               <  2:244595/251‑1 (MQ=255)                   aGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCg                                                                                                                                                                                                                                        <  1:209144/251‑1 (MQ=255)                         ttCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GGACC‑CGCGCA‑‑TTTTTTTGTGCTTTTATGGGGtgtcctagcgtggccc                                                                                                                                                                                                                                  >  1:283425/1‑236 (MQ=255)                            cGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TACAATA‑T‑TCACCCAG‑GCA‑‑TTTTTTTGAGGTTTATAGCGGGCGCTGAGCGGG‑TCCCGtt                                                                                                                                                                                                                              >  1:272680/1‑251 (MQ=255)                              ggAAGAGCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTc                                                                                                                                                                                                                            <  2:229067/251‑1 (MQ=255)                                 aGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCAAC‑CGCGCC‑‑TTTTTTTTTGCTTTTATGCGGTCGCTGAGCGGG‑TCCCGTTTctgt                                                                                                                                                                                                                         >  1:28115/1‑248 (MQ=255)                                          cggcgcgAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCCGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TCCAAAA‑T‑TCACC‑CGCGCA‑‑TTTTTTTGCGCTTTTATGCGGTCGCTGGGCGGG‑TCCCGTTTCC‑GTGTGGT‑cttt                                                                                                                                                                                                              >  1:51289/1‑248 (MQ=255)                                                            ccTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAgg                                                                                                                                                                                            <  1:331935/251‑1 (MQ=255)                                                               ggTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTAGCCTCGTAAGGTc                                                                                                                                                                                          <  1:120868/250‑1 (MQ=255)                                                                              aCTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TACAATA‑T‑GCACC‑CGCGCC‑‑TTTTTTTGCGCTTTAATGCGGGCGCTGAGCGGG‑TCCAGTTTTA‑GTGTGGG‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑GTAACCCTTTCt                                                                                                                                                                         >  1:5368/1‑251 (MQ=255)                                                                               cTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAAAA‑T‑GCACC‑AGCGCC‑‑TTTTTTTGCGGTTTTATGGGGGCGCTGAGCGTG‑TCCAGTTTTT‑GGGTGGT‑CATTGGTGTTTTCCTCGTAAGGTCGCGG‑TTAAACCTTTCTg                                                                                                                                                                        >  1:287847/1‑251 (MQ=255)                                                                                    cGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTggcag                                                                                                                                                                    <  2:5368/250‑1 (MQ=255)                                                                                     gtgtTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CATAATAAT‑GC‑CC‑CCCGCC‑‑CTTTTTTGCGCTTTTATGTGGGCGCTGGGCGGG‑TCCAGTTTTATGGGGGGT‑C‑TGGTTGTTTTCCTCGTGATGTCGCGG‑TTTAcccctttttgccggg                                                                                                                                                                  >  2:152975/1‑239 (MQ=255)                                                                                                         gATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TACAATAAT‑G‑AAC‑ACCGCA‑‑TTTTTTTGGGCTTTTATGCGGGTGCCGGGCGGG‑GCCAGATTTA‑TGGTGGT‑CATGGTTGTTATTCTCGTAAAGTCGCGG‑GTTACCCATTTTGGCCGGCCTACTTCTTGGGCaagtacg                                                                                                                                              >  1:326060/1‑246 (MQ=255)                                                                                                                cAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTTTAACAAAA‑T‑‑TACC‑ACCGCA‑‑TTTTTTTGCGCTTTTATGCGGTCGCCTAGCGGGGTCCACTTT‑A‑GCGTGGG‑GCTGGTTGTTTTCATCGTGAGGTGGCGGGTTAACC‑TTTCTGGCCGGCCTACTGGC‑GGGCAAAAACGAGCatttt                                                                                                                                      >  2:335832/1‑247 (MQ=255)                                                                                                                      aTGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TCCCATA‑T‑GCACC‑CACGCC‑‑TTTTTTTTCGCTTTTATGCGGGCGCCGGGCGTG‑TCCCGGTTCT‑GTGGGGTGCTTGGTT‑TTATTCTCGTAAGGGCGCGG‑GTATCCCTTTCTGGGAGGCATAATGCTTGGGCaaatcgcgcaagtttggcgg                                                                                                                                 >  1:367620/1‑234 (MQ=255)                                                                                                                                            cGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TACAATA‑T‑GCCCC‑CACGCA‑‑CTTTTGTGCGCTTTAATGCGGGCCCTTAGCGTG‑TCCAGTTTCC‑GTTGGGG‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑GTAAACCTTTCTGGGAGGCATACTGCT‑GGGCAAAGTCTAGCAATTTTGGCGCATCCAGCGGGTTTTTCGGGGCg                                                                                                          >  1:235543/1‑251 (MQ=255)                                                                                                                                               aCCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑ACACC‑AGCGCA‑‑TTTTTTTGCGCTTTTATGCGGGCGCTGGGCGCG‑TCCCGGTTTA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGACGAGCCGTTCTGGCGCATCAAGCGGTTTTTCCCGGGCGTc                                                                                                        >  1:344850/1‑250 (MQ=255)                                                                                                                                                                        gttgATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGatat                                                                               <  2:80815/250‑1 (MQ=255)                                                                                                                                                                             tACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑TTCAAAA‑C‑GCAAC‑ACAGCAATTTTTTGTGCGCT‑‑AATGGGGGCGCAGCGCGGG‑TCCCGTTTTA‑GTGGGGG‑GATTGTTGTTATTCCCGTAAAGGGGCGG‑TTTACCCCTTCTTGGAGGCGTAATGCT‑GGGGAAAGAAGAGCAGTTTTGGGGCCTCAAGCGGGTTTTTCGGGGCGTCGTAGGGCCGATtaagtgaggggtgttttttg                                                                         >  1:152631/1‑236 (MQ=255)                                                                                                                                                                                            tgctgcGCGCGCTCGGTACAGACTGGCGGTTT‑CACAAAA‑T‑TCACC‑AGCGCC‑‑TTTTTTTGCGCTTTAATTCGGGCGCTGAGCGGG‑GCCAGTTTCC‑GTGTGGT‑CATTGTTGTTATTCTCGTAAGGGCGCGG‑GTAACCCTTTTTGGCAGGCATTCTGCT‑GGGCGAAATCGAGCCGTTCTGGCGCATCAAACGGGTTTTTCTGGGCGGCCTAAGGCAGATTTAGGAAGGGATATTTTTTGATGCCCCGCGTGTgg                                                          >  2:9135/1‑251 (MQ=255)                                                                                                                                                                                                         cgtgaCAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGTTCGCTGAGCTTG‑TCCAGTTTCA‑GTTTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGGGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAAt                                             <  1:43555/247‑1 (MQ=255)                                                                                                                                                                                                                       ggTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAg                                                                                               >  1:183230/1‑187 (MQ=255)                                                                                                                                                                                                                       ggTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTGTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAg                                                                                               <  2:183230/187‑1 (MQ=255)                                                                                                                                                                                                                          tt‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGc                            >  1:95806/1‑251 (MQ=255)                                                                                                                                                                                                                                aaTA‑T‑GCCCC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCCCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGAGTAAGTAAGTGATATTTGTTTATGCCCAGCGTGTGGTAGTGCAGAAAAGGAATTTCTCCAACGTGCAGCtc                       <  2:234497/251‑1 (MQ=255)                                                                                                                                                                                                                                   a‑t‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCGGCtcgtc                    >  2:296472/1‑251 (MQ=255)                                                                                                                                                                                                                                   a‑t‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCtcgtc                    <  1:332518/251‑1 (MQ=255)                                                                                                                                                                                                                                           c‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCggcggc              <  2:95806/251‑1 (MQ=255)                                                                                                                                                                                                                                                gcA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCaaaa          <  1:9135/251‑1 (MQ=255)                                                                                                                                                                                                                                                          gtgCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAAt  >  1:313452/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                TTTTTCGATCCACGAGACAGTATATTCCGGGAAGATCGGCGCGGCGCGAACTTCGCTTGCCTGGTTGCCAATGATCAACTCATCGTGTTTGATCCAGATGGTGCGATTCGCCAGGTGATGCGCCAGTGCCAGCGCGCGACGTACCGGGATCGGCTTATCGAGATGTTGTTGATACATCTCGGTATAGTGCTGCGCGCGCTCGGTACAGACTGGCGGTTT‑CACAATA‑T‑GCACC‑AGCGCA‑‑TTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTG‑TCCAGTTTCA‑GTGTGGT‑CATGGTTGTTATCCTCGTAAGGTCGCGG‑TTAACCCTTTCTGGCAGGCATACTGCT‑GGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAAT  >  NC_000913/862316‑862809

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