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.4%	128.6 / ‑3.6	39	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (14/24);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAG  >  NC_002947/335981‑336264                                                                                                                                                  |                                                                                                                                             cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GA‑ttac                                                                                                                                          <  1:35714/148‑4 (MQ=255) cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GA‑ttac                                                                                                                                          <  3:73076/148‑4 (MQ=255)        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGCCCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCg                                                                                                                                    >  1:13084/1‑148 (MQ=255)        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCg                                                                                                                                    >  2:160919/1‑148 (MQ=255)            aTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACAt                                                                                                                                <  6:28521/148‑1 (MQ=255)             tctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGAc                                                                                                                                  >  8:112427/1‑145 (MQ=255)              ctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACa                                                                                                                                 <  7:112427/145‑1 (MQ=255)               tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATgg                                                                                                                              >  5:14746/1‑147 (MQ=255)                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGc                                                                                                                             >  3:142221/1‑146 (MQ=255)                    cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAgg                                                                                                                        >  6:135202/1‑148 (MQ=255)                     gcgTAGAGCCCGGCGTCCCCGTCGATGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAg                                                                                                                         <  1:140604/146‑1 (MQ=255)                     gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAg                                                                                                                         <  3:180460/146‑1 (MQ=255)                     gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGgc                                                                                                                       >  2:172325/1‑148 (MQ=255)                           agCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                  <  1:20854/147‑1 (MQ=255)                                ggCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACg                                                                                                            >  8:18144/1‑148 (MQ=255)                                  cGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGat                                                                                                          <  3:97319/148‑1 (MQ=255)                                     ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                        >  6:93551/1‑147 (MQ=255)                                       ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCga                                                                                                     <  6:14746/148‑1 (MQ=255)                                                cgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                            <  2:145333/148‑1 (MQ=255)                                                    cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                        <  4:141740/148‑1 (MQ=255)                                                              gTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAtt                                                                              <  6:127960/148‑1 (MQ=255)                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                           <  1:50328/145‑1 (MQ=255)                                                                         gATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                   >  5:56562/1‑148 (MQ=255)                                                                           tGGTCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGAGCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcgc                                                                  <  7:84447/147‑1 (MQ=255)                                                                               cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCa                                                             >  6:185134/1‑148 (MQ=255)                                                                                tACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAACC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                            >  1:98263/1‑148 (MQ=255)                                                                                     cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                       <  6:56562/148‑1 (MQ=255)                                                                                       tgctccTGGGGGGGGAGGGTGATATGGTGGACGACCTTATATGCCTGGGGGCTGCC‑GATTACGTGTGCGACATGGCCAATGCGCTGATGATCGATATCGAGATCGCCAAGGGTAAGGCGCATTAGCCATTGAAGGAGTCATTTGTggg                                                     <  6:157749/144‑1 (MQ=255)                                                                                       aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                      <  5:93551/147‑1 (MQ=255)                                                                                             tGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGc                                               >  3:70820/1‑148 (MQ=255)                                                                                             tGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAACC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGc                                               <  2:98263/148‑1 (MQ=255)                                                                                                     gAGGGCGATATGGAGGCCGGCCATGGATTCCTGGGGGCAGCCGGT‑TCCCTGTGCGACATGGCCCCGGCGCTGCTGCCGGAGGTCGAGATCGCCACGGGGCCGGGGCGGTTGGCATTGCCGCGGTCAGTTGTGGGAGCGGCCTTGCGTc                                       >  6:122475/1‑148 (MQ=255)                                                                                                            atatGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATg                                <  7:138496/148‑1 (MQ=255)                                                                                                            atatGGAGGACGACAATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGATTATGAACGATATCGAGATCGCCAATGGTAAGGCGGATTAGGCATTGACTCGGTCATTTGTTGGAGCGGCCTTGCGTCGCGAATg                                <  4:127370/148‑1 (MQ=255)                                                                                                              atGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                              <  7:18144/148‑1 (MQ=255)                                                                                                                ggaggaCGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGc                             <  4:70820/147‑1 (MQ=255)                                                                                                                ggaggaCGACCATAAATTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                            <  1:91478/148‑1 (MQ=255)                                                                                                                             aaaTTCCTGGGGGCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCTTTGCAGccca               <  3:91310/148‑2 (MQ=255)                                                                                                                                        ggCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                            >  2:28875/1‑60 (MQ=255)                                                                                                                                         gCAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAg                                                                                           <  1:28875/60‑1 (MQ=255)                                                                                                                                          cAGCC‑GATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCagag  <  5:74762/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             CCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCC‑GA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAG  >  NC_002947/335981‑336264

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