Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	27.9%	62.4 / 24.7	43	*559E (TAA→GAA)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (30/1);  new base C (0/12);  total (30/13)
Fisher's exact test for biased strand distribution p-value = 8.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGG  >  NC_000913/501426‑501700                                                                                                                                           |                                                                                                                                         ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAc                                                                                                                                         >  1:585577/1‑140 (MQ=255)     gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCg                                                                                                                                     >  1:709111/1‑140 (MQ=255)      ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGt                                                                                                                                    <  2:238142/140‑1 (MQ=255)        aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                                                                  >  1:384056/1‑140 (MQ=255)                ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                          >  1:1228935/1‑140 (MQ=255)                         cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                 >  1:1094227/1‑140 (MQ=255)                                 gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                         >  2:1071676/1‑140 (MQ=255)                                    tccccgccgcTTCTGGCTTTTGCCTCTCGATCCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTCCACCCGCCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                      <  2:1278440/136‑1 (MQ=255)                                    tccccgccgcCTCTGGCTTTGGGCTCTCGATGCCTGAGGCGACGCTGGCCCGTTTTTTTACCCTTACATTTTTCTAATTTTTTACCCCCGACACCCTCCTTTTCCCCCTTCACCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                                                      <  2:709111/136‑1 (MQ=255)                                     ccccgccgcATCTGTCTTTTGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTAAATTTTTTTCATTTTTTACACCGGGCAACCCCCTTTTCCCCCTTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                                                                     <  2:1228935/137‑1 (MQ=255)                                      cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                    >  1:858143/1‑140 (MQ=255)                                          ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                >  2:1352339/1‑140 (MQ=255)                                             cATCTGTCTTTTGCCTCCCGATGCGTTACGCGACGCTGGCGCGTCTTATCTTGCCTACATTTTTTTAATATTTTCCACCCAGCAACCACCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                             <  2:672693/140‑1 (MQ=255)                                               tCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                           <  2:946075/140‑1 (MQ=38)                                               tCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                           >  1:299325/1‑140 (MQ=255)                                                             ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgctgtctc                                                                             >  1:373743/1‑133 (MQ=255)                                                                  tGCCTGTCGCGCCGCTGGCGCCTTTTATCATCCCTACTTTTTTTTCATATTTTACTCCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                        <  2:384056/140‑1 (MQ=255)                                                                   gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                       <  1:1213230/140‑1 (MQ=38)                                                                    cttagtGCGACGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTCACACCCGCCAACCCCATTTTCCCCCAACACCCCCACCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                      <  2:589101/135‑1 (MQ=255)                                                                    ccTGATGCGCCGCTGGCGCTTCTTATCATCCCTACTTTTTTTTCATATTTTCCCCCCGCCAACCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                      <  2:689881/140‑1 (MQ=255)                                                                    ccTGATGCGCCGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                      <  2:124648/140‑1 (MQ=255)                                                                       gATGCGGCGCTGGCGCGTCTTATCATGCCTATATTTTTTTAATATTTCACATCCGGCAACCCCCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                   <  2:1269747/140‑1 (MQ=255)                                                                        aTGCGCCGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCACCTTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                  <  2:1123462/140‑1 (MQ=255)                                                                        aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                  >  1:627295/1‑140 (MQ=255)                                                                            gACGCTGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                              <  1:1293188/140‑1 (MQ=255)                                                                             aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                             >  1:84965/1‑140 (MQ=255)                                                                                    cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                      >  1:662011/1‑140 (MQ=255)                                                                                             tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                             >  1:1261602/1‑140 (MQ=255)                                                                                             tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                             >  2:1155740/1‑140 (MQ=255)                                                                                              cATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                            <  2:858143/140‑1 (MQ=255)                                                                                                 gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                         >  1:1090607/1‑140 (MQ=255)                                                                                                 cccTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                         <  2:299325/139‑1 (MQ=255)                                                                                                    tACATTTTTTTCATATTTTACACCCGAAACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                      <  2:1161299/140‑1 (MQ=255)                                                                                                    tACATATTTTTCATATTTTACATCCGGAAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                      >  1:1161299/1‑140 (MQ=255)                                                                                                     aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                     >  1:336832/1‑140 (MQ=255)                                                                                                      cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                    >  1:799132/1‑140 (MQ=255)                                                                                                      cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                    >  2:517333/1‑140 (MQ=255)                                                                                                      cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                    >  1:1268855/1‑140 (MQ=255)                                                                                                       atatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                   >  2:479520/1‑140 (MQ=255)                                                                                                          tttttCATTTTTTCCACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                <  2:84965/140‑1 (MQ=255)                                                                                                           ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                               >  1:832145/1‑140 (MQ=255)                                                                                                             ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                             >  1:160138/1‑140 (MQ=255)                                                                                                              tCATTTTTTACATCCGGCACCCACCGTTTCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                            <  1:1352339/140‑1 (MQ=255)                                                                                                                atatTTTACATCCGGCAACCCCCGTTTCCCCCTTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                          <  2:627295/140‑1 (MQ=255)                                                                                                                    tttACACCCGGCAACCCCCTTTCCCCCCGCCACCCCCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                      <  2:720369/140‑1 (MQ=255)                                                                                                                      tACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  1:713638/140‑1 (MQ=255)                                                                                                                      tACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  2:1455550/140‑1 (MQ=255)                                                                                                                      tACACCCGGCAACCCCCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  2:585577/140‑1 (MQ=255)                                                                                                                          ccgggCACCCACCGTTACCCCCGTCCCCCCCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                <  2:925176/137‑1 (MQ=255)                                                                                                                            cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa              >  1:1356687/1‑140 (MQ=255)                                                                                                                             ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGAGATATACGCCACTTCATCGTCAt             >  1:925176/1‑140 (MQ=255)                                                                                                                                     ctttttCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg     <  2:1145987/137‑1 (MQ=255)                                                                                                                                     ccGTTTCCCCCGCCCCCCCCCCAACCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg     <  2:1135772/140‑1 (MQ=255)                                                                                                                                     ccGTTTACCCCGTCACCACCTCACCCGCCGGGGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCCTCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGTgcg     >  2:165919/1‑140 (MQ=255)                                                                                                                                     ccGTTCCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg     <  2:1268855/140‑1 (MQ=255)                                                                                                                                      cGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc    <  1:517333/140‑1 (MQ=255)                                                                                                                                       ttttCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg   <  1:1512111/139‑1 (MQ=255)                                                                                                                                       gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGTAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTA‑CTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg  >  2:1094833/1‑140 (MQ=255)                                                                                                                                        tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg  >  2:560345/1‑140 (MQ=255)                                                                                                                                        tttACCCCGTCACCAACTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg  >  2:133798/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                         TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGG  >  NC_000913/501426‑501700

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