Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	3,440,331	T→G	44.5%	A19A (GCT→GCG)	CJ019_03266 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	3,440,331	0	T	G	44.5%	0.6 / 3.7	22	A19A (GCT→GCG)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base G (7/3);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGC  >  ACB143/3440201‑3440435                                                                                                                                   |                                                                                                         aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTggcgg                                                                                                     >  6:64183/1‑136 (MQ=255)             ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGGGCGGCGCCCGGGCCGc                                                                                         >  2:221492/1‑136 (MQ=255)                    gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCTGGgtg                                                                                  >  1:268228/1‑136 (MQ=255)                            cTGGAGGATCGCACGGCCCAACACACGCGAGCGGTCGGGGGGCTGCAGGGCGCGCTGGGGCCCCCCGCGTGGGGGCCGGGGCGGGGGTGGGGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                          <  6:137620/136‑1 (MQ=255)                                aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCCCTGGCGGCCGCCGGGCAGCCTGGGTGTGGCAGGccc                                                                        >  1:64006/1‑133 (MQ=255)                                       cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCtat                                                                >  8:30495/1‑135 (MQ=255)                                        aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCAGGGGGTGGGACGCCCACGCCgata                                                               >  2:66207/1‑135 (MQ=255)                                                  aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCGATAAGGCGCAGcgc                                                    >  1:239426/1‑134 (MQ=255)                                                      gcgagcgGGCGGGGTGCTGCAGTGCGGGCTTGGGCCACCCGCGGGGGCGCCGGGGCGGGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg                                                <  7:70152/136‑1 (MQ=255)                                                         agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGGGGGAGGCCCACCCCGATAAGGCGCAgcgcgcggctc                                             >  2:25689/1‑131 (MQ=255)                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAgg                                          >  7:229787/1‑136 (MQ=255)                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAAGGAGCAGATCGCTGATCatg                                          >  7:207531/1‑134 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCAATAAGGCGCCGCTCGCTGCTCAGGcc                                        >  7:90296/1‑136 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCACCGCCGATAAGGCGCCGCGCGCGGATCAGGcc                                        >  3:241591/1‑136 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCGATAAGGCGCAGCTCGCTGATCAGGcc                                        >  1:9377/1‑136 (MQ=255)                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCCCGGGCGCCGCCCCGGCCGCCAGGGGGTGGCACGCCCCCCCCCATAAGGCGCAGCGCGCGGAACAGGcc                                        >  6:65009/1‑136 (MQ=255)                                                                                 gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCCGGCCCACGCCTATAAGGCGCAGCTCGGTGGGCGGGCCGCTGACTTCGGGCAGCAGt                     >  7:264335/1‑136 (MQ=255)                                                                                 gcTTGGCCCCCCCGCGGGGGCGCGGGGGCGCGGCTGCGGCCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                     <  8:67210/136‑1 (MQ=255)                                                                                          ccccGCGTGGGCGCGGGGGCGCTGCGGCGGGCCGCCCGCGCTGGCGCCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga            <  1:66207/135‑1 (MQ=255)                                                                                              gCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCGGGCGGCGCCCCGGGCGCCGGGGTGTGGGGGGCCCACCCCCATAATGCGCAGCGCGGGGATAAGGGCGCGGACTTCGGGCAGCAGGTGGGGGGGAAggg        >  2:280985/1‑135 (MQ=255)                                                                                                    tCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc  <  7:148665/136‑1 (MQ=255)                                                                                                    gcggcGTGGCGCTGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc  <  3:283288/132‑1 (MQ=255)                                                                                                                                   |                                                                                                         AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGC  >  ACB143/3440201‑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: ‑