Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,215,016	G→C	20.0%	P296P (CCG→CCC)	mexD →	multidrug RND transporter MexD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,215,016	0	G	C	20.0%	49.8 / 10.1	25	P296P (CCG→CCC)	mexD	multidrug RND transporter MexD
Reads supporting (aligned to +/- strand):  ref base G (12/8);  new base C (2/3);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CATCGTCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGAAGGTGATC  >  NC_002947/3214875‑3215154                                                                                                                                              |                                                                                                                                           cATCGTCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCaa                                                                                                                                    >  2:358384/1‑150 (MQ=255)   tCGTCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCCGGCGCCaact                                                                                                                                  <  2:599718/150‑3 (MQ=255)     gTCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATAGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGcc                                                                                                                                >  2:399766/1‑150 (MQ=255)      tCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCa                                                                                                                               >  2:465632/1‑150 (MQ=255)      tCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCa                                                                                                                               <  1:358384/150‑1 (MQ=255)                 aaCCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACcgcg                                                                                                                    >  1:408295/1‑150 (MQ=255)                    cccGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACg                                                                                                                 >  1:268596/1‑150 (MQ=255)                             tCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAg                                                                                                        >  1:207315/1‑150 (MQ=255)                                    ttcGCCGGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTCTAACCAGTCGGCCCGACTGCCTCGCCACGCGGCGGGGGCCGGTGGGGTGCAGTTTGCGCCGGGCGCCAATGCCATCCAGCCCGCGACGCTGGTCAAGGAGCGCt                                                                                                 <  2:531724/148‑1 (MQ=255)                                    tCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCt                                                                                                 >  1:356351/1‑150 (MQ=255)                                                                           agaACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGAc                                                                                                                  <  2:231084/94‑1 (MQ=255)                                                                           agaACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGAc                                                                                                                  >  1:231084/1‑94 (MQ=255)                                                                                tACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGTCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCg                                                     >  1:586410/1‑150 (MQ=255)                                                                                   aaCCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTg                                                  >  1:70409/1‑150 (MQ=255)                                                                                                  cTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCTCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTcg                                   <  2:356351/150‑1 (MQ=255)                                                                                                      aCGGCCACCCGGAGGTGGCCGGTGCGGTGCAGTAGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCt                               >  2:434729/1‑150 (MQ=255)                                                                                                        ggCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTtc                             <  1:465632/150‑1 (MQ=255)                                                                                                        ggCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTtc                             <  2:207315/150‑1 (MQ=255)                                                                                                           cACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTtcgtc                          >  1:135126/1‑150 (MQ=255)                                                                                                             cccGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGa                        <  1:366162/150‑1 (MQ=255)                                                                                                                            tgcggtgcAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGaa         <  2:135126/150‑1 (MQ=255)                                                                                                                                     attGGCGCCCGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCgg                                                    <  2:236392/97‑1 (MQ=255)                                                                                                                                     attGGCGCCCGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCgg                                                    >  1:236392/2‑98 (MQ=255)                                                                                                                                     attGGCGCCCGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGAAGGTGATc  <  1:427450/147‑1 (MQ=255)                                                                                                                                     attGGCGCCCGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGAAGGTGATc  >  2:427450/2‑148 (MQ=255)                                                                                                                                              gggCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGaa         <  2:94418/132‑1 (MQ=255)                                                                                                                                              gggCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGaa         >  1:94418/1‑132 (MQ=255)                                                                                                                                              |                                                                                                                                           CATCGTCCTGCGCGCCAACCCCGATGGCTCCAGCGTCCGCCTGGGTGATGTGGCGCGCATGGCCATCGGTAGCGAGAACTACAACCTGTCGGCGCGCCTGAACGGCCACCCGGCGGTGGCCGGTGCGGTGCAGTTGGCGCCGGGCGCCAATGCCATCCAGACCGCGACGCTGGTCAAGGAGCGCTTGGCCGAGCTGTCGCAGTTCTTCCCCGAAGGTGTCGAATACTCGGTGCCGTACGACACCTCGCGCTTCGTCGACGTGGCCATCGAGAAGGTGATC  >  NC_002947/3214875‑3215154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑