Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,094,817	T→G	60.0%	Q183P (CAG→CCG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,094,817	0	T	G	60.0%	‑5.2 / 2.9	20	Q183P (CAG→CCG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base G (6/5);  total (13/7)
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.06e-01

CTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB122/2094702‑2094936                                                                                                                    |                                                                                                                        cTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATATGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCGGGCCGTTTTg                                                                                                      >  8:96443/1‑135 (MQ=255)      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCCAGCCCGCCCGGCCCGGCGGGGCGTTTTTGAggg                                                                                                 >  8:129301/1‑135 (MQ=255)      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGGCGTTTTTGAggg                                                                                                 >  4:133090/1‑135 (MQ=255)      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCTGGGCGTTTTGCAggg                                                                                                 >  6:88580/1‑135 (MQ=255)             ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGCCCCCCCGCGCCGGGCCGGCCAGTTTTGAGGGccggggg                                                                                          >  7:286346/1‑130 (MQ=255)                        gTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCCCCCCGGGCCGGGCAGGCCGTTTTGCACGGCCTGGTGGCCGAGcttgt                                                                               >  8:26607/1‑135 (MQ=255)                            ctcAACATCCGCCAGCTGCTCTTCAACCCCCTCGGGCGCCGGCTGGTGGCGAAAATCCTGCACGCCCCGCTCGACGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                            <  4:32776/134‑1 (MQ=255)                              cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCAGGCCGGGCGGGCCTTTTTGGCGGGCCTGGTGGCCGAGGTTGTTgaaga                                                                         >  3:225528/1‑135 (MQ=255)                                        cAGCTCCCCTCCATACCGCTGCGGCGCCGGCTGGTGCCGAAAATACTGCCCGCCGCGCGCCTCGGCCAGCCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                               <  7:32166/135‑1 (MQ=255)                                         aGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGGCTGTTGGGCGCGGGTGTTGAAGTACTCTatattc                                                              >  2:91768/1‑133 (MQ=255)                                          cctgctCTCCAACCCCGTTCGGCTCGGGCTGCTTCCCAAAAAACCCCCCGCCGCGCGCGCCGGCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgca                                                             <  6:283869/134‑1 (MQ=255)                                              ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTTCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATgcatcac                                                         >  8:249079/1‑134 (MQ=255)                                                      accccTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGGTGGCCGAGGTTGTTGAAGAACTCGATATGCCGAAGCACCAgg                                                  >  6:299298/4‑134 (MQ=255)                                                      acacCTTCGGCTCCGGCTGATTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGGCGTGTTGCAGGGCCTGTTGGCCGACGGTGTTGAACAACTCTATATGGAGACGCACCAgg                                                  >  7:53613/1‑134 (MQ=255)                                                           ttCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTTCAGGGCCGGTTGGCCGCGGTTGTTGTATAACTCTATATGCAGCAGCCCCCGGAGATgg                                            >  2:20513/1‑135 (MQ=255)                                                              ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCGGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCATCAGCACCAGGATATGGc                                           >  8:161844/1‑133 (MQ=255)                                                              ggCGCCGGCGGCTGGCGAAAATACTCCCCGCCGCGCGCGACGGCCCGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                           <  2:242385/133‑1 (MQ=255)                                                                                  attcTGCCCGCCGCGCTCGTAGGCCTGGGGGGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCtt                     <  7:96443/132‑1 (MQ=255)                                                                                                     ccggCCCGGCGCGGCGCGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  6:2352/133‑1 (MQ=255)                                                                                                     acggCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  7:248696/133‑1 (MQ=255)                                                                                                                    |                                                                                                                        CTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB122/2094702‑2094936

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