Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	325,669	A→G	7.7%	T319A (ACC→GCC)	oprQ →	outer‑membrane porin D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	325,669	0	A	G	7.7%	143.9 / 2.9	52	T319A (ACC→GCC)	oprQ	outer‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (2/2);  total (26/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 = 9.87e-01

GCGACTCCAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGTGGTGGGTGACAACATCAACA  >  NC_002947/325523‑325816                                                                                                                                                    |                                                                                                                                                       gCGACTCCAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTccacc                                                                                                                                                      <  7:202506/149‑1 (MQ=255) gCGACTCCAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTccac                                                                                                                                                       >  8:63436/1‑148 (MQ=255)      tCCAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTCAAGCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑Tc                                                                                                                                                  >  1:110894/1‑148 (MQ=255)        cAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTa                                                                                                                                                <  7:226240/147‑1 (MQ=255)         aaCACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCt                                                                                                                                             >  2:87343/1‑149 (MQ=255)             ccATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGcc                                                                                                                                         <  8:214478/149‑1 (MQ=255)               aTCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCt                                                                                                                                                 >  5:120135/1‑139 (MQ=255)               aTCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCt                                                                                                                                                 <  6:120135/139‑1 (MQ=255)                      gTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTAc                                                                                                                                <  7:63436/149‑1 (MQ=255)                         ttGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTgg                                                                                                                             >  8:343331/1‑149 (MQ=255)                             aGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTcc                                                                                                                          >  2:267452/1‑148 (MQ=255)                                   cTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTc                                                                                                                    >  3:267619/1‑148 (MQ=255)                                        cgTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCGCCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGcc                                                                                                              >  6:272817/1‑149 (MQ=255)                                        cgTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCGCCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGcc                                                                                                              >  6:267750/1‑149 (MQ=255)                                         gTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCGCCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCg                                                                                                             <  5:272817/149‑1 (MQ=255)                                         gTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCGCCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCg                                                                                                             <  5:267750/149‑1 (MQ=255)                                                 cgcACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCATTCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGa                                                                                                      <  1:267452/149‑1 (MQ=255)                                                 cgcACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAg                                                                                                     <  7:343331/149‑1 (MQ=255)                                                        gTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCt                                                                                              <  2:110894/149‑1 (MQ=255)                                                            aCCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGc                                                                                           >  4:29181/1‑148 (MQ=255)                                                             cccTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCaa                                                                                         <  1:174712/149‑1 (MQ=255)                                                              ccTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAAc                                                                                        >  1:22855/1‑149 (MQ=255)                                                              ccTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAAc                                                                                        >  3:11752/1‑149 (MQ=255)                                                                  gCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTggg                                                                                    <  1:133971/149‑1 (MQ=255)                                                                  gCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTg                                                                                      <  4:158764/147‑1 (MQ=255)                                                                    tC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGc                                                                                   >  2:270157/1‑148 (MQ=255)                                                                    tC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGc                                                                                   >  1:318221/1‑148 (MQ=255)                                                                         cAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑Cgg                                                                             >  7:117831/1‑149 (MQ=255)                                                                             gcagcaCTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTg                                                                         <  4:11752/149‑1 (MQ=255)                                                                               agcaCTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGa                                                                       >  8:38985/1‑149 (MQ=255)                                                                                      gcagcaCGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCg                                                                  <  2:230251/147‑1 (MQ=255)                                                                                        agcaCGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGc                                                              <  3:29181/149‑1 (MQ=255)                                                                                        agcaCGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGCTGc                                                              >  2:320912/1‑149 (MQ=255)                                                                                          caCGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGcc                                                             >  8:223134/1‑148 (MQ=255)                                                                                          caCGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGcc                                                             >  4:343491/1‑148 (MQ=255)                                                                                                 tACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCAt                                                     >  4:38300/1‑149 (MQ=255)                                                                                                   cAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTg                                                                         <  5:199237/127‑1 (MQ=255)                                                                                                   cAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTg                                                                         >  6:199237/1‑127 (MQ=255)                                                                                                      cTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCgg                                                <  1:270157/149‑1 (MQ=255)                                                                                                            cTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGa                                         <  8:117831/149‑1 (MQ=255)                                                                                                             ttCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGAc                                                                                                        >  2:131992/1‑87 (MQ=255)                                                                                                             ttCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGAc                                                                                                        <  1:131992/87‑1 (MQ=255)                                                                                                                    cccaCGCCGCCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAg                                 >  5:313682/4‑149 (MQ=255)                                                                                                                     gAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCT‑GCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCt                               >  2:221147/1‑149 (MQ=255)                                                                                                                      aaCGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACt                                                                     <  5:143035/112‑1 (MQ=255)                                                                                                                      aaCGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACt                                                                     >  6:143035/1‑112 (MQ=255)                                                                                                                         gCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTACCGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTgg                             <  1:87343/149‑1 (MQ=255)                                                                                                                             gcggcgTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGt                                                                <  1:122264/110‑1 (MQ=255)                                                                                                                             gcggcgTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGt                                                                >  2:122264/1‑110 (MQ=255)                                                                                                                               ggcgTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACgt                      >  5:171988/1‑149 (MQ=255)                                                                                                                                gcgTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACgtg                     <  7:38985/149‑1 (MQ=255)                                                                                                                                   tGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACgtggtg                  <  7:223134/149‑1 (MQ=255)                                                                                                                                      gTGACGGTGGTTCCACCA‑TCTACCT‑GCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGCTCTGACCTACAAGCTGGCTTACGTggtgggt               <  1:221147/149‑1 (MQ=255)                                                                                                                                            gtggtTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGTGGTGGGTGacaac         <  3:343491/149‑1 (MQ=255)                                                                                                                                                   cacca‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGtggtgg                 >  1:330210/1‑134 (MQ=255)                                                                                                                                                   cacca‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGtggtgg                 <  2:330210/134‑1 (MQ=255)                                                                                                                                                   cacca‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGTGGTGGGTGACAACATCAaca  >  2:95004/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GCGACTCCAACACCATCTGGAGTCTTGCAGCCACCTACGCGTACGGTGCGCACTCGTTCACCCTCGCTC‑ACCAGCGCAGCACTGGCAGCACGGGCTACAACTACGGCTTCTACCAGAACGCCGGCGGCGTGGGTGACGGTGGTTCCACCA‑TCTACCTGGCCAACTCGTACTGGTCCGACTTCAACGCCGAAGACGAGCGCTCCTGGCAACTGGGCTA‑CGGCCTGGACTTCGGTGCCTACGGCATTCCGGGC‑CTGACCTACAAGCTGGCTTACGTGGTGGGTGACAACATCAACA  >  NC_002947/325523‑325816

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