Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,840,662	G→A	7.1%	T286M (ACG→ATG)	PP_3390 ←	porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,840,662	0	G	A	7.1%	162.7 / 2.5	56	T286M (ACG→ATG)	PP_3390	porin
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTTGTAGTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCGGGTCTTGAACTTAA  >  NC_002947/3840514‑3840806                                                                                                                                                     |                                                                                                                                                  aaaGTTGTAGTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACg                                                                                                                                                   >  3:298355/1‑149 (MQ=255)     ttGTAGTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGtgtg                                                                                                                                               <  2:95482/149‑1 (MQ=255)       gTAGTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTgg                                                                                                                                              <  3:126682/148‑1 (MQ=255)          gTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCt                                                                                                                                          <  6:277563/149‑1 (MQ=255)             gTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGa                                                                                                                                        <  2:10905/148‑1 (MQ=255)                       aTCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTAccc                                                                                                                             >  1:180657/1‑149 (MQ=255)                        tCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTAccc                                                                                                                             >  4:254618/1‑148 (MQ=255)                          cACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGaa                                                                                                                          <  4:298355/149‑1 (MQ=255)                           aCGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGCACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGCACGTGGGGCCCTGAACGCTGTACCCGAAg                                                                                                                         >  5:289551/1‑149 (MQ=255)                           aCGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAg                                                                                                                         >  6:322173/1‑149 (MQ=255)                               ccGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTg                                                                                                                      <  4:56876/148‑1 (MQ=255)                               ccGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTg                                                                                                                      >  3:56876/1‑148 (MQ=255)                                     gttGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGa                                                                                                               <  3:254618/149‑1 (MQ=255)                                             aGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGa                                                                                                       <  5:322173/149‑1 (MQ=255)                                              ggCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGAt                                                                                                      >  2:300419/1‑149 (MQ=255)                                              ggCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGAt                                                                                                      >  3:336920/1‑149 (MQ=255)                                               gCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGAtt                                                                                                     <  1:118047/149‑1 (MQ=255)                                                 aaaTGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGt                                                                                                   >  3:121798/1‑149 (MQ=255)                                                  aaTGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTc                                                                                                  <  4:336920/149‑1 (MQ=255)                                                          gACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACt                                                                                                                 <  3:195524/126‑1 (MQ=255)                                                          gACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACt                                                                                                                 >  4:195524/1‑126 (MQ=255)                                                                   ttCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAAc                                                                                 <  1:300419/148‑1 (MQ=255)                                                                   ttCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACg                                                                                >  8:22713/1‑149 (MQ=255)                                                                   ttCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACATGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACg                                                                                >  4:177008/1‑149 (MQ=255)                                                                     cGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGcc                                                                              >  7:100361/1‑149 (MQ=255)                                                                      gCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGAtt                                                                                                     <  7:23830/126‑1 (MQ=255)                                                                      gCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGAtt                                                                                                     >  8:23830/1‑126 (MQ=255)                                                                          aaCAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGcc                                                                         >  3:336682/1‑149 (MQ=255)                                                                          aaCAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGcc                                                                         <  8:100361/149‑1 (MQ=255)                                                                          aaCAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGc                                                                          <  4:88666/148‑1 (MQ=255)                                                                               gtaggtaTCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCt                                                                    >  6:155755/1‑149 (MQ=255)                                                                                taggtaTCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTc                                                                   >  3:54269/1‑149 (MQ=255)                                                                                     aTCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACATGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTc                                                                                                  <  6:120135/114‑1 (MQ=255)                                                                                     aTCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACATGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTc                                                                                                  >  5:120135/1‑114 (MQ=255)                                                                                        ggTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCa                                                            <  4:54269/148‑1 (MQ=255)                                                                                         gTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCAtt                                                          >  6:147691/1‑149 (MQ=255)                                                                                            ccccccACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGaa                                                       <  2:180657/149‑1 (MQ=255)                                                                                              ccccACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAAc                                                      >  1:18336/1‑148 (MQ=255)                                                                                                           tatGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGTACGTGTGGCCCTGAACGCTGTACCCGTATTTGGTACTGAGGGTCTGATTGTCAGACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCgg                                         <  1:281458/149‑1 (MQ=255)                                                                                                           tatGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGt                                        <  1:13366/149‑1 (MQ=255)                                                                                                                    tCCCCACGCACTTTCTGATAGCCGCCACTGAACATGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATc                               <  3:177008/149‑1 (MQ=255)                                                                                                                      cccACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGc                                                                     <  8:263042/109‑1 (MQ=255)                                                                                                                      cccACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGc                                                                     >  7:263042/1‑109 (MQ=255)                                                                                                                      cccACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATccc                             <  2:18336/149‑1 (MQ=255)                                                                                                                        cACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATccccc                           >  6:107368/1‑149 (MQ=255)                                                                                                                         aCGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCgg                                         >  8:332916/1‑134 (MQ=255)                                                                                                                         aCGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCgg                                         <  7:332916/134‑1 (MQ=255)                                                                                                                             aCTTTCTGATAGCCGCCACTGAACGTGTGGCACTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCTCCCAGCAACGCCCGGCCCTGCTCACTTGCATTGAACAAGAGCCACTCGTTGAACACATCCCCCGccc                       >  4:287007/1‑148 (MQ=255)                                                                                                                               tttCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGccc                      >  3:223890/1‑147 (MQ=255)                                                                                                                                  cTGTTAGCCGCGACTGAAAGTGTGGCCCTGAACGCTGTACCGGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTGCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCTTTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGt                  <  3:287007/148‑1 (MQ=255)                                                                                                                                     aTAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCggg              >  1:335337/1‑149 (MQ=255)                                                                                                                                     aTAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCgg               <  7:22713/148‑1 (MQ=255)                                                                                                                                      tAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCgggg             >  3:213940/1‑148 (MQ=255)                                                                                                                                         ccgccACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCGGGTCtt          <  5:107368/149‑1 (MQ=255)                                                                                                                                                  aaCGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGa                                                        <  3:50836/94‑1 (MQ=255)                                                                                                                                                  aaCGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGa                                                        >  4:50836/1‑94 (MQ=255)                                                                                                                                                  aaCGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCGGGTCTTGAACTTaa  <  5:147691/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAAGTTGTAGTCGTAGCGCAGCATCCACGCCCGTTCGTTGGCCAAGGCAAATGTGCTGACCTGCTGTTCGCTGAACAGGTAGGTATCGGTGCCCCCCACGTAGGCATATGCCGTATCCCCACGCACTTTCTGATAGCCGCCACTGAACGTGTGGCCCTGAACGCTGTACCCGAAGTTGGTACTGAGGGTCTGATTGTC‑GACTCCGCCCAGCAACGCCCGGCCCTGCTCCCTTGCATTGAACCAGCGCCACTCGGTGAACACATCCCCCGGCCCAAGTCGGGTCTTGAACTTAA  >  NC_002947/3840514‑3840806

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