Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,416,951	T→A	7.4%	L67M (TTG→ATG)	PP_1239 →	metallo‑beta‑lactamase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,416,951	0	T	A	7.4%	162.1 / 3.4	54	L67M (TTG→ATG)	PP_1239	metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (2/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.53e-01

GTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCGCGGGTCGAG  >  NC_002947/1416811‑1417097                                                                                                                                              |                                                                                                                                                   gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCa                                                                                                                                            >  3:114026/1‑149 (MQ=255)   gACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAAGCCGACCCCATGCACGTG‑GTCGGCATgttcgtgg                                                                                                                                           >  8:50348/1‑141 (MQ=38)   gACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAAGCCGACCCCATGCACGTG‑GTCGGCATgttcgtgg                                                                                                                                           <  7:50348/148‑8 (MQ=38)      acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACggtg                                                                                                                                       >  5:279342/1‑147 (MQ=255)          tCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTa                                                                                                                                    >  3:35867/1‑148 (MQ=255)             tCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAAt                                                                                                                                >  7:211612/1‑149 (MQ=255)             tCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAAt                                                                                                                                <  2:272321/149‑1 (MQ=255)                tGGTCGATTGCCGCTTATCCCTGCCTGAAACCGAGCGGCGCCTGGCTCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTTCAGTGCTGGTCACCCACGAACATGCCGACCAAGAGCATGGG‑GTCGGCTTGCTGTCACGGAGCTACAAtct                                                                                                                              >  8:233096/1‑146 (MQ=255)                  gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTAc                                                                                                                            <  8:211612/148‑1 (MQ=255)                          cGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAcc                                                                                                                   >  8:137350/1‑149 (MQ=255)                             cTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCa                                                                                                                <  4:114026/149‑1 (MQ=255)                             cTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCa                                                                                                                <  3:52856/149‑1 (MQ=255)                                   ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          >  1:231171/1‑133 (MQ=255)                                   ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          <  2:231171/133‑1 (MQ=255)                                   ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                                                                                                          <  3:300098/149‑1 (MQ=255)                                    cTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAgg                                                                                                         <  6:279342/149‑1 (MQ=255)                                    cTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAgg                                                                                                         >  6:185434/1‑149 (MQ=255)                                        gTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACt                                                                                                     >  5:22501/1‑149 (MQ=255)                                           aaaCCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTg                                                                                                  >  6:172390/1‑149 (MQ=255)                                               cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGCGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          <  1:42323/122‑1 (MQ=255)                                               cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGAGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          <  1:37500/122‑1 (MQ=255)                                               cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          >  2:42323/1‑121 (MQ=255)                                               cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                                          >  2:37500/1‑121 (MQ=255)                                                   cggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCAt                                                                                          <  5:185434/149‑1 (MQ=255)                                                     gcgcCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGa                                                                                                       >  6:200485/1‑134 (MQ=255)                                                     gcgcCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGa                                                                                                       <  5:200485/134‑1 (MQ=255)                                                        ccTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATgcgc                                                                                      >  4:73657/1‑148 (MQ=255)                                                         cTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCaa                                                                                    <  4:35867/149‑1 (MQ=255)                                                                tgctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTg                                                                             <  2:268184/149‑1 (MQ=255)                                                                  ctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTg                                                                             >  5:268455/1‑147 (MQ=255)                                                                     ggggTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAgg                                                                         <  6:268455/148‑1 (MQ=255)                                                                              gCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAg                                                                          <  2:139690/138‑1 (MQ=255)                                                                              gCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAg                                                                          >  1:139690/1‑138 (MQ=255)                                                                               cagctcagctCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGtt                                                                <  3:73657/147‑1 (MQ=255)                                                                                agctcagctcAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTtc                                                             >  7:242874/1‑149 (MQ=255)                                                                                     agctcagTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTcgc                                                         <  8:242874/148‑1 (MQ=255)                                                                                        tcagTGCAGTGCTGGTCAACCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTATCGCgtgt                                                     >  3:16580/1‑149 (MQ=255)                                                                                           gtgcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGc                                                  >  5:100558/1‑149 (MQ=255)                                                                                            tgcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGc                                                  >  3:234304/1‑148 (MQ=255)                                                                                                      gTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTg                                         <  5:172390/147‑1 (MQ=255)                                                                                                         aCCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCATCCAAGGGACTTTGCGCTGAATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCata                                     <  7:233096/148‑1 (MQ=255)                                                                                                           ccACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGcc                                                 <  4:152204/134‑1 (MQ=255)                                                                                                           ccACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGcc                                                 >  3:152204/1‑134 (MQ=255)                                                                                                            cACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCAGGTGGAGCTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATAACgg                                 >  7:108784/1‑149 (MQ=255)                                                                                                             aCGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGTCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCTTGTGGCCAAAGACTGCATATCgg                                 >  5:64146/1‑148 (MQ=255)                                                                                                                 gcaagCCGACCCCATGCACGTG‑GTCGGCATGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGa                                                                           <  7:26834/98‑1 (MQ=37)                                                                                                                 gcaagCCGACCCCATGCACGTG‑GTCGGCATGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGa                                                                           >  8:26834/5‑102 (MQ=37)                                                                                                                  cATGCCGGCCACGTGCGTGGG‑GTCGGCTTGCTGTCACGGCGCTACAGTGTTGCGGTCTAGCTCAGCCAAGGGGCTTTGCGCGGCATGGGCGGGGCGGTGGAGGTGGCCGGGTTTCTCGGGTGTGGCCAGCGCGTGCAGATCggggggct                           >  8:257816/1‑143 (MQ=255)                                                                                                                      ccGTCCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGaa                       <  4:179112/149‑1 (MQ=255)                                                                                                                      ccGACCCCGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGaa                       <  2:178072/149‑1 (MQ=255)                                                                                                                      ccGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGaa                       <  7:137350/149‑1 (MQ=255)                                                                                                                          ccACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTaa                   >  7:54120/1‑149 (MQ=255)                                                                                                                                      g‑gTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCGCggg       >  3:236159/1‑149 (MQ=255)                                                                                                                                            gCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCGCGGGTCGAg  <  4:234304/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   GTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGG‑GTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAGTAACCGCAGCGCGGGTCGAG  >  NC_002947/1416811‑1417097

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