Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	339,836	T→G	36.1%	S354R (AGC→CGC)	CJ019_00314 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	339,836	0	T	G	36.1%	28.7 / 2.1	28	S354R (AGC→CGC)	CJ019_00314	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/7);  new base G (2/9);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 5.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTCGGCGTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAGGCCGCCGCT  >  ACB122/339706‑339968                                                                                                                                   |                                                                                                                                     cgtgcGGCGTCGAAGGTGGCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCTAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCGGGgt                                                                                                                                  <  2:215755/131‑1 (MQ=255) cGGTCGGCGTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGgt                                                                                                                                  <  4:241320/135‑1 (MQ=255) cGGGCGGCGGCGAAGGGGGCGCCCTTAAGGTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGgt                                                                                                                                  <  8:297368/135‑1 (MQ=255)  ggTCGGCGTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCCGGTCTTCCAGGGCCAGGCTGCTGGgtg                                                                                                                                 >  6:301555/1‑135 (MQ=255)      ggCGTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCGGCAGGTGCTGTCGTGCCGCTTGACCGCCCGGGGGTGCCGGGGCCGGCGGCGGGGtgtggc                                                                                                                             >  1:95748/1‑133 (MQ=255)         gTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCCCGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTGCAGGGCCCGGCTGCTGGGTGTGCCcag                                                                                                                          >  3:166803/1‑135 (MQ=255)           cGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCccgca                                                                                                                        >  7:212363/1‑135 (MQ=255)           cGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCcagca                                                                                                                        >  7:206909/1‑135 (MQ=255)                 ggcggaccTTATGTTCCAGGCGCAGGGCTTGTTGGGCTAAATCTTCGACGGCCACGCGGCCGCTGATGGGGCTGTCGTCCAGCTACACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAg                                                                                                                  <  4:239541/129‑1 (MQ=255)                       ccTTAAGTGGCTGGCGCTGGGCTTGTTTGGCTAAAGCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGGCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCAcca                                                                                                            <  8:184778/135‑1 (MQ=255)                                    gcgcAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTgg                                                                                                >  6:175758/1‑134 (MQ=255)                                    gcgcAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGc                                                                                               >  7:79906/1‑135 (MQ=255)                                                       aaaTCTTCGACGGCCACGCGGCCGCTGAACGTGCTGTCGTCCAGCTTCACACCCAGGTCTTACAGCGCCCGGCATCTGGGTGTTCCCAGGAGGAGAGGCACCCGCTCCAGCTTGGCGAAGGCGGCCGGGGCGTtg                                                                            >  5:301067/1‑135 (MQ=255)                                                                  ggCCACGCGGCCGCGGACGGTGCTGGCGTCCAGCTGCACCGCCAGGGCTTCCAGGGCCAGGCTGCTGGGGGGGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAg                                                                 <  7:241141/135‑1 (MQ=255)                                                                          ggCCGGTCAAGGTGGTGGCGACCCGCTACCCCCCCTGGTCTTCCCGGGCCAGGGTGCGGGGTGTCCCCCGCTGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGc                                                         <  6:10412/135‑1 (MQ=255)                                                                              gcggAAGGGGCGGCCGTCCAGCTGCACCGCCAGGACTTCCAGGGCCCGGCTGCGGGGTGGGCCCCGACGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAAt                                                     <  2:80500/132‑1 (MQ=255)                                                                                      tGCTGGCGGCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGGTGGGGGTGCCCAGCTGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCaa                                             <  2:211403/135‑1 (MQ=255)                                                                                          gtcgtcCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCGGGGTGTGCCCAGCAGGCGAGTCACCCGCTCCAGCTTGGCGAATGCGGCCGGGTCGTTGGTGGCCCGCAGCGGGTGGCCAATGCGTTGaaggaa                                         >  6:229071/1‑135 (MQ=255)                                                                                            ggccccgCTTCACCGCCAGGGCTTCCGGGGCCCGGGGGGGGGGGGTCCCCCGCCGGCGAGGCCCCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCaaggaagg                                       <  7:225527/129‑1 (MQ=255)                                                                                            cgacccgCTTCACCGCCAGGGCTTCCAGGGCCAGGCGGCGGGGGGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCaaggaagg                                       <  6:328895/129‑1 (MQ=255)                                                                                            cgacacgCTGCACCGCCAGGTCTCCCCGGGCCCGGCGGCTGGGTGGGCCCAGCAGGCGAGTCACCAGCTCCAGGTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCaaggaagg                                       <  8:72284/129‑1 (MQ=255)                                                                                              tcCAGCTTCACCGCCTGGTCTTCCAGGGCCAGGCGGCGGGGGGGGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGgcg                                     <  4:306296/135‑1 (MQ=255)                                                                                                   cTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAg                                >  3:193162/1‑135 (MQ=255)                                                                                                        ccGCCAGGTCTTCCAGGGCCAGGCGTCGGGGTGTCCGCTGCAGGCTAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGGGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGa                           <  8:13850/135‑1 (MQ=255)                                                                                                           ccAGGTCTTCAAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGCAGGAGCGCAGGTTGAACt                        >  3:307033/1‑135 (MQ=255)                                                                                                                       gggCCCGGCGGCGGGGGGGCCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACa            <  3:227418/135‑1 (MQ=255)                                                                                                                         gCCCGGCTGGGGGGTGGGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAgg          <  5:175758/135‑1 (MQ=255)                                                                                                                             ggCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAGgccgc      >  6:131388/1‑135 (MQ=255)                                                                                                                               ctgctgGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAGGccgccg    >  8:288115/1‑135 (MQ=255)                                                                                                                                 gctgGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAGGCCGCCGCt  >  4:374305/1‑135 (MQ=255)                                                                                                                                   |                                                                                                                                     CGGTCGGCGTCGAAGGTGTCGCCCTTAAGTTGCAGGCGCAGGGCTTGTTTGGCGAAATCTTCGACGGCCACGCGGCCGCTGAAGGTGCTGTCGTCCAGCTTCACCGCCAGGTCTTCCAGGGCCAGGCTGCTGGGTGTGCCCAGCAGGCGAGTCACCAGCTCCAGCTTGGCGAAGGCGGCCGGGTCGTTGGTGGCCGGCAGCGGGTGGCCAATGCTTTCAAGGAAGGCGCGCAGGTTGAACTGGGCGATGGACAGGCCGCCGCT  >  ACB122/339706‑339968

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