Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,718,463	T→G	34.8%	Q65P (CAG→CCG)	CJ019_01614 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,718,463	0	T	G	34.8%	14.1 / 3.6	27	Q65P (CAG→CCG)	CJ019_01614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/13);  new base G (5/5);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTG  >  ACB122/1718333‑1718591                                                                                                                                   |                                                                                                                                 aTTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGgcgc                                                                                                                             >  3:857/1‑136 (MQ=255)      tgaACCGCAACTTGCCATTCAGCTCATGGCGCGCACTTCCCGCATAGCGGGGGCCTGGGGTGCCCGGCGCGGGGTCCTCGGTCGGTTGGACCGGCTGCGGCCCGTGGCCTTTCACCGGGCGGGCCTGGCGCGCCAg                                                                                                                        <  6:18854/136‑1 (MQ=255)        acccgCAACTTGCCATGCCGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTGGCCCTGCGCGGGGTCATAGGTGAGAGGGAACGGCCGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                        <  1:43676/130‑1 (MQ=255)             cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGGGCGCCCggggcgcg                                                                                                                 >  4:189529/1‑131 (MQ=255)             cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACg                                                                                                                 >  4:135193/1‑136 (MQ=255)             cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCAGGCGCGCCCGGGTcgcc                                                                                                                 >  1:133748/1‑133 (MQ=255)                cTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACgg                                                                                                                >  6:88281/1‑134 (MQ=255)                               aaTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCa                                                                                                 <  1:3945/134‑1 (MQ=255)                               aaTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                               <  1:146/136‑1 (MQ=255)                                aTGCGCGCACTTACCGCATAGCGGTGGCCCGGGTTGACGTGCGCGGGGTCATAGGTGAGCTGGACCGGCCGCGGCAGGTAGCCTTTGCCCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                              <  3:189529/136‑1 (MQ=255)                                   cgcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                            <  2:145300/135‑1 (MQ=255)                                              cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGt                                                                                <  1:68538/136‑1 (MQ=255)                                                    gCGGGGGCCCGGGTTGCCCTGGGGGGGGGCATAGGGGAGATGCACCGGCAGCGGCAGGTTGCCTTTGCCCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                           <  1:192236/135‑1 (MQ=255)                                                       gTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCCTCAGCGACACCTCCTGCAGTTCCCCGcg                                                                        >  4:57355/1‑134 (MQ=255)                                                                  tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGGGCCAGGGTCACCGCCGGGGCGTGCAGCAGCGACCCCACCATCAGGTCCACGCTCAGGGTggcggc                                                            >  1:74718/1‑136 (MQ=255)                                                                                   tcggcgagaTGCCACGGCAGCGGCACGTGGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgc                                           <  3:135193/131‑1 (MQ=255)                                                                                     ggTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                         <  7:32241/136‑1 (MQ=255)                                                                                     ggTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                         <  4:25326/136‑1 (MQ=255)                                                                                     ggTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                         <  4:857/136‑1 (MQ=255)                                                                                     ggTGAGATGGAACGGCAGCGGCACGTGGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                         <  3:100330/136‑1 (MQ=255)                                                                                      ggtcgaTGCAACGGCAGCGGCACGTGGCGTTTGCCCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                         <  6:4831/131‑1 (MQ=255)                                                                                            tGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGGCCGGGTGACCGGCGGGGCGGGCCTCAGCCACACATCCTGAAGTTGCACGCTCACGGGGGCGGGAGGCGGCAGGGCAATGCGCGGCAGGt                                  >  4:199525/1‑136 (MQ=255)                                                                                                   gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGTGCGCGGCAGGGGCACGGGAGGGGGGTCCCGCACCGACCCATACTGCGCTTCCACGCTCAGCGTGGCGGCCGGCGGGCGGGCCAAGCGCTGCAGGGAGAAAAc                           >  2:165309/1‑136 (MQ=255)                                                                                                    cAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                          <  2:133748/136‑1 (MQ=255)                                                                                                             gTTGCCTGTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTcc                 <  8:115399/136‑1 (MQ=255)                                                                                                                     ggcccGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCCGGCTGGc         <  3:18228/133‑1 (MQ=255)                                                                                                                            gccggccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCTGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg  <  2:74718/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 ATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTG  >  ACB122/1718333‑1718591

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