Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,686,855	C→G	35.1%	L90V (CTG→GTG)	CJ019_01584 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,686,855	0	C	G	35.1%	26.1 / 6.3	33	L90V (CTG→GTG)	CJ019_01584	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/5);  new base G (2/7);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGC  >  ACB143/1686723‑1686969                                                                                                                                      |                                                                                                                   aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCGTTGGCCC‑GCCGGCGGTGCAACGCACCCCCCCCCGCcaggtgcggg                                                                                                                 >  4:209492/1‑136 (MQ=255) aaGCCGATCGCCCAGGCACCTGCCTGCGTCGCGGGCATTCTCTCCCATCGTGGCCAGATGGTTCCGGTCATCGCACTTTCCGCCTTGAGCATTTGGCA‑GCCGGCGGGGCCACGCACCAGCCCCCGCctggtgctgg                                                                                                                 >  8:289336/1‑136 (MQ=255)          gCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCCGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCCCCCGCCTGGTGCTGGTGCACTAcc                                                                                                        >  6:205153/1‑136 (MQ=255)           cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGGCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCCCCCGCCTGGTGCTGGGGCACTACCg                                                                                                       >  4:392748/1‑136 (MQ=255)                         tgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGGTGGTGCACTACCGAGGCGGCCTGCAGc                                                                                         >  5:119330/1‑136 (MQ=255)                         tgGGTCGCGGGCATTCTCGCCCAACGTGGCCATCTGGTTCCGGTCATCGACCTTTCCGCCTTGATCTTTGGCCA‑GCCCGCGGTGCAAACCACCAGCCCCCGCCTGGTGCTGGTTCACTACAGAGGCGGCCTGCAg                                                                                          >  8:87038/1‑135 (MQ=255)                             tCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCCCCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTgg                                                                                     >  2:371966/1‑136 (MQ=255)                             gggggggCATTCCCCCCCTTGGGGGCCTGGTGGGGCCGGTCAGCCACCTGGCCGCGTTGTGCTTGGGCCC‑GCGGGGGGTCCACCCCACCACCCCCCGCGGGGTGGTGGGCCACTACCGAGGCGGCCTGCAGCTTgg                                                                                     <  2:178586/132‑1 (MQ=255)                                   gCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCATTGGCCA‑GCCGGCGGTGCAACCCCCCAGCAACCCCCATGTGCTGGTGCACTACCGACGCGGCCTGCAGCATGGCCTGAt                                                                               >  8:245989/1‑136 (MQ=255)                                   gCATTCTCGCCCATAGTGGCCAGCTGGTTCCGGTCCTCGCCCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCCCCCCCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGGCTGa                                                                                >  8:382904/1‑135 (MQ=255)                                            cccATCGTGGCCCTCGGTTTCGGGGAATCGACCTTCCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCa                                                                      <  1:272263/136‑1 (MQ=255)                                              cATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGGTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAgg                                                                    >  1:163226/1‑136 (MQ=255)                                              cATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAgg                                                                    >  8:56754/1‑136 (MQ=255)                                                    ggCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCAcc                                                               >  7:6261/1‑135 (MQ=255)                                                    ggCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCg                                                              >  2:88121/1‑136 (MQ=255)                                                         gCTGGTTCCGGTCATCGACCTTGCCGCCTTTTGCTTTGGCCA‑GCGGGCGGTGCAACGCACCAGCACCCGCCTGGTGGTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAAcc                                                         <  6:396356/136‑1 (MQ=255)                                                                   gTCATCGACGTTGCCGCATTGAGCTTTGGCCC‑GCCGGGGGTGCACCGCACCAGCACCCGCGTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGcc                                               <  7:102116/136‑1 (MQ=255)                                                                        cGACATTGCCGCCTTGAGCTTTGCCCA‑GCCGGCGGTGAAACGCACCAGCACCGGCCTTGTGATGTTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCg                                          <  8:183693/136‑1 (MQ=255)                                                                           ccTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGAc                                       >  4:137355/1‑136 (MQ=255)                                                                           ccTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCCGTGCAACGCACCAGCACCCGCCTGGTGCTGGAGCCCTACCGAGGCGGCCTGCCGCTTGGCCGGGTCCTGGAGCAGGCCACCCAAACCCTGCGTTGCCAACCGGAc                                       >  5:299330/1‑136 (MQ=255)                                                                              tGCCGCCCTGCGTTTGGCCCC‑GCCGGTGGTCCATCCCACAAGAACCGGCGTGGTGGTGGTGCACTACCGAGGCGGCCTGCAGCTTGTCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACg                                      <  6:233732/134‑1 (MQ=255)                                                                               gcggccTTGAGCTGTGCCAC‑GCGGGGGGTCAACCCCACCCGCACCGGCGGGGGGGTGGGCCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAg                                    <  1:139841/132‑1 (MQ=255)                                                                               gcggccTGGCGCTTTGGCCA‑GCGGGGGGGGCAACGCACCAGCACCCGCGTGGTGGTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAg                                    <  2:127396/132‑1 (MQ=255)                                                                               gccgccTTGTGGTTGGGCCC‑GCCGGGGGTCCAACCCACCCCCCCCCGCGTGGGGGGGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAg                                    <  6:13057/135‑1 (MQ=255)                                                                                ccgctTTGCGGTGTGGCCC‑GCCGGCGGTGAAACGCACCAGCACCCGCCGGGGGTTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                  <  2:8053/136‑1 (MQ=255)                                                                                ccgccTTGAGCTTTGGCCC‑GCCGGGGGTGCAACGCACCACCACCCGCGTGGTGCGGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                  <  3:404917/136‑1 (MQ=255)                                                                                ccgccTGGTGCTTGGGCCA‑GCCGGGGGTGCAACGCACCAGCACCCGCTGGGGGGTGGGGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                  <  5:205153/136‑1 (MQ=255)                                                                                          cTTTGGACC‑GCGGGCGGGCCACCCCCCCAGCCCCCGCGTGGGGTGGGTGCACTACCGAGGCGGCATGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTAt                        <  8:371791/136‑1 (MQ=255)                                                                                           cttGGCCACGCGGGGGTTGAAACCCACCACCACCCGCGTGGTGCTGGTGCACTACCGAGGCGGCATGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTAt                        <  8:371795/135‑1 (MQ=255)                                                                                               gCCA‑GCGGGGGGTGCACCCCACAAGCACCCCCCGGGTGTGGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGt                   <  3:127765/136‑1 (MQ=255)                                                                                                ccA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTg                  >  8:356373/1‑136 (MQ=255)                                                                                                ccA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGAAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTg                  >  5:226379/1‑136 (MQ=255)                                                                                                             gCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGAt      >  2:122937/1‑136 (MQ=255)                                                                                                               accccACCAGCCCCCGCGTGGTGCGGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATgc    <  1:263237/132‑1 (MQ=255)                                                                                                                 cGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATgcgc  >  3:178808/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                   AAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCA‑GCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGC  >  ACB143/1686723‑1686969

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