Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,213,301	G→A	9.3%	R123Q (CGG→CAG)	mexC →	multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,213,301	0	G	A	9.3%	116.3 / 2.8	43	R123Q (CGG→CAG)	mexC	multidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (2/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GTGTTGCACAAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCGCACGCCTGAACCTG  >  NC_002947/3213153‑3213440                                                                                                                                                     |                                                                                                                                            gtgtTGCACAAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCCGATCGCCCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGcgcg                                                                                                                                             <  6:143052/149‑1 (MQ=255)  tggtGCACAAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCgg                                                                                                                                            <  4:201941/146‑1 (MQ=255)      gCACAAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTg                                                                                                                                         >  8:67158/1‑148 (MQ=255)       cacaAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGaa                                                                                                                                       <  5:117698/149‑1 (MQ=255)              gcTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTAt                                                                                                                                <  4:302232/149‑1 (MQ=255)              gcTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTAt                                                                                                                                <  7:67158/149‑1 (MQ=255)                    aGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCa                                                                                                                          <  6:5158/149‑1 (MQ=255)                       aaGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTg                                                                                                                       <  3:179635/149‑1 (MQ=255)                         ggggCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGt                                                                                                                     >  1:220630/1‑149 (MQ=255)                                     aaGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGc                                                                                                         <  2:229457/149‑1 (MQ=255)                                       ggCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCAGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCgg                                                                                                       >  1:139815/1‑149 (MQ=255)                                                    cTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGa                                                                                          >  2:266560/1‑149 (MQ=255)                                                       ttCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGAc                                                                                         <  1:266560/147‑1 (MQ=255)                                                             aTCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAg                                                                                 >  5:147433/1‑149 (MQ=255)                                                                 aCCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCAGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAgcgc                                                                              >  6:223092/1‑148 (MQ=255)                                                                    cGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCAcc                                                                          >  2:255438/1‑149 (MQ=255)                                                                           gTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGc                                                                   <  6:147433/149‑1 (MQ=255)                                                                               aaGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCAccgccgcaccgg                                                                >  4:163450/1‑142 (MQ=255)                                                                                 ggCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCAGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCa                                                             <  5:223092/149‑1 (MQ=255)                                                                                  gCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAg                                                            >  2:86175/1‑149 (MQ=255)                                                                                   ccGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGt                                                           >  6:38002/1‑149 (MQ=255)                                                                                           gCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTgc                                                   >  6:202421/1‑149 (MQ=255)                                                                                            cccGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTgcg                                                  <  8:45729/149‑1 (MQ=255)                                                                                             ccGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTgcgg                                                 >  7:174866/1‑149 (MQ=255)                                                                                                  ccGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCGGTGCCGGTGCGGCGGTg                                            >  3:259951/1‑149 (MQ=255)                                                                                                       gCTGACCTGGCCCGAGCCCAGCCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGctc                                       >  4:57405/1‑149 (MQ=255)                                                                                                        cTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGctct                                      >  1:279074/1‑149 (MQ=255)                                                                                                          gACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGc                                    >  3:189409/1‑149 (MQ=255)                                                                                                           aCCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGcc                                   <  1:86175/149‑1 (MQ=255)                                                                                                           aCCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGc                                    <  1:255438/148‑1 (MQ=255)                                                                                                             cTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGcc                                   >  3:34891/1‑147 (MQ=255)                                                                                                             cTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGcc                                   <  4:34891/147‑1 (MQ=255)                                                                                                             cTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGcc                                   >  7:241922/1‑147 (MQ=255)                                                                                                                gCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCAGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAg                                <  2:139815/147‑1 (MQ=255)                                                                                                                  ccGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGcc                             >  5:275913/1‑148 (MQ=255)                                                                                                                   cGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGa                           >  4:211971/1‑149 (MQ=255)                                                                                                                            gCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGc                  <  8:241922/149‑1 (MQ=255)                                                                                                                              cGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCg                >  7:289854/1‑149 (MQ=255)                                                                                                                                    gcagGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTgcggc                                                >  2:208079/1‑111 (MQ=255)                                                                                                                                    gcagGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTgcggc                                                <  1:208079/111‑1 (MQ=255)                                                                                                                                    gcagGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGt                                             >  2:250957/1‑114 (MQ=255)                                                                                                                                    gcagGAGGCCCAGGCGCGCGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGt                                             <  1:250957/114‑1 (MQ=255)                                                                                                                                         aggCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCGCACGCCTGAAc     <  3:211971/149‑1 (MQ=255)                                                                                                                                            cccAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCGCACGCCTGAACCTg  <  6:275913/149‑1 (MQ=255)                                                                                                                                            cccAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCGCACGCCTGAACCTg  <  5:202421/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            GTGTTGCACAAGCGCTTCGAGGAAGGGGCTGACGTCAAGGCGGGGGATGTGCTGTTCCAGATCGACCCGGCGCCGTTCAAGGCCGCACTGGCCCGTGCCGAGGCTGACCTGGCCCGAGCCCAGGCCGTGCAGCAGGAGGCCCAGGCGCGGGTGAAGCGTTATGAGCCACTGGTCAAGATCGAGGCGGTCAGCCAGCAGGACTTCGACAGCGCCACCGCCGAGCTGCGCAGTGCCGGTGCGGCGGTGCGCTCTGCCCAGGCCGACGTGCAGGCCGCACGCCTGAACCTG  >  NC_002947/3213153‑3213440

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