Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	3,660,726	T→G	19.4%	D334A (GAC→GCC)	CJ019_03469 ←	putative protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	3,660,726	0	T	G	19.4%	54.5 / 5.1	38	D334A (GAC→GCC)	CJ019_03469	putative protein
Reads supporting (aligned to +/- strand):  ref base T (2/28);  new base G (5/3);  total (7/31)
Fisher's exact test for biased strand distribution p-value = 2.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGTTGTGAGCCCAGACCACCG  >  ACB143/3660597‑3660861                                                                                                                                  |                                                                                                                                        cAGTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCCCCTGGCCGCGGGGGGTgg                                                                                                                                   >  7:30130/1‑136 (MQ=255) cAGTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTgg                                                                                                                                   <  5:335012/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  8:213928/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  1:4766/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  8:155984/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  5:62137/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  5:168466/136‑1 (MQ=255)   gTCGCTCAGGCAGGCTGGCACGCCGGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGcc                                                                                                                                 <  4:330818/136‑1 (MQ=255)          agggcggcTGGCACGCCTGCGAGGCGGAACTGATGCGCCCCGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGGCGCTGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGGCGGTGGCCATGcc                                                                                                                            <  1:364391/129‑1 (MQ=255)          aggcaggcTGGCACGCCGGCGAGGCGGAACTGATGCTCCCAGCTGGCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGcc                                                                                                                            <  2:381042/134‑1 (MQ=255)                       cGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGGCCCCAGGCCCCGGGCGGGGGGCAAGCGGAGCAACACcgccgg                                                                                                             >  3:354133/1‑133 (MQ=255)                               aGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGa                                                                                                     <  3:48537/136‑1 (MQ=255)                                  cGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGGCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAg                                                                                                  <  4:389550/136‑1 (MQ=255)                                             cTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCa                                                                                        <  1:257987/135‑1 (MQ=255)                                             cTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGGGGGTGGGCAATCCGATCCACCCCCCCGGGTTCCCCAAAACGCggcggcc                                                                                        >  2:43896/1‑134 (MQ=255)                                                 cAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCt                                                                                   <  6:210155/136‑1 (MQ=255)                                                 cAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCt                                                                                   <  5:203917/136‑1 (MQ=255)                                                           tGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGtt                                                                         <  6:21079/136‑1 (MQ=255)                                                             cgcgAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGtt                                                                         <  3:6882/134‑1 (MQ=255)                                                             cgcgAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGtt                                                                         <  4:51884/134‑1 (MQ=255)                                                             cgcgAAGGTCGGATTTCCTTGAGCTGCATGGCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGtt                                                                         <  1:205118/134‑1 (MQ=255)                                                                gAAGGTCGGATTTCCTTGATCTGCATTTCGCCGTCCCAGTCATCGGCGGCGGCCACCTTGCCCCGGTCGGTGGCCCAGGCGATAAACACCCCGTCGGGCCCGAAAGCGGAGACGCAAAACAGGGCGCTGTTGaaat                                                                    >  8:99277/1‑134 (MQ=255)                                                                gAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACt                                                                    <  8:236321/136‑1 (MQ=255)                                                                gAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACt                                                                    <  1:53432/136‑1 (MQ=255)                                                                 aaGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCATTCATCGGCGGCGGCCAACTGGCCGCGGTGTGTGGCCATTGCGGACAAAAACCCGTGGGGGCCGCAAGCGCTGAGGCACAGCGGGCCGATGTTTaaacg                                                                   >  2:388578/1‑133 (MQ=255)                                                                     tCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCt                                                               <  2:236649/136‑1 (MQ=255)                                                                                       cATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCAttt                                                  <  6:386426/131‑1 (MQ=255)                                                                                        aTGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTgg                                            <  7:125342/136‑1 (MQ=255)                                                                                        aTGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTgg                                            <  6:382702/136‑1 (MQ=255)                                                                                          gTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGGAGGTGGCCATGCCGATCAACACCGCGTCGGGCCCGAAAGCGCTGCGGGAAAGCTGGGCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCt                                          >  7:256684/1‑136 (MQ=255)                                                                                             gCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGGCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCg                                       <  3:271817/136‑1 (MQ=255)                                                                                               cGTCCCAGTCATCGGCGGCGGCCACCTGGGCGCGGGCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGc                                     <  4:189669/136‑1 (MQ=255)                                                                                                      gTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCg                              <  2:387426/136‑1 (MQ=255)                                                                                                              cggcggcCACCTGGCCGCGGGGGGTGGGCCTGCCGCTCAACACCCCGGGGTGCCCGAAAGCGCGGCGGCAAAGCTGGCCGATGTTTTACAGGCCATTACCGCCCAATTGCCTGGCTGCGGCCTTGGCGATATgggg                      >  2:338423/1‑135 (MQ=255)                                                                                                                gcggcCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGt                    <  7:197237/136‑1 (MQ=255)                                                                                                                gcggcCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGt                    <  5:179544/136‑1 (MQ=255)                                                                                                                gcggcCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTAGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTTGCCTTTCCAGCACATTTCCGTGGCTGCGGCATTGCCGATATGGGAGt                    <  2:32921/136‑1 (MQ=255)                                                                                                                            ccGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGTTGTGAGCCCAGa        <  1:53745/136‑1 (MQ=255)                                                                                                                                  tcggtGGCCATGCCGATCAACACCGCGTCGTGCCGGAAAGCGCTGCGGCAATGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGTTGTGAGCCCAGACCACCg  <  2:130565/136‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        CAGTCGCTCAGGCAGGCTGGCACGCCTGCGAGGCGGAACTGATGCTCCCAGCTGTCGCTGCGCGAAGGTCGGATTTCCTTGATCTGCATGTCGCCGTCCCAGTCATCGGCGGCGGCCACCTGGCCGCGGTCGGTGGCCATGCCGATCAACACCGCGTCGTGCCCGAAAGCGCTGCGGCAAAGCTGGCCGATGTTGAACTGGCCTTTCCAGCCCATTTCCGTGGCTGCGGCATTGCCGATATGGGAGTTGTGAGCCCAGACCACCG  >  ACB143/3660597‑3660861

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