Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,146,403	C→A	19.2%	R110S (CGC→AGC)	PP_2757 →	sugar‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,146,403	0	C	A	19.2%	54.7 / 2.4	20	R110S (CGC→AGC)	PP_2757	sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base A (2/2);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCGGGCGCAACGACCGCAATCTGGTCACACAACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCGCGTTGGCC  >  NC_002947/3146286‑3146507                                                                                                                      |                                                                                                         gACGCCGGGCGCAACGACCGCAATCTGGTCACACAACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTa                                                                                        <  1:237873/136‑1 (MQ=255) gACGCCGGGAGGAACGACCGCATTCTGGTCTCAGAACTGCACAGGGTGGTAACGAAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGTCCGCTGGCTCAAGCGCTTCCTGGCAGCAGGTa                                                                                        <  1:306668/136‑1 (MQ=255)           gcAACGACCGCAATCTGGTCACACAACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCtttt                                                                              <  1:39901/136‑1 (MQ=255)                   cGCAATCTGGTCACACAACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTg                                                                      <  1:331274/136‑1 (MQ=255)                                cacaACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCa                                                         <  2:18816/136‑1 (MQ=255)                                              gtggtAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTcaaca                                           <  1:76444/136‑1 (MQ=255)                                              gtggtAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTcaaca                                           <  1:352707/136‑1 (MQ=255)                                                  tAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGAGCCATCTCAGTACAGCCTCAACAacac                                       >  1:340720/1‑136 (MQ=255)                                                     cGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGAGCATCCAGGCAGCAGGt                                                                                         <  1:372326/83‑1 (MQ=255)                                                     cGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGAGCATCCAGGCAGCAGGt                                                                                         <  1:374342/83‑1 (MQ=255)                                                     cGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGAGCATCCAGGCAGCAGGt                                                                                         >  2:372326/1‑83 (MQ=255)                                                     cGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGAGCATCCAGGCAGCAGGt                                                                                         >  2:374342/1‑83 (MQ=255)                                                            gCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGAc                             >  2:213138/1‑136 (MQ=255)                                                            gCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGAc                             >  2:262665/1‑136 (MQ=255)                                                            gCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGAc                             >  2:81816/1‑136 (MQ=255)                                                              ccGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACaa                           >  2:212657/1‑136 (MQ=255)                                                                                  cTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGTTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCcgcgcgg       >  1:229219/1‑135 (MQ=255)                                                                                       tACCCTGTGCGTGATCGCCCCGTGGCTCAAGCGCATCCAGGCAGCAGGTACCCCGCTTTTCACCATTGATGTGCCATCGCATAACCGCCTCAACAACACCCCCTCGGACAACGTCGCGGCCGGCCGCGCGTTGGcc  <  2:229035/136‑1 (MQ=255)                                                                                       tACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCGCGGTGGcc  >  1:287188/1‑136 (MQ=255)                                                                                       tACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCCGCCGCGCGTTGGcc  >  2:228492/1‑136 (MQ=255)                                                                                                                      |                                                                                                         GACGCCGGGCGCAACGACCGCAATCTGGTCACACAACTGCACACGGTGGTAACGCAAAAGCCCGATGCGGTAATTCAGACCCTCGGTACCCTGAGCGTGATCGACCCCTGGCTCAAGCGCATCCAGGCAGCAGGTATCCCGCTTTTCACCATTGATGTGCCATCGCAGTACAGCCTCAACAACACCACCTCGGACAACGTCGCGGCCGGCCGCGCGTTGGCC  >  NC_002947/3146286‑3146507

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