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	86.7%	28.3 / ‑2.1	31	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base C (21/5);  total (23/8)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTG  >  ACB143/3440201‑3440445                                                                                                                              |                                                                                                                        aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                               >  8:148275/1‑136 (MQ=255)     gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCAGGCGGCGcc                                                                                                           >  7:310189/1‑136 (MQ=255)          aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCTGGCGGCGCCCGGGc                                                                                                      >  5:39905/1‑136 (MQ=255)           aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGcc                                                                                                     >  6:111278/1‑135 (MQ=255)                  tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCGGGCGGCGCCCGGGCCGCCgggg                                                                                              >  6:2421/1‑136 (MQ=255)                    gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCGGGCGGCGCCCGGGCCGCCTGGgtg                                                                                            >  1:230536/1‑136 (MQ=255)                              ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCCGGCCGCCAGGGGGTGGCAGGcc                                                                                   >  3:82275/1‑135 (MQ=255)                                        aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACCCCGATaa                                                                        >  7:357554/1‑136 (MQ=255)                                                         agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCGCCGGCAGCGGGGGGTTGGGAGGGCAACGCCTCTCAGGGCCAGCTCGCGGATc                                                       >  5:313269/1‑136 (MQ=255)                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCTATAAGGCGCAGCTCGCGGATCAg                                                     >  3:317535/1‑135 (MQ=255)                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCAGCCTCGGGGGGGCCGGCCCCCGCCCATAATGCGCCGCGCGCTGATCac                                                     >  4:163677/1‑134 (MQ=255)                                                            gcTCGAGGTGCTGCAGTGAGCGCTAGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGGAGGCCAACGCCGAAAATGCGCACCTCGCTGATCAg                                                     >  4:114214/1‑135 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGGGGCCCACGCCGAGAATGCGCAGCTCGCGGCTCAGGcc                                                  >  6:213329/1‑136 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCCGCTGGGGGTGGGGCGGCCCCCCCAATAAGGCGCCGCGCCCGGCTCAGGcc                                                  >  1:352322/1‑136 (MQ=255)                                                                 aGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCCGGCAGCCGGGGTGTGGGAGGCCCACGCCTATAAGCCGCCGCTCGCTGATAAGGCCg                                                 >  2:79747/1‑134 (MQ=255)                                                                   gTGCGGCTGTGCGCGCGTGCCCCACCCGCCTGGGCGCCGTGGCGCTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATg                                              <  4:354873/135‑1 (MQ=255)                                                                         cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGATGACTTCgg                                       >  5:301871/1‑136 (MQ=255)                                                                         cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCCGGCCCACGCCGATAAGGCGCCGCTCGCGGATCAGGCCGCTGACGTCgg                                       >  5:363578/1‑136 (MQ=255)                                                                                 gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCGGGGGGCGCCCGGGGCGCCGGGGGGGGGCACGCCCACCCCGATAAAGCGCCGCTCGCGGATCAGGGCGCGGACTTCGGCAAGCAgg                               >  6:306858/1‑135 (MQ=255)                                                                                 gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCACCGCCGATAAGGCGCAGCGCGCGGATCAGGCCGATGACTTCGGCAAGCAgg                               >  7:202927/1‑135 (MQ=255)                                                                                   ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCCAGCAGTTg                             >  1:264046/1‑136 (MQ=255)                                                                                        ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCAACGCCGATAAGTGGGAGCTGGCTGGTCCGGCCGCTGACGGCGGGCAGCAGGTGggggg                        >  6:94614/1‑135 (MQ=255)                                                                                              gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCCCCCCGATAAGGCGCCGCGCGCGGATCAGGCCGATGACGTCGGCACGCAGTTGGCGGAGCAGGa                  >  6:142087/1‑136 (MQ=255)                                                                                                    ggggggTGGGGCGGGTGCTGGCCGCCCGGGCGGGGGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc            <  5:301064/131‑1 (MQ=255)                                                                                                    gcggggTGGGGCTGCGGCGGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc            <  8:202927/131‑1 (MQ=255)                                                                                                     ggCCGTGGCGCTGCTGCGGGCCGCCCGCGGGGGGGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc            <  3:313662/134‑1 (MQ=255)                                                                                                     cGCCGGGGCGCGGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc            <  3:127263/135‑1 (MQ=255)                                                                                                         ggggCGCGGCTGCTGGCCGCCCGCGCGGGGGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg       <  7:9624/134‑1 (MQ=255)                                                                                                         gTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGGAGGCCAACGCCGATAAGGCGCAGCTAGCGGATCAGGCCGATGACTTCGGCAAGCAGTTGGGGGAGAAGGATCGCGCGCACg       >  6:252290/1‑136 (MQ=255)                                                                                                          gggCGGTGCGGCGGGCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg      <  1:253446/135‑1 (MQ=255)                                                                                                              gctgcggcgGGCCGCCAGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg  <  5:61850/136‑1 (MQ=255)                                                                                                                              |                                                                                                                        AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTG  >  ACB143/3440201‑3440445

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