Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	23.7%	79.8 / 15.8	38	*559E (TAA→GAA)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (27/2);  new base C (0/9);  total (27/11)
Fisher's exact test for biased strand distribution p-value = 3.37e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGG  >  NC_000913/501432‑501701                                                                                                                                     |                                                                                                                                          tATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                    >  1:1204527/1‑140 (MQ=255) tATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                    >  1:847542/1‑140 (MQ=255)  aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                                                                   >  1:81970/1‑140 (MQ=255)    cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc                                                                                                                                 >  1:991076/1‑140 (MQ=255)     gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacca                                                                                                                                >  1:746256/1‑140 (MQ=255)       cgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCaccacc                                                                                                                              <  2:991076/140‑1 (MQ=255)         gTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                  <  2:54427/134‑1 (MQ=255)         gTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                  >  1:54427/1‑134 (MQ=255)                           gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                          >  2:898879/1‑140 (MQ=255)                              ttCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCcag                                                                                                       <  1:70157/140‑1 (MQ=255)                               tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                      >  1:1155875/1‑140 (MQ=255)                                cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                     >  1:911191/1‑140 (MQ=255)                                cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                     >  1:55933/1‑140 (MQ=255)                                 ccgccgcATCTGGCTTTTTGCTCTCGATGCCTGAGGCGACGCGGGCGCTTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagcag                                                                                                    <  2:837729/139‑1 (MQ=255)                                     cgcATCTGGCATTTTGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                >  1:837729/1‑140 (MQ=255)                                                  tgcctctCGATGCCTGATGCGACGCTGGCGCGTCTTACCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                   <  2:998670/137‑1 (MQ=37)                                                  tGGCTCTCGATGCCTGATGCGCCGCTGGCGCTTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                   <  2:746256/140‑1 (MQ=255)                                                  tGGCTCTCGATGCCTGATGCGACGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                   <  2:1057454/140‑1 (MQ=38)                                                               cTGATGCGACGCTGGCGCGTCTTATCATGCCTACACATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                      >  1:340698/1‑140 (MQ=255)                                                                 gATGCGACGCTGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCTGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                    <  2:664793/140‑1 (MQ=255)                                                                 gATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCACCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                    <  1:898879/140‑1 (MQ=255)                                                                  aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                   >  1:450318/1‑140 (MQ=255)                                                                   tGCGCCGCGGGCGCTTTTTTCCACGCCTCCTTTTTTTTCATTTTTTACACCCGGCACCCACCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                  <  2:383758/140‑1 (MQ=255)                                                                           tGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCAACCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                          <  2:984215/140‑1 (MQ=255)                                                                                  tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                   >  2:121286/1‑140 (MQ=255)                                                                                  tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                   >  1:16627/1‑140 (MQ=255)                                                                                       tCATGCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                              <  2:911191/140‑1 (MQ=255)                                                                                          tGCCTACATTTTTTTCATTTTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                           <  1:599268/140‑1 (MQ=255)                                                                                             cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                        >  1:83543/1‑140 (MQ=255)                                                                                               aCATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                      <  2:576011/140‑1 (MQ=255)                                                                                                 ttttttttCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                    <  1:402354/137‑1 (MQ=255)                                                                                                     ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                >  1:849798/1‑140 (MQ=255)                                                                                                           tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                          >  2:1075141/1‑140 (MQ=255)                                                                                                              tttACACCCGCCACCCACCGTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                       <  2:81970/140‑1 (MQ=255)                                                                                                               ttACACCCGGCACCCACCGTTCCCCCCTCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                      <  2:83543/140‑1 (MQ=255)                                                                                                                tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                     >  1:832757/1‑140 (MQ=255)                                                                                                                   aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                  >  2:764289/1‑140 (MQ=255)                                                                                                                     ccGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                >  1:1060609/1‑140 (MQ=255)                                                                                                                        gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa             >  2:707484/1‑140 (MQ=255)                                                                                                                        gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa             >  2:1306142/1‑140 (MQ=255)                                                                                                                        gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa             >  1:869201/1‑140 (MQ=255)                                                                                                                        gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa             >  1:1011755/1‑140 (MQ=255)                                                                                                                        gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTGCATCGTCATa             >  1:952198/1‑140 (MQ=255)                                                                                                                                cGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc     <  2:869201/140‑1 (MQ=255)                                                                                                                                cGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc     <  1:764289/140‑1 (MQ=255)                                                                                                                                  tttCCCCCATCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg   <  2:450318/140‑1 (MQ=255)                                                                                                                                   ttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg  >  2:683105/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                          TATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGG  >  NC_000913/501432‑501701

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