Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,721,874	A→C	34.9%	W97G (TGG→GGG)	CJ019_01617 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,721,874	0	A	C	34.9%	38.2 / 2.1	32	W97G (TGG→GGG)	CJ019_01617	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/4);  new base C (6/6);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTG  >  ACB143/1721741‑1722001                                                                                                                                      |                                                                                                                                tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGc                                                                                                                               >  8:219954/1‑136 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGc                                                                                                                               >  3:99994/1‑136 (MQ=255)  cgcgGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGc                                                                                                                               <  8:168267/135‑1 (MQ=255)       gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCg                                                                                                                         >  3:153618/1‑136 (MQ=255)           tcccATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCTACCCCGTCTTCCTCGCCCCCGCCCCCCGgcgc                                                                                                                     >  2:275351/2‑136 (MQ=255)                  ccTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCCTGAGGGTATTCCCCGTCTTCCTGGCCCCCGCCCCCCGGCgcgccgccg                                                                                                              >  1:297735/1‑134 (MQ=255)                  ccTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCgcgcagcg                                                                                                               >  3:22875/1‑135 (MQ=255)                    tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTTATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCCCCCGGCGCCCAGCGCAc                                                                                                            >  7:137592/1‑136 (MQ=255)                    tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCCTCTTCTTACCCGTCTTCCTAGCCCCCGCCACCCGGCGCGCCGCGCAc                                                                                                            >  4:197051/1‑136 (MQ=255)                                 aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCATCCCCGTCTTCCTCGACCCCTCCACCCGGCGCGCATCGCACTTCCACACCCAgc                                                                                               >  4:257088/1‑136 (MQ=255)                                 aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACCCCCAgc                                                                                               >  2:100332/1‑136 (MQ=255)                                                tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATgcggcg                                                                                >  3:229433/1‑136 (MQ=255)                                                        gCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACATCCACACCCAGCGCGGCAAATGCGGCGCGGgcg                                                                          >  5:261268/1‑134 (MQ=255)                                                              tGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAAc                                                                  >  3:147326/1‑136 (MQ=255)                                                                aTTTCCTCTTCACCGTCGGCGCTGATGCCGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCTACCCAGCCCACCGGCGCGCAGCGAACTTCCCAACCCAGCGCGGCAAATGAGGCGCTGGCGCAGGCAACGt                                                                >  8:103617/1‑136 (MQ=255)                                                                    ccTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCCCGCCGCGCACTTCCCCACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGGCATc                                                            >  2:105246/1‑136 (MQ=255)                                                                     ctctTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCg                                                           >  5:131935/1‑136 (MQ=255)                                                                                          gcgccccccGCGGTCATCTTCTTCTCCCCCTTCTCCCCCCCCCCCCCCCCCCGGCGCCCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                       <  4:54265/130‑1 (MQ=255)                                                                                           cgcccccAGCGGGAAGAATGTTCTCCCCCGTATCCCTCCCCCCCCCAACCGGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                       <  4:216341/130‑1 (MQ=255)                                                                                                ccAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTc                                    >  5:143335/1‑132 (MQ=255)                                                                                                 ccgCGGTCAGAATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTCGTCGCTTTCCAGAt                               <  7:219954/134‑1 (MQ=255)                                                                                                     ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTCTGGGCGCGGTCCAGATACAg                           >  4:10112/1‑136 (MQ=255)                                                                                                         atcATCTTCTCCCCCTTCTCCCTCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTg                       <  1:203653/134‑1 (MQ=255)                                                                                                         aGCATCTTCTACCCGTTCTTCCTCCACCCACCCCCCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTg                       <  1:84349/136‑1 (MQ=255)                                                                                                                    ccccGTCTTCCTCCCCCCCCCCCCCCGGCGCGCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc             <  4:229433/135‑1 (MQ=255)                                                                                                                    ccccGTCTCCCGCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc             <  3:239467/135‑1 (MQ=255)                                                                                                                    ccccGTCTCCCCCCCCCCCCCCCCCCGGCGCGCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCATGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc             <  4:99994/135‑1 (MQ=255)                                                                                                                        gTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCCACGTCATCGGCCCACGGCGTCTCGTCGCTTTCCAGATAAATGCTGTTCCAATTCCCGACa        >  2:42117/1‑136 (MQ=255)                                                                                                                        gTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTa          >  2:211791/1‑134 (MQ=255)                                                                                                                        gTCTGCCTCGCCCCACCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACa        <  8:201039/136‑1 (MQ=255)                                                                                                                          cTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTGCATTTACCTACAGc      >  2:188852/1‑136 (MQ=255)                                                                                                                              cTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTg  >  4:225841/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                                TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTG  >  ACB143/1721741‑1722001

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑