Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,417,332:1	+A	13.0%	coding (580/759 nt)	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,417,332	1	.	A	13.0%	178.4 / 22.5	54	coding (580/759 nt)	PP_1239	metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base . (25/22);  new base A (3/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

GCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCA  >  NC_002947/1417190‑1417477                                                                                                                                                 |                                                                                                                                                  gCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGCCATTTGAACAACCACCAGGCTGCGC‑GCCTg                                                                                                                                              <  1:28702/148‑1 (MQ=255)   tGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGACGCTAACCACTGTCGCGACCTGCTGACACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTg                                                                                                                                           >  8:313579/1‑149 (MQ=255)     gAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGC‑GCGCAGCCTGGTgg                                                                                                                                          <  6:259189/148‑2 (MQ=255)     gAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGC‑GCGCAGCCTGGTg                                                                                                                                           <  1:56763/147‑1 (MQ=255)              cAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTggg                                                                                                                                >  5:128380/1‑149 (MQ=255)                 ggCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTg                                                                                                                             >  6:27982/1‑149 (MQ=255)                     tGGATGCACTGCTGATCGAGGCTACCCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGa                                                                                                                                      <  5:280017/136‑1 (MQ=255)                     tGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCGGGTGCACGa                                                                                                                                      >  6:280017/1‑136 (MQ=255)                       gATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGaacaa                                                                                                                       <  3:88723/149‑1 (MQ=255)                            aCTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTg                                                                                                                   >  5:84864/1‑148 (MQ=255)                              tgctgATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCaa                                                                                                                >  1:39587/1‑149 (MQ=255)                                  gATCGAGGCTAACCACTGTCGCGTCCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGGATGCTAGGGCATTTTAACAACCACCAGGCTGCGC‑GCCTGGTGCTCGAGTTGGGCTgg                                                                                                                            <  7:145620/133‑1 (MQ=255)                                  gATCGAGGATAAACACTGTCGCGAACTGCTGGCACGTGGTCACTACCAGTATTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAATGGCATTTGAACAACCACCAGGCTGCGA‑GCCTGGTTCACGATTTggtcagg                                                                                                                            >  8:145620/1‑128 (MQ=255)                                   aTCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCTTGGTGCACGAGTTGGGCTGGAACAACCTGCAacac                                                                                                             <  1:235342/147‑1 (MQ=255)                                   aTCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTgg                                                                                                                                 <  6:319109/127‑1 (MQ=255)                                   aTCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTgg                                                                                                                                 >  5:319109/1‑127 (MQ=255)                                    tCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTg                                                                                                          <  6:154460/149‑1 (MQ=255)                                       aGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTg                                                                                                       <  6:128380/149‑1 (MQ=255)                                          cTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCAGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTg                                                                                                    >  5:332878/1‑149 (MQ=255)                                           tAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAg                                                                                                                                     >  6:177591/1‑115 (MQ=255)                                           tAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAg                                                                                                                                     <  5:177591/115‑1 (MQ=255)                                           tAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTgg                                                                                                   >  4:250962/1‑149 (MQ=255)                                                    tCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTca                                                                                          <  8:202335/149‑1 (MQ=255)                                                     cgcgACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTcag                                                                                         >  7:118494/1‑149 (MQ=255)                                                       cgACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTcagca                                                                                       <  4:64190/149‑1 (MQ=255)                                                                 gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCAcc                                                                                             >  6:12110/1‑133 (MQ=255)                                                                 gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCAcc                                                                                             <  5:12110/133‑1 (MQ=255)                                                                 gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaacaa                                                                             >  1:270899/1‑149 (MQ=255)                                                                  cACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaacaac                                                                            <  2:39587/149‑1 (MQ=255)                                                                      tGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGc                                                                        <  3:250962/149‑1 (MQ=255)                                                                        gTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCca                                                                      >  2:253868/1‑149 (MQ=255)                                                                         tCACTACCCGGCTTTTCTGAAGCAGCGGAGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCt                                                                          >  1:21047/1‑145 (MQ=255)                                                                         tCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCt                                                                          <  2:21047/144‑1 (MQ=255)                                                                             tACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACtt                                                                 >  3:22955/1‑149 (MQ=255)                                                                                       ttCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGcc                                                         <  8:118494/147‑1 (MQ=255)                                                                                          tGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGc                                                    >  8:132355/1‑149 (MQ=255)                                                                                                             gCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTg                                                <  6:157602/135‑1 (MQ=255)                                                                                                             gCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTg                                                >  5:157602/1‑135 (MQ=255)                                                                                                              cATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGc                                                              <  8:177683/120‑1 (MQ=255)                                                                                                              cATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGc                                                              >  7:177683/1‑120 (MQ=255)                                                                                                                tGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACCCCTGGTGCTGGCCCACCTCAGCAGAAAGAACAACCTGACACACTTGGCCAGCAAGTGCTTTGTCGACACCCTccgatc                                >  5:366050/1‑143 (MQ=255)                                                                                                                   aaGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCCCCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGAc                            >  3:212735/1‑149 (MQ=255)                                                                                                                               acaacCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCgg                                   >  7:118260/1‑130 (MQ=255)                                                                                                                               acaacCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCgg                                   <  8:118260/130‑1 (MQ=255)                                                                                                                                 gtgcaccaGGC‑GCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTc                >  2:356289/4‑147 (MQ=255)                                                                                                                                 aacCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGaacaac                                                                                                                      <  4:113847/45‑1 (MQ=255)                                                                                                                                 aacCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGaacaac                                                                                                                      >  3:113847/1‑45 (MQ=255)                                                                                                                                  tgcaccaGGC‑GCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCg                         <  6:274691/135‑1 (MQ=255)                                                                                                                                  tgcaccaGGC‑GCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCg                         >  5:274691/3‑137 (MQ=255)                                                                                                                                  tgcaccaGGC‑GCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAg             >  7:78430/3‑149 (MQ=255)                                                                                                                                  accacgAGGCTGCGC‑GCCTGGTGGACGAGTTGGGCTGGAACCACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAg             <  4:22955/149‑1 (MQ=255)                                                                                                                                     accaGGC‑GCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTg          <  1:356289/149‑1 (MQ=255)                                                                                                                                      gcaGGCTGCTC‑GCGTTTTGCACGCTTTGTGCTGGAACATCCTGAAACTCTTGGTGCTTTACCACATCAGAAGCTAGTACAACATGCCACTCTTGGCCCGCCTGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTgg         <  2:270899/148‑1 (MQ=255)                                                                                                                                          gCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGaa     >  2:287024/1‑149 (MQ=255)                                                                                                                                              cgc‑gcCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCa  <  6:84864/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGG‑GTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGC‑GCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCA  >  NC_002947/1417190‑1417477

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