Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,882,291:1	+G	13.3%	coding (906/1416 nt)	oprN →	multidrug RND transporter outer membrane protein OprN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,882,291	1	.	G	13.3%	100.2 / 10.9	30	coding (906/1416 nt)	oprN	multidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base . (14/12);  new base G (2/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCG  >  NC_002947/3882145‑3882439                                                                                                                                                    |                                                                                                                                                     gACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCCGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTg                                                                                                                                                   <  2:334419/150‑1 (MQ=255)      gcTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCAc                                                                                                                                              <  1:692987/150‑1 (MQ=255)          aGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCAccgcc                                                                                                                                          >  1:12930/1‑150 (MQ=255)             gTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGAc                                                                                                                                       >  1:531788/1‑150 (MQ=255)                   cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCA‑CGCCAACGTTGGCGGTGGccgc                                                                                                                                              >  1:790880/1‑135 (MQ=39)                   cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTc                                                                                                                                 <  1:519074/150‑1 (MQ=255)                   cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCACCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCA‑CGCCAACGTTGGCGGTGGccgc                                                                                                                                              <  2:790880/137‑3 (MQ=39)                          ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCcgcgcg                                                                                                                          >  1:285365/1‑150 (MQ=255)                           cAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGt                                                                                                                         >  1:21338/1‑150 (MQ=255)                                     ccGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                               >  1:445745/1‑150 (MQ=255)                                        gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGc                                                                                                                                                 <  2:647823/113‑1 (MQ=255)                                        gCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGc                                                                                                                                                 >  1:647823/1‑113 (MQ=255)                                                            ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                                        >  2:828002/1‑150 (MQ=255)                                                                      ccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCg                                                                              >  2:690914/1‑150 (MQ=255)                                                                                 gctgcGCCGCCGCCCGGACATCCGCAGTGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGctc                                                                   >  2:370569/1‑150 (MQ=255)                                                                                 gctgcGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGctc                                                                   >  2:410813/1‑150 (MQ=255)                                                                                           cgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGcc                                                         <  2:445745/150‑1 (MQ=255)                                                                                              ccGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAAc                                                      >  1:51333/1‑150 (MQ=255)                                                                                                    aTCCGCAGTGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                                <  1:370569/150‑1 (MQ=255)                                                                                                     tCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                               <  1:452869/150‑1 (MQ=255)                                                                                                     tCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                               >  2:77127/1‑150 (MQ=255)                                                                                                       cGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCAGATCTTCGGCGGCCAAAGCCTAGGCa                                             >  1:601890/1‑150 (MQ=255)                                                                                                         cAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACt                                           <  2:242306/150‑1 (MQ=255)                                                                                                                 aGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGa                                   <  1:77127/150‑1 (MQ=255)                                                                                                                           ggcggtggcCA‑CGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTg                         >  2:598758/1‑150 (MQ=39)                                                                                                                            gcggtggcCA‑CGCCAACGTTGGCGGTGGCCGCCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTgg                        <  1:598758/150‑1 (MQ=39)                                                                                                                              ggcggcCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGCCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGggc                      <  1:690914/150‑1 (MQ=255)                                                                                                                                    cACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCTAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCtt                <  1:410813/150‑1 (MQ=255)                                                                                                                                        gCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGAc            <  2:531788/150‑1 (MQ=255)                                                                                                                                                  gC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTATCCAACCGCCCGCGGCTCAAAAAGCGGCGCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCg  >  1:743589/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGC‑GTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCGACCTGGGCAGCG  >  NC_002947/3882145‑3882439

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