Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,483,254	0	G	A	43.2%	50.2 / 19.9	45	S29L (TCA→TTA)	emrK	multidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base A (0/19);  total (12/33)
Fisher's exact test for biased strand distribution p-value = 3.39e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCAT  >  NC_000913/2483108‑2483403                                                                                                                                                   |                                                                                                                                                        gTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑aaa                                                                                                                                                     >  1:784830/1‑150 (MQ=255)    gAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACg                                                                                                                                                  >  2:180227/1‑150 (MQ=255)     aaCGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCa                                                                                                                                                <  2:338123/151‑1 (MQ=255)             tCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCaataaata                                                                                                                                         <  1:668071/150‑1 (MQ=255)              cGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATT‑CCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAAtaaataaa                                                                                                                                       <  2:988857/150‑1 (MQ=255)                    tttATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTAcc                                                                                                                                >  1:1176299/1‑151 (MQ=255)                      tATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCg                                                                                                                               <  2:937323/150‑1 (MQ=255)                       aTGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGc                                                                                                                              <  1:1259875/150‑1 (MQ=255)                       aTGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGc                                                                                                                              <  2:1060759/150‑1 (MQ=255)                        tGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGcc                                                                                                                             >  2:909326/1‑150 (MQ=255)                         gATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCa                                                                                                                            <  2:614448/150‑1 (MQ=255)                             aacaacAGTGACCCTACCTGAGTTTTGTGAAAAAATTGGATCTGCATCCCCCGTGACATAGGCGTCTTCTGTACAAATCATTTCTTCTAATCCCATTGACAAATAGGAATAGCCACCTG‑AAAACGAAAAAAATAAAA‑CTACCGCCAATa                                                                                                                         <  1:43546/149‑1 (MQ=255)                                 aCAGTGCCTCAACCGGAGTTTTGCGCAAAAATTGGATCTACATTCCCCGTGACATAGGCGTCATCTGTAATAATCTTGTTTTAAAATTCAATTCACCAAAAGCAATGGCCACCTA‑AAACCGCAATAAAAAAAA‑CTACCGCCAATAAAGaaa                                                                                                                   <  1:366009/151‑1 (MQ=255)                                  cAGTGACACTCCCTGAGATTTGTGAAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATTTCTTCTAATTCCATTGACAAAAAGAAATAGACCCCTA‑AAAACGCAAAAAAAAAAA‑CTACCGCCAATAAAGaaa                                                                                                                   <  1:873076/150‑1 (MQ=255)                                   aGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGaaa                                                                                                                   >  1:209034/1‑149 (MQ=255)                                                gagttttGTGCAAAAATTGTATCTGCATTCCCCGTAACATAGGCGTCATCTGTACTAATCATGTTTTCAAATCCCATTGCCAAAAAGGAATAGGCACCTA‑AAAACACAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTttctt                                                                                                        <  1:944302/142‑1 (MQ=255)                                                 agaTTTGTGCAAAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATCCCATTGACAAATAGGCATAGCCACCTA‑AAAACGAAAAAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                   <  2:856320/151‑1 (MQ=255)                                                     ttGTCCAGAAATTGGATCTGCATTCCCCGTGACATGGGCGTCATCTGTACTAATCATGTCTTCAAATTCCATGAACAAATAGGAAAAGGC‑CCTGAAAAACGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                <  2:1239895/150‑1 (MQ=255)                                                      tgtgCAAAAATTGTTTCTCCATTCCCCGTGAAAGGGTCGTTTTCTAACATAATGTTGTTTTATAATTCAATTATCCAATAGGC‑TGAGCCTCTT‑AAAACCCAAAAAAAAAAATTTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                                                                                 <  2:48982/148‑1 (MQ=255)                                                                     tCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAt                                                                               >  1:927181/1‑151 (MQ=255)                                                                       tGCATTCCCCGTAAAATAGGCGTCATCTTTACTAATCATGTTTTCAAATTCCATTGACAAAAAGGCATAGGCACTTA‑AAAACGAAAAAAAAAAAA‑TTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGaa                                                                                <  1:395563/148‑1 (MQ=255)                                                                        gCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                                                                             >  2:741534/1‑150 (MQ=255)                                                                           tcccccGTAACATGGGCGTCATTTGTACTAATAATTTCTTCAAATCCCATTGACAAATAGGCATAGACACCTA‑AAAACGAAATAAATAAAT‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                              <  1:306687/144‑1 (MQ=255)                                                                            tCCCCGTGACATGGGCGTCATCTGTACTAATCATTTCTTCAAATTCCATTGACCAATAGGCAGGGGCCCTAA‑AAAACGCAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                         <  2:1234713/150‑1 (MQ=255)                                                                             ccccGTGACATAGGCGTCATCTTTACTAATCATGTCTTCTAATCCCATTGCCCAATGGGCATAGGCACCTG‑AAAACGAAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCt                                                                         <  2:972281/149‑1 (MQ=255)                                                                             ccccGTGACAAAGGCGTAATCTTTTCAAAAATTTTTTTCAAATCCCATCGACAAATAGAAAGAGACACTTA‑AAAACGAAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAAt                                                                           <  1:404884/147‑1 (MQ=255)                                                                             ccccGTAACATAGGCGTCTTCTGTATTAATCATTTTTTCTAATCCCATTGCCAAAAAGGAATAGGCACCTA‑AAACCAAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                        <  1:555549/150‑1 (MQ=255)                                                                              cccGTGACATGGGCGTCATCTTTACAAATCTTTTTTTCAAATTCCATTGACCAATGGGCAAAGCCACCAA‑AAAACGCAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGtt                                                                      <  2:317411/151‑1 (MQ=255)                                                                                  tGAAATAGCCTTCATCTGTCAAAATCTTGTTTTCAAATCCCATAAACAAATAGCCATAGGCCCTTG‑AAAACGAAAAAAAAAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa                                                                   <  1:463403/150‑1 (MQ=255)                                                                                      aTAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat                                                               <  1:1125451/150‑1 (MQ=255)                                                                                         ggCGTCATCTGTACTAATCATGTCTTATAATTCCATTGACCAATAGGCATAGGCACTTA‑AAAACGAAAAAAAAAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtattat                                                            <  1:732837/150‑1 (MQ=255)                                                                                          gCGTCATCTGTACTAATCATTTCTTCTAATTCCATTGACCAAAAGGAATAGGCACCTA‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtct                                                          <  1:95856/151‑1 (MQ=255)                                                                                                tCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAtt                                                    >  2:259990/1‑151 (MQ=255)                                                                                                   gTACTAATCATTTCTTCTAATTCCATGGACCAATAGGCATAGGCACCTA‑AAACCGCAAAAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                                  <  2:1067825/150‑1 (MQ=255)                                                                                                    tACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtctc                                                >  1:1237372/1‑151 (MQ=255)                                                                                                     acaaaTCATGTCTTCTAATCCCATTGACCAATAGGCATAGGCACCTG‑AAAACGAAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAt                                                     <  2:297977/142‑1 (MQ=255)                                                                                                        aaTCATGTCTTCTAATCCCATGGACAAATAGGAATAGGCCCCTA‑AAACCGAAATAAAAAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCAt                                              <  1:1101816/149‑1 (MQ=255)                                                                                                             tGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTA‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGAt                                          <  2:413240/148‑1 (MQ=255)                                                                                                                     aaTTCCATTGACAAATAGGAATAGGCACCTA‑AAAACGAAATAAATAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGaaa                                      <  1:37968/144‑1 (MQ=255)                                                                                                                                ccAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatgaa                     >  1:1259213/1‑150 (MQ=255)                                                                                                                                caaaagggaagaggcaCCAA‑AAAACGCAATAAATAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatga                      <  2:913375/138‑1 (MQ=255)                                                                                                                                  aataggaataggCCCCTA‑AAAACAAAAAAAAAAAAT‑TTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                    <  1:553170/149‑1 (MQ=255)                                                                                                                                  aaTAGGCATAGGCACCTA‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaaa                   <  2:181687/150‑1 (MQ=255)                                                                                                                                     aggcataggcaCCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGc               >  2:249940/1‑151 (MQ=255)                                                                                                                                                  tG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat  >  2:300029/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                        GTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTG‑AAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCAT  >  NC_000913/2483108‑2483403

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