Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,377,746	C→A	13.9%	A176S (GCG→TCG)	grpE ←	heat shock protein GrpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,377,746	0	C	A	13.9%	70.2 / 4.1	29	A176S (GCG→TCG)	grpE	heat shock protein GrpE
Reads supporting (aligned to +/- strand):  ref base C (15/10);  new base A (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 = 7.59e-01

CCATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCA  >  NC_002947/5377609‑5377885                                                                                                                                          |                                                                                                                                            ccATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCtt                                                                                                                                 >  2:106404/1‑150 (MQ=255) ccATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCtt                                                                                                                                 >  2:106389/1‑150 (MQ=255)           tctcCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGaccacc                                                                                                                       >  1:255911/1‑150 (MQ=255)                cGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCAt                                                                                                                  >  2:131012/1‑150 (MQ=255)                    aCTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCg                                                                                                              <  1:106583/150‑1 (MQ=255)                       tGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCgggg                                                                                                           >  1:217717/1‑150 (MQ=255)                           tttGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGCCgca                                                                                                       >  2:398208/1‑150 (MQ=255)                           tttGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGAGCGCGgcgcaacc                                                                                                                                >  1:106543/1‑121 (MQ=255)                           tttGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGAGCGCGgcgcaacc                                                                                                                                <  2:106543/125‑5 (MQ=255)                           tttGGTAGCAGCGACGACTTTGGCCAACTGCGGTAATACTTGAAGGCTTCACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTGATTAATCGTGGGTCGCGCCGCGCTCGGCGCCTTCCTGACCACCACCATGGCGGGgcgcc                                                                                                       >  2:93310/1‑149 (MQ=255)                                                ggCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGt                                                                                  <  2:103722/150‑1 (MQ=255)                                                          ggTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCAAGGCCAGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTg                                                                        >  1:48306/1‑150 (MQ=255)                                                                    gAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTc                                                              >  2:140904/1‑150 (MQ=255)                                                                        gCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAaca                                                          >  1:10396/1‑150 (MQ=255)                                                                        gCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGAAACACGTTCAaca                                                          >  1:81322/1‑150 (MQ=255)                                                                          ttGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTc                                                                          <  2:34658/132‑1 (MQ=255)                                                                          ttGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTc                                                                          >  1:34658/1‑132 (MQ=255)                                                                                ccTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGttt                                                  <  1:106404/150‑1 (MQ=255)                                                                                ccTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGttt                                                  <  2:217717/150‑1 (MQ=255)                                                                                ccTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGttt                                                  <  1:106389/150‑1 (MQ=255)                                                                                 cTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTg                                                 <  1:131012/150‑1 (MQ=255)                                                                                     aTGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTc                                             <  2:81322/150‑1 (MQ=255)                                                                                                           cAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTTCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAAACCTTGTGGAACACGGTCAACTTGCAGTTTGGCTCGAACGCGGCGTTTTCCTGcatg                       >  2:140439/1‑150 (MQ=255)                                                                                                                          aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGtt                                                               >  1:236492/1‑95 (MQ=255)                                                                                                                          aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGtt                                                               <  2:236492/95‑1 (MQ=255)                                                                                                                           tCGAGGGTGTTGCCGCGCGCGGCGCCTTGCTGACCACCACCGTGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTggtggt       >  2:130196/1‑150 (MQ=255)                                                                                                                                ggTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCa  <  2:255911/150‑1 (MQ=255)                                                                                                                                   ggcgcCGAGCGCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGc                                    <  1:268841/112‑1 (MQ=255)                                                                                                                                   ggcgcCGAGCGCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGc                                    >  2:268841/2‑113 (MQ=255)                                                                                                                                          |                                                                                                                                            CCATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCA  >  NC_002947/5377609‑5377885

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