Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,363,720	C→T	13.8%	A32T (GCC→ACC)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,363,720	0	C	T	13.8%	62.9 / 4.1	29	A32T (GCC→ACC)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (15/10);  new base T (2/2);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATG  >  NC_002947/5363596‑5363860                                                                                                                             |                                                                                                                                             aTCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAg                                                                                                                     >  2:70711/1‑150 (MQ=255)       tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTc                                                                                                                   <  1:249766/146‑1 (MQ=255)       tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTc                                                                                                                   >  2:249766/1‑146 (MQ=255)       tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                                                                                                               <  1:70711/150‑1 (MQ=255)       tCGGCACCTGCCGCGGACAGGCCGACTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc                                                                                                               <  2:589804/150‑1 (MQ=255)               tGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCgg                                                                                                       <  2:81584/150‑1 (MQ=255)                         aGGCCGGCTTCGAGTGCGGACTCGAACATGTAGACGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCg                                                                                             >  1:557653/1‑150 (MQ=255)                            ccGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTTCGGAAGGCCATGCCGGCCGCCCAGCCAAGCCTCCGCATGAAGTCAGGCGTGATCGGGTATTCGCCGAcg                                                                                          >  1:409241/1‑150 (MQ=255)                             cGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGAcgc                                                                                         <  1:209434/150‑1 (MQ=255)                                    gAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACg                                                                                  >  1:254141/1‑150 (MQ=255)                                                                  aCGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCt                                                    <  2:114187/150‑1 (MQ=255)                                                                    gcgTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAGAGTATTTTCTGCTCa                                                  >  1:197093/1‑150 (MQ=255)                                                                            ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCgg                                                                       >  1:24947/1‑121 (MQ=255)                                                                            ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCgg                                                                       <  2:24947/121‑1 (MQ=255)                                                                            ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACt                                          >  2:229559/1‑150 (MQ=255)                                                                               cccACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCa                                       >  1:602362/1‑150 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg                                                                           >  1:576364/1‑109 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg                                                                           >  1:576326/1‑109 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg                                                                           <  2:576326/109‑1 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg                                                                           <  2:576364/109‑1 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGc                                                                                      >  1:551022/1‑98 (MQ=255)                                                                                    caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGGATTCGCCGACGCGGc                                                                                      <  2:551022/98‑1 (MQ=255)                                                                                                          tagccagaGGCTATGGAGGCCGACCAGCTAATCGTCAGCATGGAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGc            >  2:499997/8‑150 (MQ=255)                                                                                                           tgcGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCa           <  1:229559/150‑1 (MQ=255)                                                                                                               gAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGa       >  1:504842/1‑150 (MQ=255)                                                                                                                aaGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc                                 <  1:96999/123‑1 (MQ=255)                                                                                                                aaGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc                                 >  2:96999/1‑123 (MQ=255)                                                                                                                    ccATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATg  >  1:284408/1‑150 (MQ=255)                                                                                                                    ccATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATg  >  1:260786/1‑150 (MQ=255)                                                                                                                             |                                                                                                                                             ATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATG  >  NC_002947/5363596‑5363860

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