Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,483,268	0	T	A	36.4%	53.1 / 11.1	34	L24F (TTA→TTT)	emrK	multidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (17/4);  new base A (0/12);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 5.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGC  >  NC_000913/2483061‑2483469                                                                                                                                                                                                                 |                                                                                                                                                                                                          tGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAttt                                                                                                                                                                              >  2:748699/1‑237 (MQ=255) tGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGGAttt                                                                                                                                                                              >  2:652161/1‑237 (MQ=255)  gCGATAGTGGCATCAGTTTTATCCAGTAAAACTAAAATGTCACCTTGTCGAACGTAGTGCGTATCTTTATGATTAACGACAGTGACACTCCCTGAGACTTGTGGAGAAATTGGATCTGCATCCCCCGGAAC‑ATAGGCTTCATCTGTTCTAATCATTTTTTCTAATGCAATTGACAAATAGGCAAAGGCCCCGGAAAACGAAAAAAAAAAAACTGCCGCCAATAAAGAAAAGTATTttc                                                                                                                                                                            <  1:652161/238‑1 (MQ=255)      tggTGGAATCATTTTTATCCAGTAAAACTAAAATGTCCCTTTGTCGAACGTATTTCGTATTTTTATGATTAACGACAGTGACACTCCCTGAGATTTGTGAAAAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTTTTCTAATTCCATTGACAAAGAGGCATAGGCCCCTAAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                                                                                                                                                       <  1:748699/237‑1 (MQ=255)                      aTCCAGTAAAACTAAAATGTCCCTTTGCCGAACGTAGTTCGTATCTTTATGATAAACGACAGTGCCCCTACCGGAGATTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATGGGCGTCATCTGTACTAATCATGTCTTCTAATCCCATTGACCAATAGGCATAGGCACCTAAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                                                                                                                                                       <  1:76284/239‑1 (MQ=255)                                                                     tATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc                                                                                                                                       >  2:12882/1‑208 (MQ=255)                                                                     tATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAAAAAAACGACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc                                                                                                                                       <  1:12882/208‑1 (MQ=255)                                                                          ttAACGACGGTGACACTACCGGAGCTTTGTGCAGAAATTGTTGTTGCATCCCCGGTGAA‑ATAGGGGTCATCTGTACTAATCATTTCTCCTATTTCCTTTGACCAATGGTCATAGCCACGGAAAAACGAAATAAAAAAAACCACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                  <  1:527557/176‑1 (MQ=255)                                                                                    tGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTAAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGAATTTTCTTCTGTTAGAATGTTTTTTa                                                                                                                                                        >  2:527557/1‑176 (MQ=255)                                                                                      acacTCCCGGAGACTTGTGCAAAAATTGGATCTGCATCCCCCGTGAC‑ATAGGCTTCATCTGTACTAATCATTTCTTCTAATTCCATTGACAAATAGGAATAGGCCCCTAAAAACCAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                                                                                          <  2:629048/188‑1 (MQ=255)                                                                                           aCCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAaga                                                                                   >  1:851569/1‑238 (MQ=255)                                                                                            ccTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACAAATAGGCATAGGCACCTAAAAACGAAAAAAATAAAACTACCGCCa                                                                                                                                                                                              <  1:377278/130‑1 (MQ=255)                                                                                                agaCTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                                                                                                >  1:629048/1‑188 (MQ=255)                                                                                                      gtgCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGaaaa                                                                                                                                                                                    >  2:377278/1‑130 (MQ=255)                                                                                                      ctgCAGAAATTGGATCTGCATACCGCGCGAC‑AAATCCCTCATCTGTACTAATCATGTCTTCTAATTCCATTAACCAATAGACATGGTCAACTGACAACGCAATAAAGAAGACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCTACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATtctttctt                                                                       >  2:656649/2‑239 (MQ=255)                                                                                                       tgCAAAATTTGCTTCTGCATCCCCCG‑GACAATAGGCGTCATCTTTACTAATCTTTTCTTTAAATCCCATTAACAAATAGGCATAGGCCCCTGAAAACGAAAAAAAAAAAACCACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttcttc                                                                      <  2:710329/239‑1 (MQ=255)                                                                                                        gaaGAAATTGGATTTGCTTTCCCCGTCAC‑ATAGGCTTCATTTGTCCTAATCTTGTCTTTAAATTCCATTGACCAATAGGCATGGGCACCTAAAACCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAACCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttcttc                                                                      <  1:810110/236‑1 (MQ=255)                                                                                                               ttGCATCTGCTTACCCCGTGAC‑ATGGGCTCCATTTTTTGTATAATTGTCTTATATTTCCATTGACAAATAGGTAAAGGTACAAAAAATCGAAAAAAAAAAAACAACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAAACTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCaaa                                                              <  1:242381/239‑1 (MQ=255)                                                                                                                                                    tACTAATCTTTTCTTCTAATCCCATTCACGAATAGGAAAAGGCACCAAAAAACGAAAAAAAAAAACCACCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                                                                                                                                                       <  2:216909/114‑1 (MQ=255)                                                                                                                                                         aTCATGTCTTCTAATCCCATTACCAAATAGGCATAGCCACCTGAAAACGAAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc                                                                                                                                       <  2:552306/125‑1 (MQ=255)                                                                                                                                                              gTCTTCTAATTCCATTGTCCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTATTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCCTGTGTTGCACATTATCCCAGGATCGTTGGAGCCCTaa                >  2:776677/1‑239 (MQ=255)                                                                                                                                                              gTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAAt                                                                                                                                             >  1:216909/1‑114 (MQ=255)                                                                                                                                                                   ctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATg                                                                                                                                                               >  1:726725/1‑91 (MQ=255)                                                                                                                                                                   ctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtc                                                                                                                             >  1:552306/1‑125 (MQ=255)                                                                                                                                                                     aaTCCCATTAACAAAAAGGCATAGACAAAGAAAAACGTAAAAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAaca                                                            <  2:577699/188‑1 (MQ=255)                                                                                                                                                                            ttGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGc  >  2:367825/1‑239 (MQ=255)                                                                                                                                                                               aCCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCAtgt                                                  >  1:577699/1‑188 (MQ=255)                                                                                                                                                                                    taggcataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCAtgtg                                                 >  1:812317/1‑184 (MQ=255)                                                                                                                                                                                       gaagaggcaCCTAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCa                                                                                                                 <  2:510752/113‑1 (MQ=255)                                                                                                                                                                                         ataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGc                                                      <  1:735224/174‑1 (MQ=255)                                                                                                                                                                                         ataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGc                                                      >  2:735224/1‑174 (MQ=255)                                                                                                                                                                                           aggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAag                                                                                    >  2:11157/1‑142 (MQ=255)                                                                                                                                                                                           aggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAag                                                                                    <  1:11157/142‑1 (MQ=255)                                                                                                                                                                                                 cTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAt                                                                                                       >  1:510752/1‑117 (MQ=255)                                                                                                                                                                                                             taaataaaACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                                                                                                >  1:631110/1‑80 (MQ=255)                                                                                                                                                                                                             taaaaaaaACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                                                                                                <  2:631110/80‑1 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                          TGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGAC‑ATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGC  >  NC_000913/2483061‑2483469

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