Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,593,393	C→T	11.4%	A39T (GCG→ACG)	CJ019_04342 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,593,393	0	C	T	11.4%	97.5 / 4.4	35	A39T (GCG→ACG)	CJ019_04342	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/20);  new base T (2/2);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

ACCTCAAGTCCTCGGGCGCCCGCCCCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGAATAATCATCTCACCTTAACTCTCCC  >  ACB122/4593261‑4593524                                                                                                                                     |                                                                                                                                    aCCTCAAGTCCTCGGGCGCCCGCCCCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAg                                                                                                                                  <  3:178806/136‑1 (MQ=255)   cgcAAGTCCTCGGGCGCCCGCCCCGTCGCATGGAGCAATACGCGACGCACCGGGGGGGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGgc                                                                                                                                <  5:27091/134‑1 (MQ=255)   cTCAAGTCCTCGGGCGCCCGCCCCGTCTCTTGGAACAATCCGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCGGCGGCGTGAATCCGGGCGTACCCGCGGCCTGATCGAAACACGCCAGgc                                                                                                                                <  3:78631/136‑1 (MQ=255)    tCAAGTCCTCGGGCGCCCGCCCCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCCGCTCCGGCGCTGCCCACTGTCGCTGCGGCGGGAAACCCGGCGTACCCCCCGCCGGATAGCCACACGCCAGgcg                                                                                                                               >  3:175620/1‑136 (MQ=255)        ggccTCGGGCGCCCGCCCCGTCTCATGGATCACTACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCg                                                                                                                           <  1:10660/134‑1 (MQ=255)                cgaaaggcccCGTATCATGGATCAATACGCGACGCTCGGTGGGTGCCTCCAGCTGCTGCGCTGACCACTGTCGGTGCGGGGTGAATCCCGGCGTACCCGCCGCCTGATCTACACACGCCAGGCGCTCGATATTGGa                                                                                                                   <  3:153694/130‑1 (MQ=255)                      gcccCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAACCCCGCCAGGCGCTCGATATTGGACACGag                                                                                                             >  8:40694/1‑135 (MQ=255)                      gcccCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGa                                                                                                                   >  4:72864/1‑130 (MQ=255)                      gcccCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGa                                                                                                                   <  3:72864/130‑1 (MQ=255)                                    aTCAATACGCGACGAACCGTGGGTGCCTCCAGCTCCTTCGCTCACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGAACGAAACACCCCAGGAGCTCGATAATGGACACCATCCGTGGGcctgtct                                                                                               >  4:153694/1‑130 (MQ=255)                                                 gCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGt                                                                                  >  2:3577/1‑136 (MQ=255)                                                         ggTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAaggcag                                                                          <  3:144016/136‑1 (MQ=255)                                                                ccAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTaa                                                                   <  7:57417/136‑1 (MQ=255)                                                                  aGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGTCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAAc                                                                  >  3:87668/1‑135 (MQ=255)                                                                  aGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGTCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAAc                                                                  >  3:90921/1‑135 (MQ=255)                                                                            cTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGAGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTaa                                                       >  6:78718/1‑136 (MQ=255)                                                                             tGACCACTGTCGCTGCGGCGTGAATCCCGGCGTAACCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCTGTGGGCAGTCCGGCGATGCACGGGTCACGGCAGGCATTAACAAGAAAAGTAAt                                                      >  7:41745/1‑136 (MQ=255)                                                                                  aCTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAAAACGCCAGGCGCTCGATATTGGACACGAACCGAGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAg                                                 >  3:193998/1‑136 (MQ=255)                                                                                       cGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGTCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCa                                            <  4:87668/136‑1 (MQ=255)                                                                                       cGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCa                                            <  5:32619/136‑1 (MQ=255)                                                                                        gCTGCGGCGTTAATCCCGGCGTACCCGCCGACTGATCCAACCACGCCAGGCGCTAGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCAc                                           <  7:83545/136‑1 (MQ=255)                                                                                        gCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGTCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCa                                            <  4:90921/135‑1 (MQ=255)                                                                                          tGCGGCGTGAATCCCCGCGTACCCGCCGCCTGATCCAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCAc                                           <  1:198517/134‑1 (MQ=255)                                                                                            cggcgTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGc                                         <  1:67799/134‑1 (MQ=255)                                                                                            cggcgTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGc                                         <  2:118981/134‑1 (MQ=255)                                                                                            cggcgTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCa                                       >  8:119920/1‑136 (MQ=255)                                                                                            cggcgTGAATCCCGGCGTACCCGCCGACTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGGGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATTAAAGTCCCACGCCa                                       <  5:78718/136‑1 (MQ=255)                                                                                                 tGAATCCCGGCGTACCCGCCGCCTGATCGAAACCCGCCAGGCGCTCGATATTGGACCCGCTCCGTGGGCAGTCCGGCGATGCCCGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTg                                   >  1:33212/1‑135 (MQ=255)                                                                                                 ggAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGt                                  <  7:99800/135‑1 (MQ=255)                                                                                                        cGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGa                           <  5:98171/136‑1 (MQ=255)                                                                                                           cGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGaat                         >  4:24795/1‑135 (MQ=255)                                                                                                             tACCCGCCGCCTGATCCAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGaataat                      <  5:124679/136‑1 (MQ=255)                                                                                                               cccGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGAATAatca                    <  1:3577/136‑1 (MQ=255)                                                                                                                               aaCACGCCTGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGAATAATCATCTCACCTTAACtctc    <  7:14232/136‑1 (MQ=255)                                                                                                                                 cacGCCATGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGAATAATCATCTCACCTTAACTCTccc  <  3:60437/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    ACCTCAAGTCCTCGGGCGCCCGCCCCGTCTCATGGATCAATACGCGACGCACCGTGGGTGCCTCCAGCTCCTGCGCTGACCACTGTCGCTGCGGCGTGAATCCCGGCGTACCCGCCGCCTGATCGAAACACGCCAGGCGCTCGATATTGGACACGATCCGTGGGCAGTCCGGCGATGCACGGGTCAAGGCAGGCATTAACAAGAAAAGTAATGAAAGTCCCACGCCACCTGTATAACGAATAATCATCTCACCTTAACTCTCCC  >  ACB122/4593261‑4593524

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