Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,882,813	A→G	8.7%	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	8.7%	127.1 / 2.8	46	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 (21/21);  new base G (2/2);  total (23/23)
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.36e-01

AGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCC  >  NC_002947/3882666‑3882957                                                                                                                                                     |                                                                                                                                                 aGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAg                                                                                                                                                 >  7:16770/1‑149 (MQ=255)    gTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGc                                                                                                                                               >  2:183371/1‑148 (MQ=255)        aGTCGACCTTTATCGCGGCATCGTCTCGATGTACAAGGCCCTGGGTGGCGGCTGGCAAGCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATg                                                                                                                                           <  8:108890/148‑1 (MQ=255)              ccTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCTGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTTCCCCGCAGGCCATGGCCcgcg                                                                                                                                    <  5:299018/149‑1 (MQ=255)               cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgg                                                                                                                                   >  5:41445/1‑149 (MQ=255)                  tATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgggg                                                                                                                                 <  4:124185/148‑1 (MQ=255)                    tCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCGGGGca                                                                                                                               >  6:285877/1‑147 (MQ=255)                      gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCGGGGcaaa                                                                                                                             <  5:138532/148‑4 (MQ=255)                              gTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                    <  8:16770/149‑1 (MQ=255)                              gTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                    >  3:162041/1‑149 (MQ=255)                                      cTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCg                                                                                                             >  6:295222/1‑148 (MQ=255)                                       tACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTa                                                                                                           >  8:73010/1‑149 (MQ=255)                                         cAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTa                                                                                                           <  6:41445/147‑1 (MQ=255)                                                  gggTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                    <  6:300558/145‑1 (MQ=255)                                                  gggTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                    >  5:300558/1‑145 (MQ=255)                                                      ggcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACg                                                                                             >  6:256339/1‑148 (MQ=255)                                                      gccggcTCGGAAGCGGAGTCGGTGGCCAGCGCGCACTGATTCAACTGCCTCACCGA‑CGCCTTTGGTGGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGCGCTTTTTAATGTTTGCCATCGTATTCACCCCAATACGc                                                                                            <  8:300778/147‑1 (MQ=255)                                                            tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCAGCAGGCCATGGCCCGAGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCa                                                                                      >  7:232968/1‑149 (MQ=255)                                                              gCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATg                                                                                    <  4:97831/149‑1 (MQ=255)                                                               cAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATg                                                                                    >  8:104730/1‑148 (MQ=255)                                                                aaCCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATgg                                                                                   >  4:8888/1‑148 (MQ=255)                                                                  ccGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAAt                                                                                <  7:73010/149‑1 (MQ=255)                                                                    ggAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGc                                                                              <  2:268911/149‑1 (MQ=255)                                                                         cGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCg                                                                         <  7:104730/149‑1 (MQ=255)                                                                                cAGCGCGCACTGATTCCACTGCCTCAACGATCTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcc                                                                   >  1:192773/1‑149 (MQ=255)                                                                                      gcACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcc                                                            <  3:200234/149‑1 (MQ=255)                                                                                         cTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcctgg                                                         <  8:38110/149‑1 (MQ=255)                                                                                         cTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcctgg                                                         <  3:8888/149‑1 (MQ=255)                                                                                              tCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCaa                                                     <  5:295222/148‑1 (MQ=255)                                                                                               ccACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACg                                                   >  1:144456/1‑149 (MQ=255)                                                                                                 ccttccTCAACGA‑CTCCTTTGGGTGGCTCCTCCCCCAACCGTTTGCCCGGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGc                                                 <  8:232968/145‑1 (MQ=255)                                                                                                         aaCGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcgg                                         >  1:142292/1‑149 (MQ=255)                                                                                                           cGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                       <  4:162041/149‑1 (MQ=255)                                                                                                                  cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTg                                 >  8:47904/1‑149 (MQ=255)                                                                                                                  cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGACCGCGGGGCTTTTTTATGTTTGACATCGTATTCACCCAAATACGCAGTGCATGGAATGCCAAAGCAGGCATGGCGCCTGGCAAACGCCGTGCCcggtacggtt                                  >  8:276223/1‑148 (MQ=255)                                                                                                                      ggttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCa                               <  7:43569/147‑1 (MQ=255)                                                                                                                      ggttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCa                               >  8:43569/1‑147 (MQ=255)                                                                                                                              cctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTAcccccc                     >  1:203001/1‑149 (MQ=255)                                                                                                                                  ccccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                       >  6:201156/1‑127 (MQ=255)                                                                                                                                    ccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCt                 <  3:191576/147‑1 (MQ=255)                                                                                                                                    ccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCt                 >  4:191576/1‑147 (MQ=255)                                                                                                                                         cggttGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCaa          <  4:216245/146‑1 (MQ=255)                                                                                                                                               ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCAt                              <  3:51611/123‑1 (MQ=255)                                                                                                                                               ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCAt                              >  4:51611/1‑123 (MQ=255)                                                                                                                                                 ccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCACGCTCTGcc  >  7:151802/1‑149 (MQ=255)                                                                                                                                                 ccGCGGGCAATGGCCTGCGGGGCTTTTTTATGTTTGCAATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGcc  <  1:156443/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGA‑CTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCC  >  NC_002947/3882666‑3882957

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