Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	3,440,326	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
*	ACB143	3,440,326	0	A	C	78.1%	15.6 / 2.1	33	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (20/5);  total (26/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

ACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGC  >  ACB143/3440198‑3440435                                                                                                                                 |                                                                                                              aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTg                                                                                                         >  8:50425/1‑135 (MQ=255)   cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGggc                                                                                                       >  2:197689/1‑135 (MQ=255)   cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTggcg                                                                                                      >  6:167172/1‑136 (MQ=255)   cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                       >  8:120336/1‑135 (MQ=255)   cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCCCTggc                                                                                                       >  1:109353/1‑135 (MQ=255)              aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGc                                                                                            >  2:399224/1‑135 (MQ=255)              aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGcc                                                                                           >  5:167978/1‑135 (MQ=255)                              tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCAg                                                                            >  7:358353/1‑135 (MQ=255)                                tGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGcc                                                                         >  5:273943/1‑136 (MQ=255)                                tGGAGCAACGCACCGCCCGCTCCCCGCGCGCGCTCGCGGGGCTGCCGGGCGCGCGGGGGCCCCCCGCCGGGGCGCCGGGGCGCGGCGGCCGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                          <  4:205436/135‑1 (MQ=255)                                 ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGc                                                                          >  8:124662/1‑134 (MQ=255)                                   aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCCGCCTGGGTGTGGGAGGCCAAc                                                                      >  6:148936/1‑136 (MQ=255)                                          cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGccga                                                                 >  4:185263/1‑132 (MQ=255)                                          cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGGCGGGGTGTGGCCGGCCAACGccga                                                                 >  7:193504/1‑132 (MQ=255)                                           aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGccga                                                                 >  3:62671/1‑131 (MQ=255)                                           aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCgat                                                                >  4:303652/1‑134 (MQ=255)                                               cccTCTACCCGCGCGCGCTCGGGGTGCGGCTGGGGGCGCGGGGCCCCCCCGCGGGGGCGGCGGGGCGGGGGGGCGGCCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc                                                           <  4:367013/135‑1 (MQ=255)                                                     aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCGGGCGGCGGCCGGGCCGCCCGGGGGGTGGAGGGCAACGCCGATAAGgcggcacg                                                     >  8:142024/1‑130 (MQ=255)                                                               gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCAATAAGGCGCCGCTCGCTGAGCAg                                           >  7:338524/1‑135 (MQ=255)                                                               gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCAGATCAgg                                          >  6:118800/1‑136 (MQ=255)                                                                 tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCAGGCGGCGCCCCGGCCGCCTGGGGGTGGCCGGCCCACCCCGATAAGGCGCAGCGCGCTGATCAGGc                                         >  3:218586/1‑135 (MQ=255)                                                                 tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCGGATCAGGc                                         >  1:219811/1‑135 (MQ=255)                                                                 tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCACCGCCGATAATGCGCAGCTCGCGGATAAGGc                                         >  4:331923/1‑135 (MQ=255)                                                                 tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCGGATCAgggc                                        >  6:382196/1‑134 (MQ=255)                                                                             aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGTGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCgg                             <  3:61226/135‑1 (MQ=255)                                                                                  gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCAGGCGGGGCCCGGGCAGCCTGGGGGTGGGAGGCCCACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGTCTTCGGCAAg                         >  2:340934/1‑134 (MQ=255)                                                                                     cgtgcgCCCCCCGCGTGGCCGCGGTGGGGCGGCCGCTGGCCGCGCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                     <  3:146545/130‑1 (MQ=255)                                                                                         ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGGCCGGGCCGCCGGGGGGGGGCAGGCCAACACCAATAAGGCGCCGCGCGCGGATCAGGCCGCTGACGGAGGCCAGCAGTTggg                 >  8:129447/1‑134 (MQ=255)                                                                                                 gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGGAAGCAGTTGGCGGAGAAgg         >  7:19236/1‑135 (MQ=255)                                                                                                 gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCGCACGCCTATAAGGCGCAGCGCGCGGCTCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGaggagg         >  3:370455/1‑135 (MQ=255)                                                                                                       gcGCCGGGGCGCGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc  <  6:337016/135‑1 (MQ=255)                                                                                                        ggcggTGGCGCGGGTGCGCGCCGCCTGCGCTGGCGGCGCCCTGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc  <  7:37954/131‑1 (MQ=255)                                                                                                        cggggTGGCGGGGCGGCGGCCCGCGCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc  <  3:303652/131‑1 (MQ=255)                                                                                                                                 |                                                                                                              ACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGC  >  ACB143/3440198‑3440435

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