Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB122	3,440,343	A→C	S18R (AGC→CGC)	CJ019_03266 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,440,343	0	A	C	79.2%	20.9 / ‑3.6	24	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (17/2);  total (21/3)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440217‑3440473                                                                                                                               |                                                                                                                                   cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcg                                                                                                                           >  7:19956/1‑136 (MQ=255)           aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCCCCCggg                                                                                                                  >  1:190875/1‑135 (MQ=255)            aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGc                                                                                                                 >  1:68592/1‑135 (MQ=255)              ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                              >  3:32304/1‑136 (MQ=255)                     gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCCGCCgggggg                                                                                                       >  6:6135/1‑134 (MQ=255)                                        cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCAGCCTCGGGGTGGCGGGCCCACCCCta                                                                                      >  1:92584/1‑134 (MQ=255)                                         aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCCCGCCgat                                                                                     >  2:59441/1‑134 (MQ=255)                                             cccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCCGTGGCGGCGCCCGCGCGGGCGGGGGGTGGCGGGCCCACGCCTCTAATgcg                                                                               >  7:95979/1‑136 (MQ=255)                                                            cgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCGCGCGGCTCa                                                                 >  1:264429/1‑135 (MQ=255)                                                             gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCCCGGGGGGCGGCGCGGCCGCCGGGGGGTGGGCGGGCCCCCCCGCTAAAGCGCCGCTCGCTGATCAg                                                                >  2:327505/1‑135 (MQ=255)                                                             gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGGCCGGGCCGCCGGGGGGGGGGCGGCCCCCCCCCATAACGCGCCGCACGCGGATAAgg                                                               >  6:40642/1‑136 (MQ=255)                                                               tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGGAGCCCGGGGGTGGCCGGCCCACGGCGATAAGTCGCCGCGCGCTGCTCAgggc                                                             >  8:324138/1‑134 (MQ=255)                                                               tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCGGGGCGGTGGCAGGGCCCCCCCCATAAAGCGCCCCTCGGGGATCAGGcc                                                             >  3:301862/1‑136 (MQ=255)                                                               tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCCACGCCGATAAGGCGCCGCGCGCGGCTCAGGcc                                                             >  5:154162/1‑136 (MQ=255)                                                                          cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCCCCCCGATAAGGCGCAGCTCGCGGCTCAGGCCGCTGACTGCgg                                                  >  7:172262/1‑136 (MQ=255)                                                                          cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCTGGCGGCGCCCGGGCCGCCCGGGGGTGGCCGGCCAACGCCCAAAAAGCGCCCCTCGCGGCTCAGGCCGCTGACttag                                                   >  2:23512/1‑133 (MQ=255)                                                                                  gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTTCTGATGGCCCCCACCGCGGGCGGCGCGCGGGCCGCCTGGGTGTGGGCGGCCAACCCCACTAACGCGCCGCTCGGTGCTGAGGCCGGGGGCTTGGGCGAGCAgg                                          >  4:299170/1‑135 (MQ=255)                                                                                  gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCCGGCCAACCCCTATAAGGCGCAGCTCGCGGATCAGGCCGATGACTTGGGCACGCAgg                                          >  7:140461/1‑135 (MQ=255)                                                                                  gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCCCCGGGGTGTGGGCGGCCCACGCCGATAAGGCGCAGCGCGCGGAGCAGGCCGCTGCCTGCGGCAAGCAGt                                          >  4:168894/1‑136 (MQ=255)                                                                                       ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCCCCGCCGAGAAGGCGCCGCGCGCTGATCCGGCCGCTGGCGTCGGCAAGCAGGTGGCg                                     >  7:77227/1‑136 (MQ=255)                                                                                               gCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCGGGCGGGGCCCGGGCGGCCGGGGGGTGGGGGGGCCACGCCCATAATGCGGCGCGCGCTGAGCAGGCCGCGGCCGTAGGCCAGCAGGTGGGGGGGGAggg                             >  4:278881/1‑135 (MQ=255)                                                                                                     gggCCGTGGCGCTGCTGCGGCCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTTGCAGGCCATCGCCTATAAGTCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCTGAGAAGTATcgcgc                       <  8:5629/134‑1 (MQ=255)                                                                                                                       gggCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga     <  7:299929/135‑1 (MQ=255)                                                                                                                          ccgccCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  4:124505/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                   CAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440217‑3440473

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