Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,141,771	T→A	8.9%	M245K (ATG→AAG) ‡	PP_3645 →	L‑proline amide hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,141,771	0	T	A	8.9%	132.0 / 2.8	45	M245K (ATG→AAG) ‡	PP_3645	L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (2/2);  total (24/21)
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

ACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCC  >  NC_002947/4141624‑4141919                                                                                                                                                    |                                                                                                                                                     aCCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCatat                                                                                                                                                      >  1:299701/1‑148 (MQ=255)             acCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGaa                                                                                                                                         <  5:58572/149‑1 (MQ=255)                       tATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGCCCCGCCAGCTCGCCGCCAGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGt                                                                                                                                >  6:317651/1‑148 (MQ=255)                         tCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGaa                                                                                                                             <  7:159707/149‑1 (MQ=255)                           cAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGt                                                                                                                           <  5:265010/149‑1 (MQ=255)                            aGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGt                                                                                                                           >  3:12467/1‑148 (MQ=255)                              gTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATCATATGCCGTTgct                                                                                                                                          <  1:209568/131‑3 (MQ=255)                              gTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATCATATGCCGTTgct                                                                                                                                          >  2:209568/1‑129 (MQ=255)                                 ccGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTcaac                                                                                                                       >  6:116079/1‑147 (MQ=255)                                   gACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTc                                                                                                                          >  7:294784/1‑142 (MQ=255)                                   gACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTc                                                                                                                          <  8:294784/142‑1 (MQ=255)                                    aCGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGAt                                                                                                                  >  5:42763/1‑149 (MQ=255)                                        tGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCACCGGCATCATATGCCGTTgctc                                                                                                                                         <  6:252454/122‑4 (MQ=255)                                        tGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCACCGGCATCATATGCCGTTgctc                                                                                                                                         >  5:252454/1‑119 (MQ=255)                                           tCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                           <  4:101501/149‑1 (MQ=255)                                            cGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                           >  4:130813/1‑148 (MQ=255)                                                   ggCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCCGt                                                                                                                                             >  8:298567/1‑107 (MQ=255)                                                    gCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCt                                                                                                   >  6:16117/1‑148 (MQ=255)                                                    gCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTc                                                                                                  <  4:12467/149‑1 (MQ=255)                                                          tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                           <  3:187585/134‑1 (MQ=255)                                                          tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGctgctg                                                                                                           >  4:187585/1‑134 (MQ=255)                                                          tcgactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCAccc                                                                                              >  6:63810/1‑147 (MQ=255)                                                            gactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGT‑GAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTc                                                                                                     >  7:148033/1‑137 (MQ=255)                                                            gactcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGT‑GAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTc                                                                                                     <  8:148033/137‑1 (MQ=255)                                                              ctcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGccc                                                                                        <  3:130813/149‑1 (MQ=255)                                                               tcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCg                                                                                       >  3:223116/1‑149 (MQ=255)                                                               tcgGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCg                                                                                       >  4:181387/1‑149 (MQ=255)                                                                  ggCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCAc                                                                                    >  6:37650/1‑149 (MQ=255)                                                                   gCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCAcc                                                                                   <  5:116079/149‑1 (MQ=255)                                                                       cACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCa                                                                               <  5:16117/149‑1 (MQ=255)                                                                           ccGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCg                                                                            >  7:186669/1‑148 (MQ=255)                                                                                     cccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAg                                                                 <  5:317651/149‑1 (MQ=255)                                                                                     cccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAg                                                                 <  3:181387/149‑1 (MQ=255)                                                                                      ccGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAgg                                                                >  7:172221/1‑149 (MQ=255)                                                                                                 ccgccTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCa                                                     <  4:172109/149‑1 (MQ=255)                                                                                                             cTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCtcgttc                                         >  6:270566/1‑149 (MQ=255)                                                                                                                                 cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc                     <  5:63810/149‑1 (MQ=255)                                                                                                                                 cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc                     >  4:7781/1‑149 (MQ=255)                                                                                                                                 cTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAc                     >  1:107927/1‑149 (MQ=255)                                                                                                                                  tCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACg                    >  5:299055/1‑149 (MQ=255)                                                                                                                                         gCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATg             <  4:223116/149‑1 (MQ=255)                                                                                                                                          cAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATgg            >  1:103396/1‑149 (MQ=255)                                                                                                                                           aaCGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGgc           <  5:37650/149‑1 (MQ=255)                                                                                                                                             cGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGgcgc         <  6:42763/149‑1 (MQ=255)                                                                                                                                              ggCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGc                                                      >  7:202818/1‑103 (MQ=255)                                                                                                                                              ggCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGc                                                      <  8:202818/103‑1 (MQ=255)                                                                                                                                                    tgatgCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGccc  <  6:299055/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCGGGCTCCACCCCGGCCATGACCCGCCAGCTCGCCGCCTGCATTCCTGGCGCACGCAGTGTGGTCCTCGCCGAGCAACGGCATATGATGCCGGTGGAAGCACCGCGTGAAGTCAACAAGATGCTGCTGGACTTCCTCACCCAAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCC  >  NC_002947/4141624‑4141919

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