Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	869,107	0	A	G	46.5%	42.3 / 13.5	43	I530V (ATT→GTT)	gsiA	glutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base G (20/0);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 1.09e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGT  >  NC_000913/868960‑869256                                                                                                                                                    |                                                                                                                                                              aaTCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTa                                                                                                                                                            <  2:501379/151‑1 (MQ=255)     cAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGGGGGtt                                                                                                                                                             >  2:195867/1‑146 (MQ=255)      aaaaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGttt                                                                                                                                                            >  1:914383/1‑145 (MQ=255)         aGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCgt                                                                                                                                                    <  2:686285/151‑1 (MQ=255)         aGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCg                                                                                                                                                     >  1:334728/1‑150 (MQ=255)            gATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTACGCATtgtgtg                                                                                                                                                 >  2:101207/1‑151 (MQ=255)               cATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGATAGCGGATTTGATATTGTGGGGC‑g                                                                                                                                             >  1:490398/1‑151 (MQ=255)               cATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGGGCGGGTTTGTCATCGtgtg                                                                                                                                                 >  2:673238/1‑148 (MQ=255)                    ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGAGGGTTTGTCATCGTGTGgg‑ggggg                                                                                                                                         >  2:823388/1‑147 (MQ=255)                      gacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTGGAGCGGACTAGTTATTGTGTGGC‑ggggg                                                                                                                                         >  2:878990/1‑145 (MQ=255)                       acgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGATAGTTATCGTGgggg‑gg                                                                                                                                            >  1:106768/1‑144 (MQ=255)                         gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGACATCGTGTGGG‑GGGGATGt                                                                                                                                      >  1:408854/1‑148 (MQ=255)                             ccGTTTCGGCTGTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑gg                                                                                                                           <  2:621721/151‑1 (MQ=255)                             ccGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑gg                                                                                                                           <  2:79787/151‑1 (MQ=255)                                 ttCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑gggggtttatttccgg                                                                                                                              >  1:888292/1‑134 (MQ=255)                                 ttCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGgg‑gggggtattactctgg                                                                                                                              >  1:732104/1‑134 (MQ=255)                                      cgcTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAAttgt                                                                                                                  <  1:624214/151‑1 (MQ=255)                                           gATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGGGGGTTAGTCATCGTGTGGC‑GGGGGTGTA‑TT‑T‑C‑GGGGCAAatgtt                                                                                                                 >  2:592510/1‑147 (MQ=255)                                             tGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGC‑GGGGGTGTA‑TT‑T‑C‑GGGGGAAATTTTGaaaat                                                                                                           >  1:269260/1‑149 (MQ=255)                                                    aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGGCATCGTGGGGC‑GGGGTTGGT‑TC‑T‑C‑GGGCAAAtatttgtaa                                                                                                             >  1:94272/1‑134 (MQ=255)                                                       cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑gg                                                                                                 <  1:934451/150‑1 (MQ=255)                                                        gcgGGCAGATTATCAACTTGTTGCTCGAACTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑Ggcg                                                                                               <  1:750351/151‑1 (MQ=255)                                                          gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑Ggcgcg                                                                                             <  1:873862/151‑1 (MQ=255)                                                             cAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGCGGTAGAGCGGGTTAGTCATCGTGTGGG‑GGGGGTGTT‑TT‑T‑C‑GGGGGAAATGTTTTAATTTGtgcccaggccccgg                                                                                           >  2:729574/1‑138 (MQ=255)                                                                aTTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑t                                                                                      <  1:80261/151‑1 (MQ=255)                                                                   aTCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTc                                                                                    <  2:672906/150‑1 (MQ=255)                                                                     cAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGCGGGAGAGCGGGTTTGTCATCGTGTGGC‑GGGGGTGTA‑TT‑T‑C‑GGGGGGAATTTTGAAATTTGTGCCC‑Ggcccgccggcgt                                                                                      >  2:357282/1‑141 (MQ=255)                                                                      aaCTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGGGTGGG‑GGGGATGTG‑TC‑T‑C‑CGGGAAAATGTTTTAATTTGTTCTC‑CGCGCGCggggt                                                                                       >  2:768589/1‑142 (MQ=255)                                                                       aCTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGTGTA‑TC‑T‑C‑GGGGGAATTTTTTTAATTTGTGCAC‑CGCGCGCGGGC‑TTCGgaa                                                                                >  1:248244/1‑150 (MQ=255)                                                                          tgttgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAAc                                                                              <  2:199657/149‑1 (MQ=255)                                                                                 cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGTGAGCGGGTTAGTGATCTTGTGGG‑GGGGGGGTA‑TT‑T‑CGGGGCAAATGGTTTAAATCGTTCCCC‑GCCCGCCGGCT‑TTCGAAAACCCccac                                                                        >  1:273517/1‑148 (MQ=255)                                                                                  gATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCCTCGTGTGGCGGGGGTTGTT‑TC‑C‑C‑GGGGCAATGTTTGGATTTGTTCCCC‑GG‑GCGCGGGCTTTCGAAAAcccccgactc                                                                     >  2:286720/1‑144 (MQ=255)                                                                                  gATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGTGTA‑TT‑T‑C‑GGGGAAATTTTTTAAATTTGGCCCC‑GGCGCGCGGGG‑TTCGAAAAacccc                                                                          >  1:24757/1‑143 (MQ=255)                                                                                         aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGTGGA‑TCTT‑C‑GGGAAAATGTTTTGAATAGGT‑CAC‑AGCGCGCGGGC‑TTCGTAAAACCgccaccacc                                                                   >  2:905605/1‑138 (MQ=255)                                                                                          gCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCCTCGTGTGGC‑GGGGGTGTT‑TT‑T‑CGGGGGAAATGTTTTAAATTGGGCCCC‑GGGGGCGGGGC‑TTTGAAAAcccccgcctcccttt                                                                >  2:745609/1‑136 (MQ=255)                                                                                                            gTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCa                                         <  1:236273/151‑1 (MQ=255)                                                                                                                 tGTTTATCTCCCACGATATGGCGGCGGGAGAGCGGATTAAGCATCGTGTGGC‑GGGGGAGTATTT‑T‑C‑GGGGAAATTTTGGAAATTGGT‑CAA‑GGGGCGCGGGC‑TTTGAAAAAACCCCGCAACCCTATATGC‑GGAAATTaatagggggggtt                                      >  2:278646/1‑137 (MQ=255)                                                                                                                  gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGtt                                      <  1:770079/148‑1 (MQ=255)                                                                                                                   tttATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGTGTA‑TT‑TCC‑GGGGAAATTTTTAAAATT‑GTTCAC‑CGCGCGCGGGCTTTCGAAAACCC‑CCGCCTCCCTTTACGC‑GGAAAATAATtgcggccggt                                      >  1:894184/1‑143 (MQ=255)                                                                                                                     tATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGG‑GGGGGTGTA‑TT‑TCC‑GGGAAAATTTTTGAAATTGGTCCCCGGGCCCGGGGT‑‑TTGGAAAACCC‑CAGCAACCTTTTAAGCGGAAAATTATTGGGGGCGGTCCCg                                   >  1:628671/1‑149 (MQ=255)                                                                                                                                    ggCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTccc                   >  1:633448/1‑149 (MQ=255)                                                                                                                                             agagCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGAcc         >  1:596834/1‑150 (MQ=255)                                                                                                                                                   gATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCgt  <  2:281015/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                              AATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGATGTA‑TC‑T‑C‑GGGCAAATTGTTGAAATTGGTCCAC‑GGCGCGCGGTC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGT  >  NC_000913/868960‑869256

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: ‑