Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,423,272	0	C	T	60.0%	‑1.4 / 11.3	15	S174L (TCG→TTG) ‡	ydaW	putative uncharacterized protein YdaW
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (3/6);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGAGCTGGCTACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCG  >  NC_000913/1423040‑1423501                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      gatgaGCTGGCTACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAgcgc                                                                                                                                                                                                                     <  1:37754/251‑1 (MQ=255)                                                                                                                     cATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGAATCGGTATCTCAGCTATGGTTATCTGCATCGTCGTCTGTCTCTCATGGGGTGTTCCTCATTACCGTGATAACGCCATCGCCTCCACAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGc                                                                                                                                                                                        >  1:159021/1‑164 (MQ=37)                                                                                                                        aCCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTCCCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAATAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGTTGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATTCAAAAGACACGAAGGAGTTAGCTGAGGCGAAAg                                                                                               <  1:54888/250‑1 (MQ=255)                                                                                                                                                      agcagcATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTTATGTTGCTGCGCTCGAGGCAAAATACACGAAGGAGTTAGCGGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGttcc                                                                <  1:160278/251‑3 (MQ=255)                                                                                                                                                          gcaTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGcn                                                            <  1:103186/251‑2 (MQ=255)                                                                                                                                                                 tNNCTCCGCNNTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTCCCGTGATAACGCCACTACCTACAAAGCCCAGCGCGACaaaaa                                                                                                                                                                                                             <  2:8790/99‑1 (MQ=14)                                                                                                                                                                    ctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTc                                                                                                                                    <  1:51657/169‑1 (MQ=11)                                                                                                                                                                      ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACTCGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcgtc                                                  <  1:54199/249‑1 (MQ=11)                                                                                                                                                                                          atcgtnNGCNNNTCNNGGGCTGTTAATCATTACCGTGATNACGCCATTACCTACAAAGCCCAGCGCGACNANANTGCCNNAGANCTGAAGCTGGNGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCNCNCGATGCAACATACACGCCGGAGTTAGCTGAGNCTAAAGCTGACCATGATGCTCTGcnggagnnggnt                                                                  >  2:127479/2‑202 (MQ=11)                                                                                                                                                                                          atcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTgg                                                          <  1:127479/220‑1 (MQ=11)                                                                                                                                                                                            cgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGCTGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtc                          >  1:11756/1‑251 (MQ=11)                                                                                                                                                                                                   cTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGAGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGCTGCGAAAGCTGAAAATGATGCTCTGCGtg                                                                          >  1:161316/1‑196 (MQ=255)                                                                                                                                                                                                     ntCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATATTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGc                 >  1:30814/2‑251 (MQ=11)                                                                                                                                                                                                            gctgtTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGAGGCTCTGCTTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGc          <  1:145835/249‑1 (MQ=11)                                                                                                                                                                                                                     cATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCTTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCg  <  1:124050/250‑1 (MQ=11)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      GATGAGCTGGCTACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCG  >  NC_000913/1423040‑1423501

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