Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,969,760	A→C	51.6%	A251A (GCT→GCG)	CJ019_03750 ←	putative protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,969,760	0	A	C	51.6%	‑1.2 / 5.5	15	A251A (GCT→GCG)	CJ019_03750	putative protein
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (3/5);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCA‑GGGCAATGTCGGCCAAGGCCTGGCGC  >  ACB122/3969635‑3969874                                                                                                                              |                                                                                                                    ggCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCCGGCACAGTGGGCGCgccccgcgcgcagc                                                                                                            >  4:55840/1‑135 (MQ=255) ggCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCgcgccgcgcgccgc                                                                                                            >  8:252062/1‑135 (MQ=255)          ggcggGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAgg                                                                                                   >  8:287000/1‑135 (MQ=255)                 gggcggggggCGACATCCGGCGGCGTGGGCCCTGGGCGAATGCGCCCACCGCCGCGACACCCACTTCTGCCATCGCCCCCCCCCCCCGCCCGCCACAGGGGGGGCGCGCCGCGCGCAGCCACACCAGGTCGgcgc                                                                                            <  2:392349/129‑1 (MQ=255)                               aTACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACg                                                                              >  3:429271/1‑135 (MQ=255)                                tccggcggCGAGGGCGCCGGGCTAATGCGGCCCCCGCCTCGACACCCACTTCTTCGAACGCGGCACCCCCCCGCCCGGCCCCGGGGGCGCGCGCCGCGCGCAGCCACACCAGGGCGGCGCGCAGGCCGTCCACGc                                                                             <  4:429812/133‑1 (MQ=255)                                tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCAGCCACACCACGTAGGCGCGCAGGCCGTCCCCGc                                                                             >  6:365171/1‑135 (MQ=255)                                aaCGGCGGCGTGGGCGCTGGGCTAATCCGCCCCCCCCCGCCACACCCCCTTCTACGAACGCGCCCCCCCCCCCCCCGCCCCCGGGGGCGCGCGCCGCGCCCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                             <  8:132978/134‑1 (MQ=255)                                       gcgATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCCGCCCCCCCAGGTCGGCGCGCAGGCCGTCCCCGCCggcggc                                                                      >  1:38861/1‑135 (MQ=255)                                                  gggCTAATCCGCCCCCCGCCGCCACAACCACTTCTGCTATCGCGCCCCCCCCCCGCCCGCCACTGGGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCt                                                           <  6:182934/135‑1 (MQ=255)                                                                 ccgcggcGACACCCACTTCTTCGATAGCTGCCGCCCCCCGCCTGGCACTGGGGGCGCGCGCAGCGCGCAGCCACACCAGGGCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCaa                                            <  7:287000/129‑1 (MQ=255)                                                                           acccACTTCTTCTATCGCGGCCCCCCCCCCCCCGGCACCGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTc                                  <  4:369575/133‑1 (MQ=255)                                                                               aCTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCCCGCAGCCCCCCCACGTCGGCGCGCACGCCGTCCACGCCGGCGGCCTAGCAGCGCGCGGTAATACACGCCAGCGCCAGGTCAtac                              >  4:60993/1‑133 (MQ=255)                                                                                     tcGATCGCGGCCCCCCCCCCCCCGCCACCGGGGGCGCGCCCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCA‑GGGCa                       <  4:298245/135‑1 (MQ=255)                                                                                       gATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGACCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCA‑GGGCAAt                     >  5:444460/1‑135 (MQ=255)                                                                                              gCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCCCGGAGGCGCCCCCGCCCGCCACGCCCGCGGCGTTGCCGCGCGCGGGTATACCCGCAAGCGCCGGGTGATGCACGGGCaagtaggc                >  8:254073/1‑128 (MQ=255)                                                                                                          cAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCACCAGGTCGGCGCGCAGGCCGTCCACGCCCGCGGCGGAGCCGCGCTGGGTAATCCAGGCCAGCGCCGGGGGCTCCC‑GGGCCATGTCGGCCAACGCCCGgcgc  >  8:355125/1‑135 (MQ=255)                                                                                                                              |                                                                                                                    GGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCA‑GGGCAATGTCGGCCAAGGCCTGGCGC  >  ACB122/3969635‑3969874

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