Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,721,870	A→C	30.9%	V98G (GTC→GGC)	CJ019_01617 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,721,870	0	A	C	30.9%	21.7 / 2.4	33	V98G (GTC→GGC)	CJ019_01617	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/2);  new base C (7/4);  total (27/6)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCC  >  ACB122/1721741‑1721998                                                                                                                                  |                                                                                                                                 tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATTCGGACCCAGCGGTGAGCATGATATGCCCCGTCTTCCGAGCCCCAg                                                                                                                             >  5:25611/1‑135 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTGCGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                             >  8:448408/1‑135 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGcccccg                                                                                                                             >  7:372506/1‑133 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                             >  1:269877/1‑135 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                             >  1:44838/1‑135 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCACCGTCGTCCTCGcccccg                                                                                                                             >  7:256189/1‑133 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCACCGGTCCGCATCGTCGTCCCCGTCTTCCACGcccccg                                                                                                                             >  8:180453/1‑133 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACAGTCGGCGCTGATGCGGACCCCGCGGTCCGCATCATCGTCCCCGTCATCCGAGACCCAg                                                                                                                             >  5:122899/1‑135 (MQ=255)     gggTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCAc                                                                                                                         <  8:235631/135‑1 (MQ=255)       gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCCGCATCTTCTTCCCCGTCTTCCTCGCCCCAGcccccc                                                                                                                       >  7:332077/1‑135 (MQ=255)       gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGcccccc                                                                                                                       >  2:60221/1‑135 (MQ=255)       gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCAccc                                                                                                                       >  2:13465/1‑135 (MQ=255)                  ccTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCCGCATCTTCTTCCCCGTCTTCCTCGACCCACCCCCCCCgcgcgccccc                                                                                                            >  6:275441/1‑131 (MQ=255)                    tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGGCCCAGCGGGCAGCATCTTCTTACCCGTCGTCCTCGCCCCCGCCCCCCGGcgcgccgcgc                                                                                                           >  1:154924/1‑134 (MQ=255)                    tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCgcgcagcgc                                                                                                           >  3:149132/1‑134 (MQ=255)                    tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCCGCgcgcagcgca                                                                                                          >  1:120772/1‑135 (MQ=255)                         aaaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCCCGGTCAGCATCTTCTTCACCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCCCATcc                                                                                                     >  2:250656/1‑135 (MQ=255)                         aaaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCGTCCTCGACCCCGCCCCCCGGCGCGCCGCGCCCGTcc                                                                                                     >  3:171641/1‑135 (MQ=255)                          aaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCAGATGCGGACCCAGCGGTCAGCATCTTCTTACCCCGCTTCCTCGACCCAGCCACCCGGCGCGCACCCCACATCca                                                                                                    >  3:77347/1‑135 (MQ=255)                              aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCccac                                                                                                 >  4:430677/1‑134 (MQ=255)                                                     tACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCATCATCCCCGTCATACTCGCCCCAGCCACCCCGCGCGCCGCGCCCTTCCACACCCACCCCCGCAAATGAGGCCCggg                                                                         >  2:109518/1‑135 (MQ=255)                                                        gCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCCTCTTCCTCCACCCCGCCACCCCGCGCGCCGCGCCCTTCCACACCCAGCGCGGCAAATGCGGCGCGGgc                                                                        >  8:72596/1‑133 (MQ=255)                                                        gCCAGGTGATTTCCTCTTCACCGTCGGCGCAGATGCGGACCCAGCGGTCAGCATCTTGTTCCCCGTCTTCCTCGCCCCACCCCCCCGGCGCGCAGCGCCCTTGCACACCCACCGCGGCAAATGCGgcgcgcgc                                                                        >  2:308189/1‑133 (MQ=255)                                                                aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCCCTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACg                                                              >  8:305033/1‑135 (MQ=255)                                                                          cACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCCGCGCGCCGCGCACTTCCCCACCCAGCGCGGCAAAATCGGCGCGGGCGCCGGCCACGTCATCGGccc                                                    >  4:355167/1‑135 (MQ=255)                                                                          cACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCGTCTTCCCCGTCTTCCTCGCCCCAGCCCCCCGGCCCGCCCCGCACGTCCACCCCCCGCGCGGCCAATGCGGCGCGGGCGCCGGCCACGGCATCGGccc                                                    >  6:39635/1‑135 (MQ=255)                                                                             cGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACg                                                 >  7:77733/1‑135 (MQ=255)                                                                                     cTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCGCGCCGCGCACTTCCACACCCACCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCCTCGGCCCCCcgcggcggg                                         >  5:466022/1‑130 (MQ=255)                                                                                                       tCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAGCCACCCGGCGCGCAGCGCCCTTCCACCCCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTCGGGGCGCGTTCCCGCTACAgg                       >  2:427771/1‑135 (MQ=255)                                                                                                                  ccccccGTCTCCCCCCCCCCCCCCCCCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCa               <  3:446372/130‑1 (MQ=255)                                                                                                                    ccgtgTCTGCCGCGCCCCCGCCACCGGGGGCTCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc          <  3:355167/131‑1 (MQ=255)                                                                                                                     ccggTCTTCCTCCCCCCCCCCACCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc          <  4:149132/131‑1 (MQ=255)                                                                                                                     cccGTCTGCCTCCACCCAGCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc          <  8:392972/134‑1 (MQ=255)                                                                                                                     cccGTCTCCCGCCCCCCCCCCCCCCGGCGCGCCCGCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGTTTTCCAGATACAGGCTGTTCCATTTcc          <  6:375133/134‑1 (MQ=255)                                                                                                                           ttCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGc   >  5:179144/1‑135 (MQ=255)                                                                                                                            tCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGcc  >  6:202137/1‑135 (MQ=255)                                                                                                                                  |                                                                                                                                 TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCC  >  ACB122/1721741‑1721998

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