Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,339,540	G→A	15.4%	G61D (GGC→GAC)	dctA‑II →	C4‑dicarboxylate transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,339,540	0	G	A	15.4%	57.1 / 4.1	26	G61D (GGC→GAC)	dctA‑II	C4‑dicarboxylate transport protein
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base A (2/2);  total (13/13)
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.46e-01

GCCGTCCTGCTGGGCGTGGTGGTCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATCGATCCTGCCACCTT  >  NC_002947/2339398‑2339687                                                                                                                                               |                                                                                                                                                    gcCGTCCTGCTGGGCGTGGTGGTCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGAc                                                                                                                                               <  2:58145/149‑1 (MQ=255)  ccGTCCTGCTGGGCGTGGTGGTCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGAcc                                                                                                                                              <  8:134867/149‑1 (MQ=255)       ctgctgGGCGTGGTGGTCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACa                                                                                                                                          <  5:55100/148‑1 (MQ=255)                    tggtCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTgggg                                                                                                                            <  4:147618/149‑1 (MQ=255)                            cacTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTtgat                                                                                                                    <  6:49302/149‑1 (MQ=255)                                      cgcAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGAcc                                                                                                          <  4:149002/149‑1 (MQ=255)                                                                   ggTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACggtggt                                                                             >  4:19126/1‑149 (MQ=255)                                                                       aTGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTtg                                                                                                                      <  5:6481/104‑1 (MQ=255)                                                                       aTGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTtg                                                                                                                      >  6:6481/1‑104 (MQ=255)                                                                            ttCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGt                                                                    >  6:275754/1‑149 (MQ=255)                                                                                     cTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACt                                                           >  1:124813/1‑149 (MQ=255)                                                                                          cAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTggc                                                        <  2:300610/147‑1 (MQ=255)                                                                                          cAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTggc                                                        >  1:300610/1‑147 (MQ=255)                                                                                           aaGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTggcggc                                                     >  3:270145/1‑149 (MQ=255)                                                                                            aGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGCCCCTGCTGCTGTTCCTGGGGCTGAAGGTGGTCAGCCTGGTGATCTGCCTGgcggcg                                                    >  3:197326/1‑149 (MQ=255)                                                                                            aGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGgcggcg                                                    <  5:115302/149‑1 (MQ=255)                                                                                                    atCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTggcg                                                       >  7:48220/1‑138 (MQ=255)                                                                                                    atCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTggcg                                                       <  8:48220/138‑1 (MQ=255)                                                                                                        cgcCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTc                                        <  2:124813/149‑1 (MQ=255)                                                                                                                   ttCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATGGGCCTGGCGGCGGTCTTCATGTGCCAGCCAGGTGc                             >  8:273743/1‑149 (MQ=255)                                                                                                                      tGCACCATCGTCACCGGCATTGCCGACATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCg                           >  7:156461/1‑148 (MQ=255)                                                                                                                      tGCACCATCGTCACCGGCATTGCCGACATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCg                           >  1:168074/1‑148 (MQ=255)                                                                                                                        cACCATCGTCACCGGCATTGCCGACATTACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCg                           <  8:156461/146‑1 (MQ=255)                                                                                                                        cACCATCGTCACCGGCATTGCCGACATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCg                           <  2:168074/146‑1 (MQ=255)                                                                                                                         aCCATCGTCACCGGCATTACCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCAt                       >  2:45381/1‑149 (MQ=255)                                                                                                                                tCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATc                >  2:32575/1‑149 (MQ=255)                                                                                                                                              ggCCTGACCGACAAGCGTGCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATCGATCCTGCCACCtt  <  7:135722/149‑1 (MQ=255)                                                                                                                                              ggCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATCGATCCTGCCACCtt  <  5:275754/149‑1 (MQ=255)                                                                                                                                               gCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATCGATCCTGCCACCtt  >  7:160786/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                    GCCGTCCTGCTGGGCGTGGTGGTCGGACACTTCTGGGCGCAACAGGCCGTGGCCCTCAAGCCGCTGGGTGATGCGTTCATCAAGCTGATCAAGATGATGATCGCGCCGGTGGTGTTCTGCACCATCGTCACCGGCATTGCCGGCATGACCGACAAGCGTTCACTGGGGCGCTTGATGAGCAAGACCCTGCTGCTGTTCCTGGGTCTGACGGTGGTCAGCCTGGTGATCGGACTGGCGGCGGTCTACCTGTTCAAGCCAGGTGCCGGCATGAACATCGATCCTGCCACCTT  >  NC_002947/2339398‑2339687

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