Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	5,994,584	G→A	20.2%	G220G (GGC→GGT)	CJ019_05610 ←	Arylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,994,584	0	G	A	20.2%	125.6 / 23.7	56	G220G (GGC→GGT)	CJ019_05610	Arylesterase
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (6/5);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

CAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGGTCTGCGTCTGCACA  >  ACB143/5994458‑5994719                                                                                                                               |                                                                                                                                        cAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGAt                                                                                                                                 >  7:163911/1‑135 (MQ=255)      gctgcGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCaa                                                                                                                           >  5:49546/1‑136 (MQ=255)       ctgcGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCaa                                                                                                                           <  2:19461/135‑1 (MQ=255)       ctgcGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAAc                                                                                                                          <  4:83750/136‑1 (MQ=255)       ctgcGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAAc                                                                                                                          >  5:233999/1‑136 (MQ=255)       ctgcGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTCCTTGCCAGTGGTCTCGAACGGCACACTCTGGTCGTCATCGCCATGGATCACCAAc                                                                                                                          >  3:229944/1‑136 (MQ=255)               gAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCCCGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCAGCCCCATGGATCACCAACGTCGGcac                                                                                                                  >  8:261012/1‑136 (MQ=255)                    aCTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGa                                                                                                             >  3:225182/1‑136 (MQ=255)                       gCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGCCCt                                                                                                          >  2:207419/1‑136 (MQ=255)                       gCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCt                                                                                                          >  6:42132/1‑136 (MQ=255)                       gCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCt                                                                                                          >  7:4721/1‑136 (MQ=255)                       gCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCt                                                                                                          >  7:87340/1‑136 (MQ=255)                         aaaaCCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTg                                                                                                        >  2:253122/1‑136 (MQ=255)                         aaaaCCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTg                                                                                                        >  7:35914/1‑136 (MQ=255)                              cGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCAt                                                                                                   <  5:129111/136‑1 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  5:201010/1‑136 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  5:110108/1‑136 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  4:160561/1‑136 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  7:167457/1‑136 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  7:164056/1‑136 (MQ=255)                                   ggcgCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTAATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcg                                                                                              >  8:199851/1‑136 (MQ=255)                                          cGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcggacgga                                                                                        >  7:29666/1‑135 (MQ=255)                                          cGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCcggacgga                                                                                        >  2:84672/1‑135 (MQ=255)                                           ggCATACCCCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCCACGTCGGCACATCGACCTTGGCCATGTCCGGACGGaaa                                                                                      >  3:268041/1‑136 (MQ=255)                                              aTACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCg                                                                                   <  8:222675/136‑1 (MQ=255)                                                    cTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGGCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTcc                                                                             <  5:184441/136‑1 (MQ=255)                                                     ttCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTcc                                                                             >  5:143065/1‑135 (MQ=255)                                                       cAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCga                                                                           <  8:4721/135‑1 (MQ=255)                                                         gCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCgag                                                                          >  4:286993/1‑134 (MQ=255)                                                         gCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCgag                                                                          >  1:87432/1‑134 (MQ=255)                                                          cTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAg                                                                       >  1:7103/1‑136 (MQ=255)                                                             ggCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAgg                                                                      <  4:68310/134‑1 (MQ=255)                                                               cgcCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGt                                                                  <  6:143065/136‑1 (MQ=255)                                                                gcCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTc                                                                 <  3:64014/136‑1 (MQ=255)                                                                 ccACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTcc                                                                >  5:253602/1‑135 (MQ=255)                                                                          aGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATcc                                                        <  7:193983/135‑1 (MQ=255)                                                                          aGCTCCGCCACCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCa                                                       <  8:243958/136‑1 (MQ=255)                                                                               cgccgcCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAggg                                                    <  4:267178/134‑1 (MQ=255)                                                                                   gccTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCtt                                              <  7:105007/136‑1 (MQ=255)                                                                                    ccTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCttt                                             >  1:260193/1‑136 (MQ=255)                                                                                           gCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATgg                                        <  3:39877/134‑1 (MQ=255)                                                                                              aGTGGTCTCGAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGc                                    >  3:284363/1‑135 (MQ=255)                                                                                                       gAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAgcg                            <  7:172741/134‑1 (MQ=255)                                                                                                       gAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAgc                             <  5:181877/133‑1 (MQ=255)                                                                                                       gAACGGCACAATCTGGTCGTCATCACCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAgc                             <  5:237756/133‑1 (MQ=255)                                                                                                       gAACGGCACAATCTGGGCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAgcg                            <  1:190929/134‑1 (MQ=255)                                                                                                           ggCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGt                      >  4:32442/1‑136 (MQ=255)                                                                                                                aaTCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAgg                 <  8:20745/136‑1 (MQ=255)                                                                                                                aaTCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAgg                 >  3:134581/1‑136 (MQ=255)                                                                                                                        cgtcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgt         >  2:86894/1‑136 (MQ=255)                                                                                                                         gtcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCGTTGTCCGGACGGAAATATGTCTACGAGAAGGCCGTCACACAATCCAGTGTGACTTTGATGGACGCCATCAGCGAAATGTTCATGttcttctt         <  8:46177/135‑5 (MQ=255)                                                                                                                          tcatcaCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  1:84672/136‑1 (MQ=255)                                                                                                                          tcatcaCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  5:62531/136‑1 (MQ=255)                                                                                                                          tcatcaCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  4:284363/136‑1 (MQ=255)                                                                                                                          tcatcaCCATGGATCACCAACGTCGGCACATAGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAAGGTTCAGGgtcggcgtct       >  4:169196/1‑136 (MQ=255)                                                                                                                          tcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  1:207419/136‑1 (MQ=255)                                                                                                                          tcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  5:271460/136‑1 (MQ=255)                                                                                                                          tcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  5:42132/136‑1 (MQ=255)                                                                                                                          tcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  7:112158/136‑1 (MQ=255)                                                                                                                          tcATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGgtctgcgtct       <  6:268167/136‑1 (MQ=255)                                                                                                                               gCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGGTCTGCGTCTGcaca  >  4:257664/1‑136 (MQ=255)                                                                                                                               gCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGGTCTGCGTCTGcaca  >  1:221178/1‑136 (MQ=255)                                                                                                                               gCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGGTCTGCGTCTGcaca  >  6:84991/1‑136 (MQ=255)                                                                                                                               |                                                                                                                                        CAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCCGCCGCCTGCTTGCCAGTGGTCTCGAACGGCACAATCTGGTCGTCATCGCCATGGATCACCAACGTCGGCACATCGACCTTGGCCATGTCCGGACGGAAATCGGTCTCCGAGAAGGCCGTCACACAATCCAGGGTGCCTTTGATGGACGCCATCAGCGCAATGTTCAGGGTCTGCGTCTGCACA  >  ACB143/5994458‑5994719

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