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	98.0%	166.4 / ‑3.1	50	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (24/25);  total (25/25)
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.36e-01

ATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGAT  >  NC_002947/335984‑336266                                                                                                                                              |                                                                                                                                                aTCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCtgt                                                                                                                                         <  1:167977/149‑1 (MQ=255) aTCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCtgt                                                                                                                                         >  6:23115/1‑149 (MQ=255) aTCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCtgt                                                                                                                                         >  3:172467/1‑149 (MQ=255)    ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                       <  4:172467/148‑1 (MQ=255)     cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                       >  8:183977/1‑147 (MQ=255)     cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                       <  8:136315/147‑1 (MQ=255)     cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                       <  7:183977/147‑1 (MQ=255)     cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                       >  7:136315/1‑147 (MQ=255)     cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCg                                                                                                                                      >  6:5524/1‑148 (MQ=255)         aTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATg                                                                                                                                 <  5:5297/149‑1 (MQ=255)              cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGc                                                                                                                               <  5:5524/146‑1 (MQ=255)               tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAg                                                                                                                           >  1:174888/1‑149 (MQ=255)                     tAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCt                                                                                                                      >  5:343078/1‑148 (MQ=255)                       gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                    <  5:23115/148‑1 (MQ=255)                        agCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                   >  4:133952/1‑148 (MQ=255)                         gCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtgatga                                                                                                                 <  3:133952/149‑1 (MQ=255)                          cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCtgt                                                                                                                                         <  5:187617/124‑1 (MQ=255)                          cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCtgt                                                                                                                                         >  6:187617/1‑124 (MQ=255)                                 tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGTCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                         <  3:246121/149‑1 (MQ=255)                                  ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                         >  3:63089/1‑148 (MQ=255)                                            gcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGcc                                                                                                >  8:78430/1‑147 (MQ=255)                                                  tACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                          >  1:219639/1‑147 (MQ=255)                                                     gAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgc                                                                                     >  1:120254/1‑149 (MQ=255)                                                       gCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                    >  4:115884/1‑148 (MQ=255)                                                         aTGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                                 >  3:316009/1‑149 (MQ=255)                                                         aTGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                                 <  7:78430/149‑1 (MQ=255)                                                           gTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGata                                                                                                           <  2:340517/121‑1 (MQ=255)                                                           gTCTCGATCCTGATGGGCTACATAAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGata                                                                                                           >  1:340517/1‑121 (MQ=255)                                                              tcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                             >  3:288353/1‑148 (MQ=255)                                                                 aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                         <  1:128271/149‑1 (MQ=255)                                                                     tGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                     >  8:358811/1‑149 (MQ=255)                                                                       aTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcgc                                                                    <  4:288353/148‑1 (MQ=255)                                                                        tGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGAc                                                                                                                                    >  5:23552/1‑83 (MQ=255)                                                                        tGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGAc                                                                                                                                    <  6:23552/83‑1 (MQ=255)                                                                            cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                              >  5:7081/1‑149 (MQ=255)                                                                                aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgt                                                          <  6:343078/149‑1 (MQ=255)                                                                                  cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                        >  1:255454/1‑149 (MQ=255)                                                                                   aaGCACCTTGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                        <  4:63089/148‑1 (MQ=255)                                                                                   aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                        >  6:289816/1‑148 (MQ=255)                                                                                        ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAgg                                                                                                                          >  7:165858/1‑77 (MQ=255)                                                                                        ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAgg                                                                                                                          <  8:165858/77‑1 (MQ=255)                                                                                                           ttGGAGTACGACGATAAATTCCTGGGTGCAGCCGATT‑ACGTGTGCGAGATGGGCAAGGCGGTGATGTACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAGGCGTTGATTTGTGGGAGCGGCCTTGCGTCGCGAATgggg                               <  2:89912/148‑2 (MQ=255)                                                                                                            tGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                              <  2:120254/149‑1 (MQ=255)                                                                                                             ggaggaCGACCATAAATTCCTGGGGGCAGCCGA‑TAACCTGTGCGACATGGCCAAGGCCCTCATCAACGCTATCGAGATCGCCAAGGTTAAGGCGCATTAGGCATTGACGCGGTCATATGAGGGAGAGGCCTTGCGTCGCGAATGGGc                               >  1:89912/1‑147 (MQ=255)                                                                                                              gaggaCGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                             <  6:7081/148‑1 (MQ=255)                                                                                                              gaggaCGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                             <  2:219639/148‑1 (MQ=255)                                                                                                                           aaTTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCgg               <  2:174888/149‑1 (MQ=255)                                                                                                                             ttCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCg             <  3:115884/149‑1 (MQ=255)                                                                                                                                   ggggCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt         <  7:358811/147‑1 (MQ=255)                                                                                                                                        aGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGAt  >  3:329996/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                ATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGAT  >  NC_002947/335984‑336266

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