Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,712,080	A→T	17.6%	intergenic (+71/‑187)	ppc → / → adk	phosphoenolpyruvate carboxylase/adenylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,712,080	0	A	T	17.6%	75.9 / 4.4	42	intergenic (+71/‑187)	ppc/adk	phosphoenolpyruvate carboxylase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base A (21/11);  new base T (2/6);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 5.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGC  >  NC_002947/1711938‑1712225                                                                                                                                                |                                                                                                                                                  gACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGgcgc                                                                                                                                              >  7:14663/1‑148 (MQ=255) gACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGgc                                                                                                                                                <  2:74547/146‑1 (MQ=255) gACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGgc                                                                                                                                                <  4:35744/146‑1 (MQ=255)   cAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAAc                                                                                                                                           <  6:191885/149‑1 (MQ=255)    aGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCTGGCGCAAcc                                                                                                                                          <  4:381411/149‑1 (MQ=255)    aGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAAcc                                                                                                                                          <  7:34313/149‑1 (MQ=255)    aGCCCGGTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAAcc                                                                                                                                          <  6:52078/149‑1 (MQ=255)      cccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGc                                                                                                                                        >  6:2041/1‑149 (MQ=255)      cccGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGc                                                                                                                                        >  6:363486/1‑149 (MQ=255)             gAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTTCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAATGACGGGGAGGGCCGTTTCGCCAGGCGCCACCGCCTGCCCg                                                                                                                                 >  1:209514/1‑149 (MQ=255)               gCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGt                                                                                                                               >  5:12973/1‑149 (MQ=255)                cAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGt                                                                                                                               >  7:109145/1‑148 (MQ=255)                            gTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCg                                                                                                                         >  1:229837/1‑142 (MQ=255)                             tgcgcGTGGGGGTTTTGGCCGCGGGGGTGCCCACCACGGTCTATGGCCACGCCGGGGTCTTGTGTGGT‑CCGGCCAATACCGGCCGGTCAAATGGCGGGGGCGGGGGGTGGCCCCCGGCCAAACCGCCTGCCCGGTCAGCCGCAAGTGgc                                                                                                                 <  8:109145/145‑1 (MQ=255)                                 cGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCAt                                                                                                             >  3:251319/1‑149 (MQ=255)                                       gggTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTg                                                                                                        >  6:26598/1‑148 (MQ=255)                                                 gccgGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTTCGCCCGGCGCAACCGCCAGCCCGGTCAGCCGCCAGTGGCGCATTGCTGCAACTTTtggg                                                                                             >  6:217725/1‑148 (MQ=255)                                                        tGCGCAACACCGGCTAAGGGCAAGCCCGGGACATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGt                                                                                      >  1:31635/1‑149 (MQ=255)                                                           ccAACACCGGATAAGGGCAAGCCGGGTTCATGTCTGGTGCCGGCA‑ATACCGGCCCGGAACAAGTCGGGGGTGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTTGGACGCTTGTCTg                                                                                   <  5:363486/148‑1 (MQ=255)                                                                    gCTACGGCCAAGCCGGGGTCAGGTCGGGGGACGGCC‑ATACCGGCCCGGAAAAAGGAGGGGGAGGGCGGTGGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCggg                                                                          <  5:177857/149‑1 (MQ=255)                                                                     cTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTACTGCAACTTTGGGACGCTTGTCTGGCTGCCGGgc                                                                         >  5:57483/1‑149 (MQ=255)                                                                               gCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCTAGGCGCAACCGCCTGCCCGGTCAGCCGCGAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATc                                                               >  3:243039/1‑149 (MQ=255)                                                                                cccGGGTCTTGTCTGGGGCCGGCC‑ATACCGGCCCGCCACAAGGAGGGGTCGGGCGGTGGCGCCAGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCt                                                              <  3:38584/149‑1 (MQ=255)                                                                                  cGGGTCATGTCGGGTGACGGCC‑ATACCGGCCCGGCACAAGGCCGGGGAGGGCGGGTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTg                                                            <  7:355449/149‑1 (MQ=255)                                                                                   gggTCTTGTGGGGTGCCGCCC‑ATCCGGGCCGGCCAAACGGCGGGGGGGGGGGGTGGCCCCTGGCCCACCCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGa                                                           <  5:2041/149‑1 (MQ=255)                                                                                           tCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                   >  1:229873/1‑149 (MQ=255)                                                                                                    ggCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGAt                                          >  4:80136/1‑149 (MQ=255)                                                                                                       c‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTAGGGACGCTTGTATGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTa                                       >  1:169275/1‑149 (MQ=255)                                                                                                             aGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAAGTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGAATTATGGt                                   >  7:164074/1‑148 (MQ=255)                                                                                                                     cacaAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGa                          <  1:30358/149‑1 (MQ=255)                                                                                                                      acaAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGAt                         >  7:189667/1‑149 (MQ=255)                                                                                                                          ggACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTAGCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTc                     <  2:31635/149‑1 (MQ=255)                                                                                                                           gcgggggggggggGTGGCCCCTGGCCAAACCCCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTc                     <  5:26598/145‑1 (MQ=255)                                                                                                                              ggggAGGGGGGTTGGGCTTGCCGAAACGGCCTGCCCGGTCAGCCGCGAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCgg                  <  4:243039/148‑1 (MQ=255)                                                                                                                                ggAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTAAGCCGCAAGTGGCGCATTGCTGCCACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGCTTTTCCGGACt               >  2:360244/1‑149 (MQ=255)                                                                                                                                 gcgggCGGTGGGCCCTGGGGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCGTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACtt              <  8:164074/147‑1 (MQ=255)                                                                                                                                   aggCGGTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTgg            >  5:364358/2‑149 (MQ=255)                                                                                                                                     gCGGTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGcc          >  6:282974/1‑149 (MQ=255)                                                                                                                                         ttGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTgg      >  4:45465/1‑149 (MQ=255)                                                                                                                                           gcgcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTc                                  <  3:98143/119‑1 (MQ=255)                                                                                                                                           gcgcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTc                                  >  4:98143/1‑119 (MQ=255)                                                                                                                                             gcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc  >  6:178949/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                  GACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCC‑ATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGC  >  NC_002947/1711938‑1712225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑