Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,882,813	A→G	9.2%	intergenic (+56/‑169)	oprN → / → PP_3428	multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,882,813	0	A	G	9.2%	145.4 / 5.8	54	intergenic (+56/‑169)	oprN/PP_3428	multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base G (3/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGC  >  NC_002947/3882669‑3882956                                                                                                                                                  |                                                                                                                                                gTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGc                                                                                                                                              >  3:178809/1‑148 (MQ=255) gTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGcc                                                                                                                                             <  6:214630/149‑1 (MQ=255)         gACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCGTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGccc                                                                                                                                      >  1:63678/1‑149 (MQ=255)          aCCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCTCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCagc                                                                                                                                    >  4:7381/1‑149 (MQ=255)            cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgg                                                                                                                                  <  3:46292/149‑1 (MQ=255)               tATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                 <  7:206655/147‑1 (MQ=255)                   gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                           >  8:98922/1‑149 (MQ=255)                   gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                           <  4:139966/149‑1 (MQ=255)                   gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtgtt                                                                                                                           >  5:227771/1‑149 (MQ=255)                     ggCATCGTCTCGATCTACAAGGCCCTGGGTCGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCACCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCttttt                                                                                                                          >  1:24873/1‑148 (MQ=255)                           gTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                   <  7:98922/149‑1 (MQ=255)                             ctcGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGc                                                                                                                 >  1:88915/1‑149 (MQ=255)                                     aCAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                   >  8:159793/1‑139 (MQ=255)                                     aCAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                   <  7:159793/139‑1 (MQ=255)                                      cAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATc                                                                                                             <  5:30177/144‑1 (MQ=255)                                           cccTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                   <  3:7381/149‑1 (MQ=255)                                           cccTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                   <  2:63678/149‑1 (MQ=255)                                            ccTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgggg                                                                                                                                <  8:39993/119‑1 (MQ=255)                                            ccTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgggg                                                                                                                                >  7:39993/1‑119 (MQ=255)                                                        c‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGc                                                                                      <  2:88915/149‑1 (MQ=255)                                                            gCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCa                                                                                          >  8:192573/1‑142 (MQ=255)                                                            gCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCa                                                                                          <  7:192573/142‑1 (MQ=255)                                                              aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                           >  7:145861/1‑107 (MQ=255)                                                              aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                           <  8:145861/107‑1 (MQ=255)                                                               aCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGaa                                                                                >  7:177185/1‑149 (MQ=255)                                                               aCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGaa                                                                                <  6:97986/149‑1 (MQ=255)                                                                        ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGc                                                                                                                               <  5:62349/93‑1 (MQ=255)                                                                        ggTGGCCAGCGCGCACTGATTACACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGc                                                                                                                               >  6:62349/1‑93 (MQ=255)                                                                              cAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTc                                                                                                       <  3:158888/111‑1 (MQ=255)                                                                              cAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTc                                                                                                       >  4:158888/1‑111 (MQ=255)                                                                                    gcACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcc                                                           >  3:91055/1‑149 (MQ=255)                                                                                      aCTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCAt                                                                                                                                           <  3:161413/67‑1 (MQ=255)                                                                                      aCTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCAt                                                                                                                                           >  4:161413/1‑67 (MQ=255)                                                                                            tCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAAAGCCACCGCAGGCCTGGCGCCTGGCCAAc                                                   >  6:102941/1‑149 (MQ=255)                                                                                             ccACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACACCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACg                                                  >  3:146415/1‑149 (MQ=255)                                                                                              cACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACgg                                                 >  4:231046/1‑149 (MQ=255)                                                                                                   ccTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGc                                            >  5:117830/1‑149 (MQ=255)                                                                                                   ccTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGc                                            >  4:86135/1‑149 (MQ=255)                                                                                                     tCAACGACTCCTTTGGTTGGCTCCTACCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATg                                                                                   >  8:8809/1‑108 (MQ=255)                                                                                                     tCAACGACTCCTTTGGTTGGCTCCTACCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATg                                                                                   <  7:8809/108‑1 (MQ=255)                                                                                                         cGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                      <  5:15184/149‑1 (MQ=255)                                                                                                           aCTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcg                                    <  3:231046/149‑1 (MQ=255)                                                                                                             tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGACGTGCCcggttcggt                                  >  4:243828/1‑149 (MQ=255)                                                                                                               cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGt                                                                                        <  6:210824/93‑1 (MQ=255)                                                                                                               cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGt                                                                                        >  5:210824/1‑93 (MQ=255)                                                                                                               cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTg                                >  8:223040/1‑149 (MQ=255)                                                                                                               cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGACGTGCCCGGTTCGGTTg                                <  3:243828/149‑1 (MQ=255)                                                                                                                     ttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTAATCACCCCAATACCCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACCCCGTGCCCGGTTCGGTTGCAttt                           >  6:191796/1‑148 (MQ=255)                                                                                                                               ccccAACCGTTTGCCCCGCATGCCATTGCCCGCGGGTCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGGCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCt                <  6:227771/149‑1 (MQ=255)                                                                                                                                           gCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGctc     <  2:24873/148‑1 (MQ=255)                                                                                                                                            ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGctct    <  4:146415/148‑1 (MQ=255)                                                                                                                                            ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg   >  4:83179/1‑149 (MQ=255)                                                                                                                                            ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg   <  4:91055/149‑1 (MQ=255)                                                                                                                                             cccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc  >  3:205416/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                GTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGC‑TGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGC  >  NC_002947/3882669‑3882956

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