Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,912,296	A→T	15.6%	S244C (AGC→TGC)	CJ019_04608 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,912,296	0	A	T	15.6%	50.7 / 3.4	25	S244C (AGC→TGC)	CJ019_04608	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base T (2/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCACCGCCAGCGTCGCCACCG  >  ACB122/4912171‑4912413                                                                                                                              |                                                                                                                      aaCAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGa                                                                                                             >  7:146273/1‑136 (MQ=255)   cAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACg                                                                                                           <  3:148391/136‑1 (MQ=255)      ctgctgACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCCCCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGccc                                                                                                        >  6:13222/1‑136 (MQ=255)      ctgctgACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGCCCTGCGCGCCAGCGCCGCACGCCTGCTGGCAAGCCTGCCGGACGccc                                                                                                        >  2:204888/1‑136 (MQ=255)                aGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTggg                                                                                              <  3:99174/136‑1 (MQ=255)                   aaaaCCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGa                                                                                           >  8:6973/1‑136 (MQ=255)                          ccAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGAGGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAg                                                                                    <  7:11604/136‑1 (MQ=255)                                    cTACCCGGCCTGCTCACCGCCCTTCAACAGCTCACCCGCTACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGccc                                                                          <  7:47340/136‑1 (MQ=255)                                                             aaCAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCACGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCtt                                                 >  1:192322/1‑136 (MQ=255)                                                             aaCAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCtt                                                 >  6:87806/1‑136 (MQ=255)                                                                    cccccGCGACCCAGGCGGCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTcgaa                                          <  5:169526/134‑1 (MQ=255)                                                                    cACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTcgaa                                          <  7:6973/136‑1 (MQ=255)                                                                      cccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTcgaagc                                        >  1:198041/1‑136 (MQ=255)                                                                       ccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCATGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGaa                                                                <  7:152692/111‑1 (MQ=255)                                                                       ccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCATGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGaa                                                                <  7:155363/111‑1 (MQ=255)                                                                       ccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCATGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGaa                                                                >  8:155363/1‑111 (MQ=255)                                                                       ccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCATGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGaa                                                                >  8:152692/1‑111 (MQ=255)                                                                       ccGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTcgaagcg                                       <  3:196027/136‑1 (MQ=255)                                                                                      cGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTggg                        <  8:146273/136‑1 (MQ=255)                                                                                          ggCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGc                                                      >  7:31763/1‑102 (MQ=255)                                                                                          ggCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGc                                                      <  8:31763/102‑1 (MQ=255)                                                                                          ggCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGAAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCa                    >  2:74798/1‑136 (MQ=255)                                                                                                      gcCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCACCGCCAGCGTCg        >  8:122882/1‑136 (MQ=255)                                                                                                        cAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCACCGCCAGCGTCGcc      <  2:198041/136‑1 (MQ=255)                                                                                                            gcCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCACCGCCAGCGTCGCCACCg  >  7:72053/1‑136 (MQ=255)                                                                                                                              |                                                                                                                      AACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAACGCCTGCTGGCAAGCCTGCCGGACGCCCGACAACTGGGCGACGCCAAGGCAGTGGCCCAGGCCCTGAACAACAGCGGCAGCTTTCTCGAAGCGAAGCTGCTGAGTGGGTCCACCGCCAGCGTCGCCACCG  >  ACB122/4912171‑4912413

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