Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,712,080	A→T	21.2%	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	21.2%	45.5 / 4.2	27	intergenic (+71/‑187)	ppc/adk	phosphoenolpyruvate carboxylase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base A (16/4);  new base T (2/4);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 5.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

GCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCG  >  NC_002947/1711932‑1712227                                                                                                                                                       |                                                                                                                                                    gCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCa                                                                                                                                                     <  2:40304/149‑1 (MQ=255)               gTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCCACCGCCTg                                                                                                                                       >  8:175745/1‑149 (MQ=255)                ttGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGc                                                                                                                                      >  4:77093/1‑149 (MQ=255)                ttGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCAGc                                                                                                                                      >  2:367439/1‑149 (MQ=255)                                 ggTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTg                                                                                                                     >  5:314647/1‑149 (MQ=255)                                  gTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCCGGCGCAACCCCCTGCCCGGTCAGCCGCAAGTgg                                                                                                                    >  8:133208/1‑149 (MQ=255)                                  gTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTgg                                                                                                                    >  4:322184/1‑149 (MQ=255)                                       cGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCAt                                                                                                               >  1:271547/1‑149 (MQ=255)                                                tATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAAc                                                                                                      >  3:171968/1‑149 (MQ=255)                                                                   aaCACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGCA‑‑GGCGGTTGCGCCTGGCGCAAcc                                                                                                                                            >  8:260902/1‑91 (MQ=38)                                                                   aaCACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGCA‑‑GGCGGTTGCGCCTGGCGCAAcc                                                                                                                                            <  7:260902/91‑1 (MQ=38)                                                                        cGGCTAAGGGCAAGCCCGGGTCTTGTCTGGTGCCGGCCTTACCGGCCCGGCACATGGCG‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCg                                                                              <  7:175745/149‑1 (MQ=255)                                                                                ggCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCGACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTc                                                                                       <  1:280078/132‑1 (MQ=255)                                                                                ggCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCGACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTc                                                                                       >  2:280078/1‑132 (MQ=255)                                                                                       ccGGGTCTTGTGTGGTTCGGGCCATACCGGCCCGTCACATGGCC‑GGGGT‑GGGGGGTTGCGCCTGGGGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCtt                                                               <  3:77093/149‑1 (MQ=255)                                                                                       ccGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCtt                                                               >  5:115463/1‑149 (MQ=255)                                                                                              aTGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGagcagc                                                        >  5:360899/1‑149 (MQ=255)                                                                                                  gggggtccGGCCAAACCGGCCCGGCAAAAGGGCGGGGGC‑GGGGTGTGGCGCCTGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCtt                                                     <  3:322184/142‑1 (MQ=255)                                                                                                              aTACCGGCCCGGCAAATGGAG‑GGGGT‑GGGGTGTTGCGCCTGGCGAAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTAt                                        <  6:360899/149‑1 (MQ=255)                                                                                                               tACAGGCCCGGCACAAGGAC‑GGGGAGGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCt                                                                                                           >  1:193003/1‑82 (MQ=255)                                                                                                               tACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATg                                       >  3:293682/1‑149 (MQ=255)                                                                                                                       cGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATAcc                               >  7:259396/1‑149 (MQ=255)                                                                                                                                       gA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACtt                <  8:259396/149‑1 (MQ=255)                                                                                                                                          aggCGGTTGCGCCTGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTgg              >  6:196597/2‑149 (MQ=255)                                                                                                                                            gccggtGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCACGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGcc            <  4:293682/144‑1 (MQ=255)                                                                                                                                                    gcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc    >  6:62448/1‑149 (MQ=255)                                                                                                                                                    gcCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGc    >  2:150310/1‑149 (MQ=255)                                                                                                                                                      cAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGcc   >  6:106545/1‑148 (MQ=255)                                                                                                                                                       aGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCg  >  5:333036/1‑148 (MQ=255)                                                                                                                                                       aGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCg  >  7:301079/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGAC‑GGGGA‑GGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCG  >  NC_002947/1711932‑1712227

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