Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,718,463	T→G	36.3%	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	36.3%	32.7 / 10.1	57	Q65P (CAG→CCG)	CJ019_01614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/27);  new base G (14/8);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 1.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGG  >  ACB122/1718335‑1718596                                                                                                                                 |                                                                                                                                      tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTCCCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgc                                                                                                                                  <  3:141394/134‑1 (MQ=255) tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgc                                                                                                                                  <  8:354927/134‑1 (MQ=255) tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgc                                                                                                                                  <  7:42117/134‑1 (MQ=255)    tgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCTGTGGCCAGGCTTGACCTGCGCGGGCTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGCCCGGgcggggcg                                                                                                                              >  5:53902/1‑130 (MQ=255)     gaccAGCAACTTGCCATTCAGCTCAATGCGCGCAATTCCCGCATAGCGGTGGCCAGGCTTGACCGGCGGGGGGGAATCGGTGCGATGCAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAg                                                                                                                             <  8:414941/132‑1 (MQ=255)     gaACAGCAACTTGCCATTCCGCTCAAGGCGCGCACTTCCCGCATCGCGGTGGCCAGGCGTGACCTGCGCGGGGGCATCGGTGAGAGGCACCGGCAGCGGCACGTGGCCTTTGACCGGGCCGGCCGGGCGCGCCAg                                                                                                                             <  4:340632/135‑1 (MQ=255)     gaACAGCAACTTGCCATTCAGCTCAATGCGCGCCCTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATTGGGGAGTTGGAACGCCAGCGGCAGGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                             <  2:150732/135‑1 (MQ=255)     gaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                             <  4:408543/135‑1 (MQ=255)     caaCACCAACTTCCCATTCAGCTCAATGCGCGCACTTACCGCATCGCGGTGGCCCGGCTTGCCCGGCGCGGGGTCATCGGTGAGAGGGACCGGCAGCGGCACGTGGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                             <  4:328661/134‑1 (MQ=255)       aCAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAggg                                                                                                                           <  3:191492/135‑1 (MQ=255)           cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCTCGCGGTCATAGGTGAGATGGAACGGCAGCGGCACGTGGCCATTGACCGGGACGGCCAGGCGCGCCAGGGTCAc                                                                                                                       >  8:449905/1‑135 (MQ=255)           cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGCCCGGGGCAc                                                                                                                       >  5:25983/1‑135 (MQ=255)              cTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                             <  6:100813/126‑1 (MQ=255)              cTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGcccg                                                                                                                             >  5:100813/1‑124 (MQ=255)              cTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTcgaggc                                                                                                                    >  8:198371/1‑130 (MQ=255)               ttGCCATGCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCa                                                                                                                   <  8:243284/135‑1 (MQ=255)                        aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCGGGGGGCGCCCGGGGCCCGGCCGGGGCGGcc                                                                                                          >  3:204678/1‑135 (MQ=255)                        aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTGCCCGCAGGGGCGGcc                                                                                                          >  1:152479/1‑135 (MQ=255)                        aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCGCGGCCGGGGCGGcc                                                                                                          >  3:212373/1‑135 (MQ=255)                              ttGCGCGCACTTACCGAATAGCGGTGGCCAGGCTTGCCCGGCGGGGGGTCATAGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                    <  6:14009/134‑1 (MQ=255)                               tGCGGGCTCTTCCCGAATCGGGGTGGCCTGGGTTGCCCTGGGCGGGGGCATCGGGGAGATGTAACGGCAGCGGCACGTTGCCTTTGACCGGCCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                                   <  6:398313/135‑1 (MQ=255)                                  gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCAGGGGCGGCCATCAGCgcc                                                                                                 >  8:250932/1‑132 (MQ=255)                                  gcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCCCGGCAGGGGCGGCCATCAGCGcca                                                                                                >  4:320145/1‑135 (MQ=255)                                      aCTTACCGCATCGCGGTGGCCAGCCTGGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGTCAGTTTCTTACATCAGCGACACATc                                                                                            <  5:344491/135‑1 (MQ=255)                                      aCTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGACCGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGTGGCGTCCATCAGCGACACATc                                                                                            <  3:224116/135‑1 (MQ=255)                                       cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTGCATCAGCGACACATcc                                                                                           >  4:72717/1‑135 (MQ=255)                                        ttACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCt                                                                                          <  2:59382/135‑1 (MQ=255)                                                   cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                                <  6:362281/134‑1 (MQ=255)                                                         ccAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACCGCAGGGGCGTCCCTCACCCGCACCTCCGCCACTTCCCCGCACAgg                                                                         >  6:167588/1‑135 (MQ=255)                                                               ttGACCTGCGCGGGGCCATAGGTGAGATGGAACGGCAGCGGCACGTAGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                   <  7:435782/135‑1 (MQ=255)                                                               tggACCTGCGCGGGGTCATAGGGGCGAGGGAACGGCCCCGGCACGTTGCCTTTGCCCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                   <  7:404664/133‑1 (MQ=255)                                                                tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTGACGGCAGGGGCGTGCCTCAGCGCCCCCACCTGCAGTTCCCCGCTCAGGGTggcgg                                                                  >  3:328897/1‑135 (MQ=255)                                                                    cTGCGCGGGGTCATCGGGGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                              <  2:1395/135‑1 (MQ=255)                                                                    cTGCGCGGGGTCATAGGTGAGAGGCAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                              <  1:45856/135‑1 (MQ=255)                                                                     tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGCGCGGCAGGGGGGGGCCGCCGCGCCACAAACAGCAGGTGCACGCGCACGGGggcgggcggc                                                             >  5:359794/1‑130 (MQ=255)                                                                          ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGGCCCGGCCGGGGCGGCCCGCAGCGACACATCCAGCCGTTCCACGCGCACGGTggcggcaggcggccg                                                        >  1:308084/1‑133 (MQ=255)                                                                          ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCCCGGCCGGGGCGTCCATCAGCGACACATCCTGCAGTTGCACGCGCAGGGTGgcggcaggcggcag                                                        >  2:43023/1‑135 (MQ=255)                                                                               cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                                <  3:116889/106‑1 (MQ=255)                                                                               cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                                >  4:116889/1‑106 (MQ=255)                                                                                 tAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCAGGCGCGCCCGGGTCCCGGCCGGGGGGGCCATCAGCGCCACCTCCAGCAGTTCCACGCGCAGGGTGGCGGCACGCGGGAGGGCCATg                                                 >  4:48219/1‑135 (MQ=255)                                                                                  aGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgc                                                <  6:401671/135‑1 (MQ=255)                                                                                  aGGTGAGAGGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgc                                                <  6:228628/135‑1 (MQ=255)                                                                                     ggagaTGGACCGGCCGCGGCACGTGGCCGTGGCCCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                              <  3:48219/133‑1 (MQ=255)                                                                                     ggagaTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                              <  2:246188/133‑1 (MQ=255)                                                                                          tGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACAACCTGCAGTTCCACCCGCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAgg                                        >  2:378136/1‑135 (MQ=255)                                                                                                 gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTagaaga                                 <  4:330672/135‑1 (MQ=255)                                                                                                 gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCCCGGCCGGGGCGGGCCTCAGCGCCCCATCCTGCAGTTGCCCGCGCAGGGTGGGGGGCGGCGGCCGGGCCATGCGCTGCAGGGagaaga                                 >  2:252745/1‑135 (MQ=255)                                                                                                   cgCGGCTCGTTGCCTTTGCCCGGCGCGGCGTGGCGCGCCAGGGCCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                               <  1:200195/134‑1 (MQ=255)                                                                                                   aGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                               <  6:156182/135‑1 (MQ=255)                                                                                                            tggCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTcc                      <  5:121340/133‑1 (MQ=255)                                                                                                               ccTTTGACCGGGCCGGCCGGGCGCGCCCCGGGCCCGGCAGGGGCGGGCCTCAGCCCCCCCTCCAGCAGTTCCACGCACACGGGGGCGGCAGGCGGCCGGGCAATGCGCGGCCGGGAGAAGACTTAACCGTcccgc                   >  8:31413/1‑132 (MQ=255)                                                                                                                cTTTGCCCGGGCCGGCCGGGCGCGCCAGGGCCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGc                  <  5:363382/135‑1 (MQ=255)                                                                                                                  tggcccGGGCGGGCCGGGCGCGCCTGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCGCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTg                <  1:188052/131‑1 (MQ=255)                                                                                                                           ggggggtGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGTCAGGGCAATGCGCTGCAGTGAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg       <  3:323270/129‑1 (MQ=255)                                                                                                                           cgggctTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg       <  1:256655/133‑1 (MQ=255)                                                                                                                            cggccTGGCGCGCCAGGGTCACGGCAGGGGCGTGCCTAAGCGCCACAACCAGCAGTTCCACGCGCAGGGGGGCGGCAGGCGGCAGGGCAATGAGCTGCAGGTAGAAGACTTAACCGTCCAGGCGTGGCTGGGtgg      >  6:107955/1‑135 (MQ=255)                                                                                                                               ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGGTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtg   <  2:152479/135‑1 (MQ=255)                                                                                                                                cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtgg  <  3:176911/135‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGG  >  ACB122/1718335‑1718596

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