Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	40.7%	20.4 / 14.7	27	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (13/3);  new base A (0/11);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 3.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTG  >  NC_000913/1222456‑1222737                                                                                                                                             |                                                                                                                                              gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaacaaac                                                                                                                                      >  1:211565/1‑150 (MQ=255)   attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAgg                                                                                                                                   <  2:197244/151‑1 (MQ=255)           ttaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCa                                                                                                                            <  1:114115/150‑1 (MQ=255)                           tGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtc                                                                                                           >  1:1052118/1‑151 (MQ=255)                             cacccTGTGCGCTGTCATTTTTTGCCTTTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCTGACGGAAAACCTTTTTTAATATTGTAGCCCTGCGCCCAGGTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTtca                                                                                                          <  1:206774/147‑1 (MQ=255)                                   gtgcgtgttCATTTTTGGCCATTTCTTTGTTTTCTGCCAAATTAAAATAGGTTCCGGACGACAACCGTTTTTTAATTTGGTAGCCCGCCGCCCAGTTTATGGGTGAAAAAAAAACCAGGCCCGCCAGTAAAGCAAAtttt                                                                                                              <  2:490819/132‑1 (MQ=255)                                       gcTTTCATTTTTTGCCTTTTCGTTGGTTTCTGCGCAAGTAATATGGGTTTCTGACGGCAACCCTTTTTTAATTTTGTAGTCCGGCCCCCAGGTCATTGGGGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCAtt                                                                                                  <  1:445153/148‑1 (MQ=255)                                         tttCATTTTTTGCCTTTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGAAACCCGTTTTTAATATTGTAGTCCTGCGCCCAGTTATTTGGTGCAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                <  1:657449/146‑1 (MQ=255)                                           ttattttttGCCTTTTCTTTGGCTTCGGCGCAAGTAATATAGGTTTCTGACGCCAAACCTTTTTAAATTTTGTAGTCCTGCGCCCAGTTCTTTGGGGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                <  2:967936/142‑1 (MQ=255)                                            cattttttGCCATTTCTTTGTTTTCTGCGAAAGTAATATAGGTTTCTGACGGAAACCCTTTTTAAATATTGTAGCCCTGCCCCCAGTTTATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTa                                                                                             <  2:525556/145‑1 (MQ=255)                                                tttcccctttttGTGGGTTTCTGCAAAGGAAATATGGGTTTCCGAGGACAACCGTTTTTTAATTTTTTAGTCTGGCGCCATGGTCTTTGTTGCGAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAt                                                                                       <  2:1042615/142‑1 (MQ=255)                                                  ttGCCTTTTTGTTGTCTTCTGCGCAAATAATATAGGTTTCGGAAGGAAACCCTTTTTTAATATTGTAGCCCTGCCCCCGGGTTTTTGGTGAAAAAACAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                      <  2:788879/149‑1 (MQ=255)                                                     cattttCGTTGTTTTTGGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATTTTGTAGCCCTGCGCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                 <  1:34475/148‑1 (MQ=255)                                                       aTTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                >  2:1028600/1‑150 (MQ=255)                                                          tCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGa                                                                             >  2:297146/1‑150 (MQ=255)                                                             ttGGCTTCTGCCCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTCCCCCCAGTTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                            <  1:421786/148‑1 (MQ=255)                                                                             tAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                         >  2:385363/1‑151 (MQ=255)                                                                              aaTAGAGGTTTCGGCCGGCAAACTTTTTTTAATATTGTGGCCCTGCCCCCAGGTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCa                                                        <  2:572446/151‑1 (MQ=255)                                                                                           gACGGCAAACCGTTTTTAATATTGTAGCCCGGCGCCCAGGTCATTGGGGAGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCCCGGCACCt                                            >  2:983585/1‑150 (MQ=255)                                                                                                 aacccttttttAATTTTGTAGTCCTGCCCCCAGGTCATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                      <  1:141603/144‑1 (MQ=255)                                                                                                      gTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGCCACCTTTAATCGCCCCt                                >  2:817885/1‑151 (MQ=255)                                                                                                           tAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                             >  1:802322/1‑149 (MQ=255)                                                                                                              tatTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                        >  1:764072/1‑151 (MQ=255)                                                                                                                ttGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAt                      >  1:889128/1‑151 (MQ=255)                                                                                                                  gTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTg                    >  2:1051811/1‑151 (MQ=255)                                                                                                                             cccAGGTAATTGGAGCGAAAAAAAACAGACCCGCCAGTGAAGCAACTTTTTTCATCATCATTCCTTATTTCACTTTACCCAGAATTGCACAACCAGCCCCGCCAATCACGGCACCTTTAATCGCCCCTCAGAGGCAATTGCCGGTCTGa           >  2:427036/1‑149 (MQ=255)                                                                                                                                    cATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg  >  1:906475/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                              GCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTG  >  NC_000913/1222456‑1222737

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