Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	336,125	+T	coding (210/237 nt)	PP_0278 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	336,124	1	.	T	97.1%	111.3 / ‑3.3	35	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (18/16);  total (19/16)
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.93e-01

TTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTT  >  NC_002947/335977‑336259                                                                                                                                                      |                                                                                                                                        ttAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑t                                                                                                                                        >  7:330982/1‑149 (MQ=255)      ccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑ttacct                                                                                                                                   >  6:348105/1‑144 (MQ=255)           ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGc                                                                                                                                <  8:209500/148‑1 (MQ=255)            cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATAGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGTAGGACGACCAAAAATTACTGGGGGAAGCCG‑ATTACCTGTGCGa                                                                                                                              >  4:41759/1‑149 (MQ=255)             tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGCT‑TACCTGTGCGa                                                                                                                              >  2:120636/1‑148 (MQ=255)                       aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAgg                                                                                                                   >  5:204840/1‑149 (MQ=255)                               agCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                            >  3:314184/1‑148 (MQ=255)                                 cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGaa                                                                                                         <  8:330982/149‑1 (MQ=255)                                   cGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACg                                                                                                       <  5:3610/149‑1 (MQ=255)                                        tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  <  6:204840/149‑1 (MQ=255)                                             gTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  >  4:239435/1‑144 (MQ=255)                                                 aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGc                                                                                          >  8:12118/1‑148 (MQ=255)                                                  cgcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCa                                                                                        <  4:314184/149‑1 (MQ=255)                                                        cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                       <  3:15232/144‑1 (MQ=255)                                                        cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                       >  4:15232/1‑144 (MQ=255)                                                                           cTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAc                                                               >  5:179848/1‑149 (MQ=255)                                                                               tGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCgg                                                           >  4:243556/1‑149 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  >  8:210053/1‑106 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  <  7:215620/106‑1 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  <  7:210053/106‑1 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                  >  8:215620/1‑106 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcgc                                                                            <  2:200796/128‑1 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCCATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcgc                                                                            >  1:200796/1‑128 (MQ=255)                                                                                         cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCATCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                 <  3:41759/149‑1 (MQ=255)                                                                                           aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGa                                               <  3:239435/149‑1 (MQ=255)                                                                                            gCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAg                                              >  1:125871/1‑149 (MQ=255)                                                                                               ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGAGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCgg                                           >  5:307277/1‑149 (MQ=255)                                                                                                cTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGc                                          >  8:278323/1‑149 (MQ=255)                                                                                                         gAGGGCGATATGGAGGACGACCGTAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                         <  5:80153/125‑1 (MQ=255)                                                                                                         gAGGGCGATATGGAGGACGACCGTAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                         >  6:80153/1‑125 (MQ=255)                                                                                                           gggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgcg                               <  7:12118/149‑1 (MQ=255)                                                                                                               gATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGaagg                           <  5:348105/149‑3 (MQ=255)                                                                                                                 tatGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                         <  3:243556/149‑1 (MQ=255)                                                                                                                     gaggaCGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                      >  4:152085/1‑148 (MQ=255)                                                                                                                        gacgacCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAt                            <  1:41494/139‑1 (MQ=255)                                                                                                                        gacgacCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAt                            >  2:41494/1‑139 (MQ=255)                                                                                                                                         gggggCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt  >  5:358412/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                        TTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTT  >  NC_002947/335977‑336259

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