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.3%	119.9 / ‑3.2	37	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (18/18);  total (19/18)
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.76e-01

CCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCG‑AATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335981‑336261                                                                                                                                                 |                                                                                                                                            cccATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑ttac                                                                                                                                         <  4:97897/148‑4 (MQ=255)     tCCCTCAATCTCCTACGCGTAGAGACCGGCGTCACCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGTGCTACATCAAGCACCTGTTGCGCGAGGGCTATATGGAGGACGACCATAACTTCCTGGGTGCAGCCGATT‑ACCtgt                                                                                                                                     <  4:39005/148‑1 (MQ=255)        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                   <  3:41611/147‑1 (MQ=255)        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCg                                                                                                                                  >  1:168648/1‑148 (MQ=255)        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑TGACCTGTGCg                                                                                                                                  >  3:129650/1‑148 (MQ=255)             tctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATg                                                                                                                             <  4:55053/148‑1 (MQ=255)                               cGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAAc                                                                                                           <  7:3464/148‑1 (MQ=255)                                     ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                     <  1:120236/148‑1 (MQ=255)                                           cGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                               <  5:134077/148‑1 (MQ=255)                                                cgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                          <  3:67282/148‑1 (MQ=255)                                                     tACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAggg                                                                                       >  4:24738/1‑146 (MQ=255)                                                        gAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAgg                                                                                  >  5:75170/1‑148 (MQ=255)                                                            aTGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCa                                                                              >  6:105937/1‑148 (MQ=255)                                                              gTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCa                                                                              >  7:38501/1‑146 (MQ=255)                                                               tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTa                                                                           >  3:44957/1‑148 (MQ=255)                                                                   gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCaa                                                                                                                        >  1:67229/1‑99 (MQ=255)                                                                   gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCa                                                                       >  1:166410/1‑148 (MQ=255)                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAg                                                                                                                       <  2:67229/99‑1 (MQ=255)                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                         <  3:38164/145‑1 (MQ=255)                                                                         gATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                      >  1:145348/1‑111 (MQ=255)                                                                          aTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                     <  2:145348/111‑1 (MQ=255)                                                                               cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCa                                                           >  6:14092/1‑148 (MQ=255)                                                                                tACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                          >  1:166907/1‑148 (MQ=255)                                                                                   aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                          >  6:105911/1‑129 (MQ=255)                                                                                   aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                       >  7:10702/1‑148 (MQ=255)                                                                                   aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                       >  5:56048/1‑148 (MQ=255)                                                                                    tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                         <  5:105911/129‑1 (MQ=255)                                                                                     cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                     <  7:119929/148‑1 (MQ=255)                                                                                          aCCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGc                                                >  3:157638/1‑148 (MQ=255)                                                                                                tGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCACGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCt                                           <  1:5836/147‑1 (MQ=255)                                                                                                   gcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCg                                       <  8:10702/148‑1 (MQ=255)                                                                                                       gggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgc                                   <  6:65080/148‑1 (MQ=255)                                                                                                                  aggaCGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCG‑AATGGGCTGc                       <  5:120284/148‑1 (MQ=255)                                                                                                                         ccATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCGTGCGTCGCG‑CATGGGGTGCAAagacg                >  7:70263/1‑145 (MQ=255)                                                                                                                          cATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCG‑AATGGGCTGCAAagcagc               <  3:61416/148‑1 (MQ=255)                                                                                                                          cATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCATTCGTCGCGCAATGGG‑TGCAAagcacc               >  1:17721/1‑146 (MQ=255)                                                                                                                            tAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCATTCGTCGCG‑CAAGGGGGGCAAACCAGccc             >  4:67603/1‑148 (MQ=255)                                                                                                                                       gggCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATAGCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGGCATTTGGGGGAGCGGCCTTGCGTCGCT‑AATGGGCTGCAAAGCAGCCCCGGCGATTTCa  <  5:14092/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            CCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCG‑AATGGGCTGCAAAGCAGCCCCGGCGATTTCA  >  NC_002947/335981‑336261

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