Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,626,701	G→A	17.4%	N195N (AAC→AAT)	PP_4940 ←	glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,626,701	0	G	A	17.4%	52.8 / 4.7	23	N195N (AAC→AAT)	PP_4940	glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base A (2/2);  total (10/13)
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.48e-01

TTGCAGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGCGGCCGCTGTAGGC  >  NC_002947/5626566‑5626802                                                                                                                                        |                                                                                                      ttGCAGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCGTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACgg                                                                                                       >  2:9116/1‑136 (MQ=255)   gCAGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGtt                                                                                                     >  2:414531/1‑136 (MQ=255)     aGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTaa                                                                                                   <  1:190457/136‑1 (MQ=255)     aGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTaa                                                                                                   <  2:394570/136‑1 (MQ=255)     aGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTaa                                                                                                   <  1:441556/136‑1 (MQ=255)     aGACCCAGCGCACGGTTCTCTTCCTAGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGACATGGCTTCGAAGTTGTTGATCATGTATTCGCCCATGCCGACATCGGCACTGTCTAAGAAACGAATACGGTTaa                                                                                                   <  1:412769/136‑1 (MQ=255)       aCCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGc                                                                                                 <  1:309473/136‑1 (MQ=255)                                   tGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTgg                                                                     <  1:414531/136‑1 (MQ=255)                                   tGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTgg                                                                     <  1:465602/136‑1 (MQ=255)                                        gTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACgg                                                                >  2:451561/1‑136 (MQ=255)                                              cAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGATTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTAcc                                                          >  2:310266/1‑136 (MQ=255)                                              cAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGATTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTAcc                                                          >  2:310493/1‑136 (MQ=255)                                                         aGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGATTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTtctc                                               <  1:310493/136‑1 (MQ=255)                                                         aGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGATTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTtctc                                               <  1:310266/136‑1 (MQ=255)                                                                 gCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGcc                                       >  1:230870/1‑136 (MQ=255)                                                                   cATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGcccc                                     <  1:297000/136‑1 (MQ=255)                                                                        cTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTc                                >  1:363682/1‑136 (MQ=255)                                                                              aGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCcag                          >  1:399030/1‑136 (MQ=255)                                                                              aGTTCTTGATCATGTATTCGCCCATGCCGACAAAGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCcag                          >  1:160463/1‑136 (MQ=255)                                                                                     gATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGcc                   >  2:441654/1‑136 (MQ=255)                                                                                         aTGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGc               <  1:9116/136‑1 (MQ=255)                                                                                         aTGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGCACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGc               <  1:142442/136‑1 (MQ=255)                                                                                                         cGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGCGGCCGCTGTAGGc  >  1:367935/1‑133 (MQ=255)                                                                                                         cGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGCGGCCGCTGTAGGc  <  2:367935/133‑1 (MQ=255)                                                                                                                                        |                                                                                                      TTGCAGACCCAGCGCACGGTTCTCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGTTAAGCGTGTTGCAGTAGTCCTCGCCCGACTTGGTACGGGTTACCAACAGGTTCTCGACCTGCCCCAACTCGTCCAGCAGTGCCTTGCGGCCGCTGTAGGC  >  NC_002947/5626566‑5626802

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