Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	55.9%	28.8 / 41.7	34	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (11/4);  new base C (0/19);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-05

TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501426‑501710                                                                                                                                           |                                                                                                                                                   ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCAGATTTTACATCCGGCAACCACCGTTTACCCCGTcaccac                                                                                                                                        <  1:789931/151‑1 (MQ=255)            tcgTTTCGTAGGCCAGAAAAGGCGTTCACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTCCATCCGGAAACCCCCGTTTACCCCGTCCCCCCCTCACccgccg                                                                                                                             <  1:536279/151‑1 (MQ=37)            tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                              <  1:173858/150‑1 (MQ=255)                      cgCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATCCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTTACCCGCCGGTGGCGttt                                                                                                                    >  1:211254/2‑150 (MQ=255)                            aTAGGGCTTTCCCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCcagcag                                                                                                             <  2:634663/151‑1 (MQ=255)                                   tttCCCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATTTTTTACATCCGGCAACCCCCGTTTCCCCCTTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                      <  1:749421/150‑1 (MQ=255)                                   gTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                       >  1:386002/1‑150 (MQ=255)                                       ccgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                  <  2:613848/150‑1 (MQ=255)                                                        tGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCGTGCCTACATTTTTTTCATTTTTTCCATCCGGCAACCACCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:409765/151‑1 (MQ=255)                                                        tGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTCCATCCGGCACCCCCCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:448934/151‑1 (MQ=39)                                                             ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                             >  2:185273/1‑150 (MQ=255)                                                                  tGCCTGATGCGCCGCGGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                       <  1:509150/151‑1 (MQ=255)                                                                    acTGATGCGCCGCTGGCGCGTTTTATGATGCCTACATATTTTTCAGTTTTTACATCGGGCAACCACCTTTTCCCCCTTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                      <  1:404909/149‑1 (MQ=255)                                                                      tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                   <  2:308551/151‑1 (MQ=255)                                                                       gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                   >  2:135067/1‑150 (MQ=255)                                                                        aTGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTCTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                  <  1:554287/150‑1 (MQ=255)                                                                         tGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                   >  1:386636/1‑148 (MQ=255)                                                                                                  ccTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                       <  1:191428/151‑1 (MQ=255)                                                                                                          tttttCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                               <  1:66018/151‑1 (MQ=255)                                                                                                               cATATTTTACATCCGGCACCCTCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                          <  2:233773/151‑1 (MQ=255)                                                                                                                  atTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  2:276933/1‑151 (MQ=255)                                                                                                                   ttttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                      >  1:444861/1‑151 (MQ=255)                                                                                                                      tccacccGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                   <  1:430933/146‑1 (MQ=255)                                                                                                                      tACATCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                   <  2:863002/151‑1 (MQ=255)                                                                                                                               caccccccGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc          <  1:518029/145‑1 (MQ=255)                                                                                                                                  ccaccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat       >  1:204320/1‑151 (MQ=255)                                                                                                                                     ccGTTTACCCCGTCACCCCCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    <  1:408135/151‑1 (MQ=255)                                                                                                                                     ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATCATGGGCGCGGGCaataat    >  1:500085/1‑151 (MQ=255)                                                                                                                                     ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    >  1:811230/1‑151 (MQ=255)                                                                                                                                     ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    >  1:209594/1‑151 (MQ=255)                                                                                                                                      cttttCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc   <  1:678505/149‑1 (MQ=255)                                                                                                                                      cttttCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc   <  1:202875/149‑1 (MQ=255)                                                                                                                                      cGTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc   <  2:468021/151‑1 (MQ=255)                                                                                                                                       ttttCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct  <  2:576406/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                   TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501426‑501710

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