Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	43.6%	64.5 / 10.2	39	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (19/3);  new base C (0/17);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 2.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066                                                                                                                                                      |                                                                                                                                                      aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACAACCGCTGCGGCccacc                                                                                                                                                     >  1:207759/1‑150 (MQ=255)   cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCaccacc                                                                                                                                                  >  2:777485/1‑151 (MQ=255)    aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                                 >  1:744716/1‑151 (MQ=255)     gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGc                                                                                                                                                <  1:65323/151‑1 (MQ=255)             tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                        >  2:62849/1‑151 (MQ=255)                    aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                                 >  2:399297/1‑151 (MQ=255)                          tcggcgTTTCTTCCAGATTTTCAGCCGGGGGTGGGGTCCCCGCTAACCCCCCGCGGGGGGTG‑GGCTTTCCCAGCCCGCCCACAGGGCGCAGCAGGGTTCTTTTG‑CCCCCCCCGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGa                                                                                                                             <  1:485757/146‑1 (MQ=255)                           tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                           >  1:496982/1‑150 (MQ=255)                                tgtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                      >  1:449379/1‑150 (MQ=255)                                tgtCTCCCTAATTTCCACCCGGGGGTGTGGTCCCCGCTACCCCACCCCCGGCGTTG‑GGCTTTCCCAGCCCGGCCAGAAGGCGCAGCGGGTTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                           <  1:463936/145‑1 (MQ=255)                                   cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                  >  2:201286/1‑151 (MQ=255)                                      cctaattttCACCCGGGGGTGTGGT‑CCCCCCAACCCCCCCCCTGGGGTTGGGCTTTCCCAGCCCCCCCAGCAGGCGCAGCAGGGTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                          <  1:830134/134‑1 (MQ=255)                                       ctaattttCGGCCGGCGGTGTGGCGCCCCCTAACACCCGGCCGGGGTTT‑GGCTTTTCAAGCCCGGCGAGGGGGCAAAGCGGGGT‑CTTTTGCCCCCCCCCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                                <  2:571470/144‑1 (MQ=255)                                             ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                        >  2:206380/1‑151 (MQ=255)                                              tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                        >  2:201661/1‑150 (MQ=255)                                              tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                        >  2:310960/1‑150 (MQ=255)                                                    agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                  >  2:801188/1‑150 (MQ=255)                                                     ggggTGGGGTCCCCGCAAACACACCGCCGGCGTTG‑GGGGTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTCCTTTTC‑CCCCCCCCGCTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                        <  1:34236/141‑1 (MQ=255)                                                     gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                                >  1:659463/1‑151 (MQ=255)                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              >  1:107216/1‑151 (MQ=255)                                                         tgtgGTGCCCGCTAACACATCGCCTGGGGTT‑GGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTACTTTTC‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                             <  2:123035/150‑1 (MQ=255)                                                             gtgCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                        >  1:685477/1‑151 (MQ=255)                                                                 ccGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                     >  2:656850/1‑150 (MQ=255)                                                                 ccGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                     >  1:112426/1‑150 (MQ=255)                                                                      aaCACACCCCCTGCGGTT‑GGCGTTCCCACCCCTGCCGGCAGGCGCAGCAGGGTCCTTTTC‑CCCCCCCCGCTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                                 <  2:44754/149‑1 (MQ=255)                                                                                         ggCTTTCCCAGCCCTGCCAGCGGGCGCAGAGGGTTACTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                              <  2:135413/150‑1 (MQ=255)                                                                                         ggCTTTCCCAGCCCTGCCAGCAGGCGCAGGAGGTTACTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                               <  2:258417/149‑1 (MQ=255)                                                                                            tttcccAGCCCTGCCAGCGGGCGCAGCGGGTTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                             <  2:731017/144‑1 (MQ=255)                                                                                                cgagcccTGCCGGAAGGCGCGGAAGGGTCCTTTTC‑CCGCCCCGCCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                         <  2:832390/143‑1 (MQ=255)                                                                                                    aCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                  >  1:535414/1‑151 (MQ=255)                                                                                                         ccAGCAGGCGCAGCGGGGTATTTTGG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                             <  2:774204/151‑1 (MQ=255)                                                                                                               gccgcAGAGGGTTATTTTTC‑CCCCCCCCGGTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                             <  1:2851/143‑1 (MQ=255)                                                                                                                   cagcggGTTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                      <  2:472732/148‑1 (MQ=255)                                                                                                                    agcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                    >  1:721172/1‑149 (MQ=255)                                                                                                                      cggggTTCTTTCG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                                 <  2:286033/148‑1 (MQ=255)                                                                                                                         gttttttttC‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                             <  2:166658/146‑1 (MQ=255)                                                                                                                         ggttcttgcg‑CCACCCCCGCTGCCACCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                               <  2:412476/140‑1 (MQ=255)                                                                                                                            aCTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                          >  2:197166/1‑151 (MQ=255)                                                                                                                               tttC‑CCCCCCCCGCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                        <  1:857060/150‑1 (MQ=255)                                                                                                                                   cccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                     <  1:79046/143‑1 (MQ=255)                                                                                                                                    ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   >  2:915512/1‑151 (MQ=255)                                                                                                                                                 gcccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa        <  2:769669/147‑1 (MQ=255)                                                                                                                                                   ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa      <  2:68774/149‑1 (MQ=255)                                                                                                                                                      accaccGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  >  2:961941/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066

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