Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,718,463	T→G	32.9%	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	32.9%	30.0 / 2.3	42	Q65P (CAG→CCG)	CJ019_01614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base G (12/3);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 9.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGG  >  ACB122/1718338‑1718593                                                                                                                              |                                                                                                                                   tgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                          <  8:81136/136‑1 (MQ=255)  caaCAGCAACTCGCCATTCAGCTCACTGCGCGCACTTACCGCATAGCGGGGGCCCGGCTGGACCTGCGCGGGGTCAGAGGTGAGAGGGAACGGCCGCGGCACGTTGCCGTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                          <  3:335929/134‑1 (MQ=255)   aaCAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCgg                                                                                                                         >  5:148502/1‑135 (MQ=255)   aaCAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGgcgcggcccg                                                                                                                         >  1:39233/1‑130 (MQ=255)    ccAGCAACTTGCCATTCTGCTCAATGCGCGCACTTCCGGCATAGCGGTGGCCAGGCTTGACCTGCGGGGGGTCATAGGTGAGATGGAACGGCCGCGGCAGGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAggg                                                                                                                        <  2:257261/134‑1 (MQ=255)             tGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCACGGTGACGGCAg                                                                                                               >  2:396090/1‑135 (MQ=255)                  ttCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCCCGGCAGGGGCg                                                                                                          >  2:78993/1‑135 (MQ=255)                          aaTGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGGCCGGGTCCCGGCAGGGGCGTGCAACAg                                                                                                  >  5:278655/1‑135 (MQ=255)                            tGCGCGCACTTCCCGAATTGCGGTGGCCAGTCTTGACCTGCGGGGGGTCTTGGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                                <  1:131552/135‑1 (MQ=255)                               gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCCCCAGGGGCACGGCAGGGGCGTGCATCAGCGaca                                                                                             >  2:370536/1‑135 (MQ=255)                               gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCCGGGGCGGCCCTAAACGacc                                                                                             >  2:131552/1‑134 (MQ=255)                               gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCACGGCCGGGGCGTCCATCAGCGccac                                                                                            >  7:295451/1‑136 (MQ=255)                                   aCTTCCCGCATAGCGGGGGCCAGGCTTGACCTGCGGGGGGTCTTAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                         <  6:78396/135‑1 (MQ=255)                                    cTTACCGCATAGCGGTGGCCAGGCTTGCCCTGCGCGGGGACATAGGTGAGATGGCACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTAGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGAGACACATCCt                                                                                       >  7:32197/1‑136 (MQ=255)                                    cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGGCCGGGGCCCGGGCGGGGCGGCCCTCAGCgccccct                                                                                          >  4:354117/1‑127 (MQ=255)                                         ggCATTGCGGTGGGCTGGATTGAGAGGCGCGGGTGAATCGGTGAGATTGACGGTCCGCGGCAGGTGGGCTTTGCCCGGGCAGCGCTTGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGt                                                                                  <  8:346687/135‑1 (MQ=255)                                         cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCCCGGCAGGGGCGTCCATCAGCGGCACATCCTGCCGt                                                                                  >  8:286264/1‑136 (MQ=255)                                          gCATCGCGGGGGCCAGGCGTGACCTGCGGGGGGTCATAGGGGCGATGGACCGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGt                                                                                  <  6:53248/135‑1 (MQ=255)                                                cGGTGGCCAGGCTTGCCCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCCCGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACCGCAGGGGCGTCCCTCCACGCCCCAACCAGCAGttgcc                                                                              >  6:159830/1‑130 (MQ=255)                                                cGGTGGCCAGGCTTGACCTGCGCGGGCTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCAGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCAACACCCACACATCCAGCAGTTCCACGc                                                                           >  7:346687/1‑136 (MQ=255)                                                  gTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGGCCATCAGCGCCCCCTCCTGCAGTTCCACGCt                                                                          >  4:66322/1‑135 (MQ=255)                                                      ccAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCGGCACCTCCTGCAGTTCCCCGCGCAggg                                                                     >  7:188612/1‑136 (MQ=255)                                                            ttcaCCTGCGCGGGGTCATAGGTGAGATGTAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                <  1:370991/132‑1 (MQ=255)                                                             tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGCGCGCGCCAGGGGCCCGGCACGGGCGGCCCGCAGCGCCCCCACCGGCCGGTACACGCTCCCGGGggcgg                                                               >  1:108519/1‑135 (MQ=255)                                                                 cggcgcGGGGCCATCGGTGCGAGGGAACGGCAGCGGCACGTGGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                           <  5:386290/133‑1 (MQ=255)                                                                 cTGCGCGGGGGCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                           <  4:93202/135‑1 (MQ=255)                                                                  tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCCGGGGCGTCCAGCAGCGACACCTCCTGCAGTTCCACGCGCCGGGTggcggccggc                                                          >  4:9898/1‑135 (MQ=255)                                                                       ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCCGGTCCCGGCCGGGGCGTGCATCAGCGACCCCTCCAGCCGTTCCACCCGCAGGGTGgcggcaggcggcag                                                     >  6:74018/1‑135 (MQ=255)                                                                               aGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgc                                             <  6:148502/135‑1 (MQ=255)                                                                                 gTGAGATGGAACGGCAGCGGCACATTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACAACCTGCAGTTCCACGCGCAGGGTGGCGGCAGGCCGCAGGGCAATgcgc                                           >  6:205861/1‑135 (MQ=255)                                                                                      aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCCGGGGCGGCCATCAGCGCCCCATCCTGCCGTTCCACGCTCCGGGTGGCGGCCGGCGGCAGGGCAATTGGCTGCAg                                      >  4:134409/1‑135 (MQ=255)                                                                                       tGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCGAGGCGCGGCAGGGTCACCCGCGGGGCGTGCCACATCGCCCCCACCAGCACATCAACGCGCAGGGTGGCGGCGGGGGGGCGGGCCAATCGCTGCAggg                                    >  8:125880/1‑135 (MQ=255)                                                                                         gAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTa                                   <  8:359972/135‑1 (MQ=255)                                                                                                aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                            <  2:220314/135‑1 (MQ=255)                                                                                                       cGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGGGACGGCCCGGGCGGGCCTCAGCGACACCAACTGCCGTTGCACGCTAAGGGTGGGGGCCGGCGGCAGGGCAAGGGGCGGCAGGGTGAAGACGTCAccctg                    >  8:208685/1‑133 (MQ=255)                                                                                                       cGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGGCTCGGCCGGGGCGTGCATCAGCGCCACAACCTGCAGTTCCACGCGCCGGGTGGGGGGCGGCGGCCGGGCCATTGGCGGGAGGTAGGAGCCGTCACCGt                     >  4:90628/1‑135 (MQ=255)                                                                                                       cGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTGACGGCAGGGGCGGCCATCAGCGCCCCATCCTGCCGTTCCCCCCGCAGGGTGGCGGCCGGCGGCCGGGCAATGCGCTGCCGGTTGAAGACTTAACCGt                     >  6:111806/1‑135 (MQ=255)                                                                                                        ggtGCCGTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTcc                   <  7:180152/134‑1 (MQ=255)                                                                                                          tGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCa                  <  4:218746/135‑1 (MQ=255)                                                                                                            ccTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGGCCATCAGCGACACCTCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCCGGc               >  8:95393/1‑136 (MQ=255)                                                                                                              tttGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCt              <  4:304848/135‑1 (MQ=255)                                                                                                                       gccggccgGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt     <  5:155091/135‑1 (MQ=255)                                                                                                                       gccggccgGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt     <  5:74018/135‑1 (MQ=255)                                                                                                                         cggccGGGCGCGCCCGGGTGACCGGAGGGGCGTCCCTCAGCGAGACATCCAGCACTTCCCCGCTCAGGGGGGCGGCAGGCGGCACGGCAATTAGCTGGAGGTAGAAGACTTCACCGGCCAGGCTGGCCTGGTtggg  >  8:118738/1‑136 (MQ=255)                                                                                                                              |                                                                                                                                   TGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGG  >  ACB122/1718338‑1718593

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