Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,608	0	G	A	50.0%	8.5 / 12.5	20	E142K (GAA→AAA)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base A (8/2);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGT  >  NC_000913/2151467‑2151755                                                                                                                                               |                                                                                                                                                    gACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATTTCATTGTTTTTTTTTTTGTTGACTATA‑CAAaaaaaga                                                                                                                                                 >  2:649154/1‑146 (MQ=255) gACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCcg                                                                                                                                            <  1:25567/151‑1 (MQ=255)    ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTTATTCTTTTTTTGTTACCTCTA‑CAAaaaaagaa                                                                                                                                                >  2:174271/1‑144 (MQ=255)                 cAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCaaa                                                                                                                            <  1:402971/151‑1 (MQ=255)                    aCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAtt                                                                                                                         <  2:203077/151‑1 (MQ=255)                        aGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCCTTTATTCTTTTTTTGTCACCCCTC‑CCAAAAAAAAAACCGACGAGCTAAAAGGAAAATTAGt                                                                                                                      >  1:1826/1‑150 (MQ=255)                          aaTATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTAACCTTTTTTGGGCACCCTTT‑CTAAAAAAAAAGACCACGAGCTAAGAGCaaaa                                                                                                                           >  1:65024/1‑143 (MQ=255)                                       gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTTTCCTTTTTTTGTTACCTATCTCAAAAAAAAAAGCAG‑CGAGGTAAGAGGAAAAATATTTTTCTTATCaagc                                                                                                         >  2:120044/1‑145 (MQ=255)                                       gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTTATCCTTTATTTGTCACCTATA‑CAAAAAAAAAAACCGACGAGCTAAGAGAAAAAATAATTTGCTAATCAgac                                                                                                         >  1:15042/1‑146 (MQ=255)                                                     aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTTATCCTTTTTTTGGTTCCCATC‑CCAAAAAAAAAAGCGACGAGGTAAAAGGAAAAATAGTTTGCTAATTAACCGGTATAttttttttttt                                                                                        >  1:6788/1‑150 (MQ=255)                                                       gaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCTTTATCCTTTTTTTGTTCCCTAAC‑CAAAAAAAAAAAACGCCAAGCTAAAAGGAAAAATAATTCGCTTAACAACCGTTAtttttttttttt                                                                                         >  2:254754/1‑148 (MQ=255)                                                             atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGc                                                                                 <  1:56382/150‑1 (MQ=255)                                                              tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCATTTTTCCTTTTTTTGTCACCTATC‑CAAAAAAAAAAAACCACGAGCTAAGAGCAAAAATAATTTGCTAATTAGCCGTTATATTTATTTTTTTAAATTGCgg                                                                               >  1:116974/1‑151 (MQ=255)                                                                         ctGTATACGCGCGGGGATAAAATTTACGGCGGTTATTCCATTTATTCTTTTTTTGTCACCTATC‑CAAAAAAAAAAAACGACCAGCCAAGAGAAAAAATAATTTGCTAATAAGCCGTTATATTTATTTTTTTAAATAGCgg                                                                               >  2:684878/2‑140 (MQ=255)                                                                                                          tataCCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAAATCACTCTTccc                                    >  1:85841/1‑150 (MQ=255)                                                                                                             aCCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACCGGCGGTTGGATTGTTTTCTTCAACCGCCCTGAGAAAATCACTCTGccc                                    >  1:142770/1‑147 (MQ=255)                                                                                                                 ttGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGcc                            >  1:92791/1‑151 (MQ=255)                                                                                                                    accttttATTGGTCACTTACC‑AAAAAAAAGAACCCGACGGGTTAAGAGAAAAATTATTTCGATAACCACCCTTTTTTTTTATTTTTAAACCGGGGGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                           <  2:640132/145‑1 (MQ=255)                                                                                                                    aTCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTCTTCAACCCACCTGAGAAAATCACCCTGCCCACAGCGCCAGt                         >  2:334850/1‑151 (MQ=255)                                                                                                                          tATTGGTCACCTATC‑AAAAAAAAGAAGCCGACGAGCTAAGAACAAAATTAGTTCGCTAATCAGCCGTTATTTTTTTTTTTTAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                    <  2:248619/150‑1 (MQ=255)                                                                                                                                 cACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCCCCAGTGCCCACGTTTTat            >  2:65190/1‑151 (MQ=255)                                                                                                                                 cACCTACC‑CAAAAGAAGAGGCCGAAGACGAAAGAGAAAAATAAGTTCGCTAATCAGCCTTTATATTTTTTGTTTAAGCTGGGGGTTAATTTGTTGTCTTCATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat            <  1:120046/151‑1 (MQ=255)                                                                                                                                             aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGt  >  1:295116/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    GACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATC‑CAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGT  >  NC_000913/2151467‑2151755

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