Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,754,185	G→T	22.2%	W68L (TGG→TTG)	PP_2405 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,754,185	0	G	T	22.2%	41.8 / 2.8	18	W68L (TGG→TTG)	PP_2405	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/10);  new base T (2/2);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GCTACGGCGCCATGCCATCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCATGATGTGGCGACTTTGCG  >  NC_002947/2754077‑2754317                                                                                                             |                                                                                                                                     gCTACGGCGCCATGCCATCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCt                                                                                                           <  2:36780/136‑1 (MQ=255)   ttcgtcgcCATGCCATAGCTGGCGGGGGCAAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCGCCGGGGCTCAATGTCAAGGTGGTGCTGGTGGTCTTCTGGAACTGGGACTTCATCAACTTCCGTCGCATCCTGc                                                                                                         <  2:251135/131‑1 (MQ=255)          ccatgccatCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTTGGACTGCATCAACTGCCGGCGCAGCCTGCCGCAtgt                                                                                                  >  2:134150/1‑136 (MQ=255)          ccatgccatCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTTGGACTGCATCAACTGCCGGCGCAGCCTGCCGCAtgt                                                                                                  >  2:1483/1‑136 (MQ=255)               ccatCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATg                                                                                             <  2:397001/136‑1 (MQ=255)                    gCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTgg                                                                                        <  2:264952/136‑1 (MQ=255)                         ggggcgCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGggcgcg                                                                                   <  1:399510/135‑1 (MQ=255)                                     agTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGa                                                                       <  2:490426/136‑1 (MQ=255)                                           tGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAgg                                                                 <  2:356288/136‑1 (MQ=255)                                             aaTTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGc                                                               <  1:20355/136‑1 (MQ=255)                                              aTTCGCCGCCCCTTGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGTGACTGCATCAACTGCCGGTGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCt                                                              <  2:347426/136‑1 (MQ=255)                                                                    ggggCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACAcc                                        >  1:112536/1‑136 (MQ=255)                                                                                     tggtgCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGa                       >  2:60314/1‑136 (MQ=255)                                                                                       gtgCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCTACTACGAGc                     <  1:427267/136‑1 (MQ=255)                                                                                           tggtggACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCatga                 >  2:386315/1‑136 (MQ=255)                                                                                                    tCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTCAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCATGATGTGGCGAc        >  1:420519/1‑136 (MQ=255)                                                                                                          ccTTGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCATGATGTGGCGACTTTGCg  <  1:1483/136‑1 (MQ=255)                                                                                                          ccTTGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCATGATGTGGCGACTTTGCg  <  1:134150/136‑1 (MQ=255)                                                                                                             |                                                                                                                                     GCTACGGCGCCATGCCATCGCTGGCGGGCGCCAGCCAGTGGTTGAATTCGCCGCCCCTGGATGCCCCGGGGCTCAAGGGCAAGGTGGTGCTGGTGGACTTCTGGACCTGGGACTGCATCAACTGCCGGCGCAGCCTGCCGCATGTGAATGATTGGGCGCGGCGCTATGCCGATCAAGGGCTGGTGGTGGTAGGGGTGCACACCCCTGAGTTCGACTACGAGCATGATGTGGCGACTTTGCG  >  NC_002947/2754077‑2754317

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