Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,688,198	A→G	14.3%	A105A (GCA→GCG)	PP_3253 →	carboxylate‑amine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,688,198	0	A	G	14.3%	63.5 / 4.0	29	A105A (GCA→GCG)	PP_3253	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (2/2);  total (13/16)
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.53e-01

GCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTT  >  NC_002947/3688051‑3688341                                                                                                                                                     |                                                                                                                                                gCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGca                                                                                                                                                 <  2:60428/148‑1 (MQ=255)  ccGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGcac                                                                                                                                                >  2:103676/1‑148 (MQ=255)     gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGTGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAgg                                                                                                                                            <  1:166485/149‑1 (MQ=255)     gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCCCAAgg                                                                                                                                            >  3:20918/1‑149 (MQ=255)     gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAgg                                                                                                                                            >  6:233231/1‑149 (MQ=255)     gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAgg                                                                                                                                            <  5:146808/149‑1 (MQ=255)                    tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                               <  3:182162/147‑1 (MQ=255)                       gAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCTTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGCGCAggggtggctt                                                                                                                                     <  5:96433/138‑7 (MQ=25)                       gAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCTTGTGCCGCG‑GGCTTCTGCAGCAGCCACGCGGCGCAggggtggctt                                                                                                                                     >  6:96433/1‑132 (MQ=25)                        aaGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTc                                                                                                                         >  4:252940/1‑149 (MQ=255)                                               aaCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                                                  >  4:125151/1‑149 (MQ=255)                                                    tCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCg                                                                                               >  7:213675/1‑147 (MQ=255)                                                                cgTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATg                                                                                 <  1:179133/149‑1 (MQ=255)                                                                      aTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGc                                                                           <  1:5706/149‑1 (MQ=255)                                                                                   gggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTgg                                                              <  3:125151/149‑1 (MQ=255)                                                                                    ggATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGc                                                             <  3:252940/149‑1 (MQ=255)                                                                                      aTGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGt                                                           <  1:103676/149‑1 (MQ=255)                                                                                             cATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgccgc                                                      >  7:82648/1‑147 (MQ=255)                                                                                                      cGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGAc                                           <  8:213675/149‑1 (MQ=255)                                                                                                              cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCAGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGc                                    >  5:27711/1‑148 (MQ=255)                                                                                                                ccccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGCCCGCATGCAAc                                 >  5:254621/1‑149 (MQ=255)                                                                                                                   ttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGCAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                                  >  1:153418/5‑145 (MQ=38)                                                                                                                   ttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGCAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                                  <  2:153418/141‑1 (MQ=38)                                                                                                                        ccgcGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGTCCGCATGCAACTGATCAAc                         <  5:233231/149‑1 (MQ=255)                                                                                                                            gTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAtgtg                                                                      >  8:181417/1‑100 (MQ=255)                                                                                                                                      cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTg           <  6:254621/149‑1 (MQ=255)                                                                                                                                      cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTg           >  7:137222/1‑149 (MQ=255)                                                                                                                                      cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTg           >  1:43007/1‑149 (MQ=255)                                                                                                                                        gAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGt         <  2:43007/149‑1 (MQ=255)                                                                                                                                            cc‑cGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTAcc     <  8:82648/149‑1 (MQ=255)                                                                                                                                               gatctcacaAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCgt   <  5:7370/144‑1 (MQ=255)                                                                                                                                                gcgGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCgtt  <  4:20918/149‑1 (MQ=255)                                                                                                                                                     acaAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCgt   >  6:7370/1‑143 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTT  >  NC_002947/3688051‑3688341

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