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	94.7%	118.5 / 1.0	38	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (14/22);  total (15/23)
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.54e-01

ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335980‑336261                                                                                                                                                  |                                                                                                                                            aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGGCGACCATAAATTCCTGCGGGCAGCCGA‑T‑a                                                                                                                                          >  1:51098/1‑147 (MQ=255)  cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tta                                                                                                                                          <  1:97213/147‑3 (MQ=255)  cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tta                                                                                                                                          >  7:165864/1‑145 (MQ=255)   ccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑ttac                                                                                                                                         >  2:172455/1‑144 (MQ=255)       cccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCtgt                                                                                                                                    <  1:381085/147‑1 (MQ=255)          tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCg                                                                                                                                 >  4:213360/1‑147 (MQ=255)          tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATACTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACAATCAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGc                                                                                                                                  >  8:134799/1‑146 (MQ=255)               ctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATg                                                                                                                            >  3:177157/1‑147 (MQ=255)               ctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATg                                                                                                                            >  3:308201/1‑147 (MQ=255)                  cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTTCATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGc                                                                                                                          >  6:285735/1‑146 (MQ=255)                   tACGCGTAGAGCCCGGCGTCCGCGTCGACGCGGCCTACGAGCATGTATCGATCCTGATTGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATgg                                                                                                                           <  8:228514/144‑1 (MQ=255)                    aCGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCACGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCaa                                                                                                                       >  3:330196/1‑147 (MQ=255)                     cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAg                                                                                                                      <  3:213360/147‑1 (MQ=255)                      gcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAg                                                                                                                      <  2:94521/146‑1 (MQ=255)                          agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCt                                                                                                                 <  1:172455/147‑1 (MQ=255)                            agCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCtg                                                                                                                <  8:165864/146‑1 (MQ=255)                                      ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                     <  3:18868/147‑1 (MQ=255)                                        ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                     >  5:275256/1‑145 (MQ=255)                                        ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                    <  6:4750/146‑1 (MQ=255)                                                        cGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTaa                                                                                   >  4:28720/1‑147 (MQ=255)                                                            cATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                               <  4:177157/147‑1 (MQ=255)                                                            cATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgc                                                                                >  7:68046/1‑146 (MQ=255)                                                                tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATCAACGATATCGAGATCGCCAAGGGTAAGGCGCa                                                                             >  4:174293/1‑145 (MQ=255)                                                                     aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                       <  4:308201/146‑1 (MQ=255)                                                                          gATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAc                                                                 >  6:318750/1‑147 (MQ=255)                                                                           aTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                     >  3:311667/1‑110 (MQ=255)                                                                           aTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                     <  4:311667/110‑1 (MQ=255)                                                                                    aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                       <  8:68046/147‑1 (MQ=255)                                                                                      cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgt                                                     <  3:110629/147‑1 (MQ=255)                                                                                            ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGc                                               <  3:28720/147‑1 (MQ=255)                                                                                                gTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑AGCGTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCCCCAAGGGTAAGGCGCATTAGGCATTCACGCGGGCATTTGTGGGAGCGGc                                            >  1:294083/1‑147 (MQ=255)                                                                                                    gcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTg                                        <  6:275256/146‑1 (MQ=255)                                                                                                    gcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTg                                        <  1:78691/146‑1 (MQ=255)                                                                                                         ggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATCTCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgc                                  <  8:353322/147‑1 (MQ=255)                                                                                                                         acCATAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagc                  <  1:51181/147‑1 (MQ=255)                                                                                                                             tAAATTCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcc              <  5:318750/147‑1 (MQ=255)                                                                                                                                  tCCTGGGGGCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCgg          >  4:164552/1‑146 (MQ=255)                                                                                                                                        gggCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt    <  1:98056/146‑1 (MQ=255)                                                                                                                                        gggCAGCCGATT‑A‑CCTGTGCCACATGGCCCAGGCGCTGATTAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGAATTGCCTCGGTCTTTTTTGGGAGCGGCCTTGCGTCGCGAATTGGCTGCAAAGCAGCCCCGGCGAtt     <  7:239289/145‑1 (MQ=255)                                                                                                                                        gggCAGCCGA‑TAA‑CCTGTGCGACATGGCCAAGGCGCTCACGAACCATATCGAGATCGCCCAGGGTAAGGCGCATTAGGCATTGACGCGGCCATTTGGGGGTGCGGCCTTGCGTCGCGAATGGGCTGCAACGCAGCCCCGGCGAtt     <  4:330196/145‑1 (MQ=255)                                                                                                                                         ggCAGCCGATT‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCa  <  2:294083/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑A‑CCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335980‑336261

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