Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	52.4%	20.9 / 31.4	21	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/11);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04

GATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAT  >  NC_000913/501452‑501708                                                                                                                            |                                                                                                                                      gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                            >  1:89639/1‑151 (MQ=255)                                     cGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                        >  2:313278/1‑151 (MQ=255)                                       aTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                      <  2:391596/151‑1 (MQ=39)                                          ccTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                    <  2:176589/150‑1 (MQ=255)                                          ccTGATGCGACGCGGGCGCGTCTTATCATGCCTACATATTTTTAATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                    <  2:434355/150‑1 (MQ=255)                                                 cgccGCTGGCGCGTCTTTTCATCCCTACATATTTTTCTTATTTTCCATCCGGAAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                            <  1:330376/148‑1 (MQ=255)                                                              tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                >  1:454144/1‑150 (MQ=255)                                                                 tATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                            >  1:178215/1‑151 (MQ=255)                                                                      tGCCTACATATTTTTCATATTTTCCATCCGGCACCCCCCGTTTCCCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                       <  2:462558/151‑1 (MQ=255)                                                                      tGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                       <  2:437613/151‑1 (MQ=255)                                                                           aCATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                  <  2:354976/151‑1 (MQ=255)                                                                           aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                  <  1:36031/151‑1 (MQ=255)                                                                            cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                 >  2:350963/1‑151 (MQ=255)                                                                             atttttttCATATTTTCCCTCCGGCAACCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                <  1:391609/148‑1 (MQ=255)                                                                                tttttCATATTTTCCACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCAGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                             <  2:251561/151‑1 (MQ=255)                                                                                  tttCATATTTTACATCCGGCACCCACCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                           <  1:110315/151‑1 (MQ=255)                                                                                        atTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                     >  1:184919/1‑151 (MQ=255)                                                                                          tttACATCCGGCAACCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATaa                   <  2:284143/151‑1 (MQ=255)                                                                                                   ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg          >  2:361402/1‑151 (MQ=255)                                                                                                           ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat  >  2:321888/1‑151 (MQ=255)                                                                                                           ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat  >  1:432454/1‑151 (MQ=255)                                                                                                                            |                                                                                                                                      GATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAT  >  NC_000913/501452‑501708

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