Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,949,878	A→G	9.3%	M199V (ATG→GTG)	ppx →	exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,949,878	0	A	G	9.3%	118.1 / 3.3	43	M199V (ATG→GTG)	ppx	exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (2/2);  total (22/21)
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.27e-01

AATTCATCATTGGCCAGCGCTTCGAGCCGCTGCTGCGCGAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCAAACGC  >  NC_002947/5949732‑5950024                                                                                                                                                   |                                                                                                                                                    aaTTCATCATTGGCCAGCGCTTCGAGCCGCTGCTGCGCGAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACtgatg                                                                                                                                                   <  4:143147/149‑1 (MQ=255)         atTGGCCAGCGCTTCGAGCCGCTGCTGCGCGAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGa                                                                                                                                           <  2:24810/149‑1 (MQ=255)                           ccGCTGCTGCGCGAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCt                                                                                                                        >  3:52448/1‑149 (MQ=255)                           cagctgctGCGCGAAAGCCTGCCGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCt                                                                                                                        <  8:230619/147‑1 (MQ=255)                                     gcgAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGAt                                                                                                              >  3:34332/1‑149 (MQ=255)                                       gAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACtgatg                                                                                                                                                   <  5:31661/111‑1 (MQ=255)                                       gAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACtgatg                                                                                                                                                   >  6:31661/1‑111 (MQ=255)                                                 aGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGc                                                                                                  <  4:52448/149‑1 (MQ=255)                                                   aTGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAg                                                                                                           >  4:27833/1‑138 (MQ=255)                                                   aTGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAg                                                                                                           <  3:27833/138‑1 (MQ=255)                                                   aTGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCt                                                                                                 >  7:87417/1‑148 (MQ=255)                                                       gTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATcc                                                                                            <  6:174930/149‑1 (MQ=255)                                                        ttGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGGTGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCg                                                                                           >  4:261473/1‑149 (MQ=255)                                                        ttGCGTGAGCTATACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGGAGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTAACGGCTCATCCg                                                                                           >  2:6536/1‑149 (MQ=255)                                                            gTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCAc                                                                                       >  2:103734/1‑149 (MQ=255)                                                             tGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCgg                                                                                                   <  5:97902/136‑1 (MQ=255)                                                             tGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCgg                                                                                                   >  6:97902/1‑136 (MQ=255)                                                                 cTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTc                                                                                  <  7:160253/149‑1 (MQ=255)                                                                             cTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGgcg                                                                      >  4:71184/1‑149 (MQ=255)                                                                                     gcgACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCAt                                                               >  5:203838/1‑148 (MQ=255)                                                                                      cgACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCACCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATc                                                              >  4:134452/1‑148 (MQ=255)                                                                                                     cccGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑cc                                                                                                                                   <  1:32710/64‑1 (MQ=255)                                                                                                     cccGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑cc                                                                                                                                   >  2:32710/1‑64 (MQ=255)                                                                                                     cccGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGTCATCAAGGCCGGCGGCCTgg                                              <  3:71184/149‑1 (MQ=255)                                                                                                      ccGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCACCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTggg                                             <  3:134452/149‑1 (MQ=255)                                                                                                           ccGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATg                                        <  8:87417/149‑1 (MQ=255)                                                                                                            cGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATgg                                       >  1:75383/1‑149 (MQ=255)                                                                                                                 tGCCCAGGCCTACACTGCCGCACGCCTGGAACTGGTGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGa                                    <  3:261473/147‑1 (MQ=255)                                                                                                                 tGCCCAGGCCTACACTGCCGCACGCCTGGAACTGGTGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATACAGGCCGGCGGCCTGGGCAATGGCGa                                    <  1:6536/147‑1 (MQ=255)                                                                                                                    ccAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTc                                                                                  <  5:275274/98‑1 (MQ=255)                                                                                                                    ccAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTc                                                                                  >  6:275274/1‑98 (MQ=255)                                                                                                                        gCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGc                           >  8:26805/1‑149 (MQ=255)                                                                                                                              acTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGa                                    >  1:198347/1‑134 (MQ=255)                                                                                                                              acTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGa                                    <  2:198347/134‑1 (MQ=255)                                                                                                                                tGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAggcc                   >  4:78997/1‑149 (MQ=255)                                                                                                                                  ccGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAggcct                  <  2:75383/148‑1 (MQ=255)                                                                                                                                    gCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAggcctggc               >  1:20455/1‑149 (MQ=255)                                                                                                                                      aCGCCTGGAACTGATGAGCATCGAAAATGCCCCTGCACCGCCTGACCTGGGGTGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAggcctggcc              >  1:28494/1‑149 (MQ=255)                                                                                                                                       cGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTc                                <  3:70988/129‑1 (MQ=255)                                                                                                                                       cGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTc                                >  4:70988/1‑129 (MQ=255)                                                                                                                                           tGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGcctggcctggc          >  5:32038/1‑146 (MQ=255)                                                                                                                                           tGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTggg          <  6:32038/147‑1 (MQ=255)                                                                                                                                            ggAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTg                                               <  6:16354/109‑1 (MQ=255)                                                                                                                                            ggAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTg                                               >  5:16354/1‑109 (MQ=255)                                                                                                                                              aaCTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCaaa     <  1:250088/149‑1 (MQ=255)                                                                                                                                                cTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCAAACg   <  4:34332/149‑1 (MQ=255)                                                                                                                                                   atgaGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCAAACGc  >  2:90813/1‑147 (MQ=255)                                                                                                                                                   atgaGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCAAACGc  >  2:93934/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AATTCATCATTGGCCAGCGCTTCGAGCCGCTGCTGCGCGAAAGCCTGCAGATGGGTTGCGTGAGCTTTACCCAGCGCTATTTCCGCGACGGCAAGATCACCCCGGCCCGCTATGCCCAGGCCTACACTGCCGCACGCCTGGAACTGATGAGCATCGAAAATG‑CCCTGCACCGCCTGACCTGGGATGAGGCTATCGGCTCATCCGGCACCATTCGCGCCATCGGCGCGGCCATCAAGGCCGGCGGCCTGGGCAATGGCGAGGTCAACGCCGAAGGCCTGGCCTGGGTCAAACGC  >  NC_002947/5949732‑5950024

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