Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	5,507,093	T→G	54.5%	L36R (CTC→CGC)	CJ019_05179 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,507,093	0	T	G	54.5%	‑2.4 / 9.2	25	L36R (CTC→CGC)	CJ019_05179	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/8);  new base G (6/6);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

ACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGCC  >  ACB143/5506965‑5507226                                                                                                                                 |                                                                                                                                      acCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCCGCCAGcgcgcaagc                                                                                                                                >  4:25667/1‑131 (MQ=255)        tCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCCGCCGGCCGGATCGAAGGCTGGgcag                                                                                                                         >  5:222273/1‑136 (MQ=255)           aGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTCCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGGAGGCCGGCCGGCGCGCAGGCGGGgcagccc                                                                                                                      >  3:274811/1‑134 (MQ=255)                     gCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCGCGAAGGCGGGGCAGCCCAAcccccttc                                                                                                            >  6:139928/1‑131 (MQ=255)                     gCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCCAGCCAGCCCGCCTGCTAGCAAGGTGGGGCGCACAAcccccttc                                                                                                            >  2:217707/1‑131 (MQ=255)                                            ctcatGCTGCTAATCCCCGTGGCCACCCACGCGAGCGAGGGGTTTTTCTTTCACCCCGCCAAGCCCGCAGCGCAGGCCGGCGTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGaaa                                                                                     <  4:261737/136‑1 (MQ=255)                                            cgcctGCTGCTAACCCCCCTGGCAACCCCCCCCAGCGCGGGCTAATTCTTTCTCCCCCCCACGCCCGCCGCGCGGGCCGGCCTGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGaaa                                                                                     <  8:190128/134‑1 (MQ=255)                                                 gctgctAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCGGGCGCGCCCGCGCGGACCCTGGGGGGGACCACCCCATGCCCCGCCCCGGCCCCGGCGGGGAAAccag                                                                                >  5:273986/1‑134 (MQ=255)                                                   tgctAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCCGCCGGCCTGCGAGAACTCAGGGCAGCACAACCCGCAGCCGCCCCCCACCCGCAGCTGGaaaaacaacg                                                                              >  5:101516/1‑130 (MQ=255)                                                          cccccTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCGCGAAGGCGGGGCCGCCCAACCCCCTGCCGCGCGCCGGCCCGGGCTGGGAAACCTGCGCCCAGGc                                                                       >  7:281250/1‑136 (MQ=255)                                                              cggccAACCCCCCCCTGCGTGGGCTAATTCTTCCCCCCCGCCACCCCCGCCGCGCCGGCGGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCa                                                                   <  5:195872/132‑1 (MQ=255)                                                                       cACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCGCGCAGGCGGGGGCGCCCCACCCCCTGCCCCCCGCCAGCCGCGGCGGGGAAAACAGCGCCAACGCCCCCTCGccccc                                                           >  6:195073/1‑135 (MQ=255)                                                                        aCGCCAGCGAGGGCTAATTATTTCACACCGCCAACCCCGCCGCCCCGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCaa                                                         <  1:217707/136‑1 (MQ=255)                                                                                gtgggCTAATTCTTTCTCACCGCCAAGCCCGCCGAGCTGGCCGGCCGGCCCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAgc                                                 <  6:124950/134‑1 (MQ=255)                                                                                   ggctatttcttcccCACCCCCACGCCCGCCGAGCCGGCCGGCCGGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCa                                              <  2:21654/130‑1 (MQ=255)                                                                                   gcctaaTGCTATCGCACCGCCACCCCCGCCGCGCAGGCCGGCCTGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCa                                              <  4:34170/132‑1 (MQ=255)                                                                                               ccccccgccccgccCGCCGCGCGGGCGGGCCTGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTg                                   <  5:151562/131‑1 (MQ=255)                                                                                                         cccccGCCGCGCCGCCCGGCCGGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCt                        <  4:251288/134‑1 (MQ=255)                                                                                                             cgccgcgccgCCCGGCCGGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTgcgc                    <  4:272036/136‑1 (MQ=255)                                                                                                             cgccgcGCCCGCCGGCCCGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTgcgc                    <  5:267095/136‑1 (MQ=255)                                                                                                                       gccggccgGCTCGAAGGCTTGGCAGCACACCCCGATGCGGCGCGCATGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCTGCCTGATCg          <  7:171930/136‑1 (MQ=255)                                                                                                                       gccggccgGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCg          <  3:245102/136‑1 (MQ=255)                                                                                                                       gccggccgGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCg          <  5:276509/136‑1 (MQ=255)                                                                                                                       gccggccgGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCg          <  1:39194/136‑1 (MQ=255)                                                                                                                       gccggccCGCCCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCg          <  3:280071/136‑1 (MQ=255)                                                                                                                               gCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGcc  >  2:66943/1‑136 (MQ=255)                                                                                                                                 |                                                                                                                                      ACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCTGCGCGGCCTGATCGACAACGCC  >  ACB143/5506965‑5507226

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