Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,094,817	T→G	37.0%	Q183P (CAG→CCG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,094,817	0	T	G	37.0%	10.8 / 3.2	26	Q183P (CAG→CCG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (8/7);  new base G (3/5);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB143/2094701‑2094936                                                                                                                       |                                                                                                                         ccTGGCCATTGCCCGGGGCAAACGGTTTCCCATCATACTTCAGCGGCTCTCCATACACCTTCGGCGCCGGCTGCTGGCGAAAAAA‑CTGCCCGCCGCGCTCGACGGCCAGGCG‑CGCGTGGCCGGGC‑TGGCCTTTTTg                                                                                                      <  2:241841/136‑1 (MQ=255)       cttgggCCGGGGCCAACAGTTTCTCAAGATCCGCCAGCTGCTATTCATACACCTTCGGCTCGGGCTGGTTGCGCAAATA‑CTCCCCGCCGCGCGCGACGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGc                                                                                                <  4:268737/132‑1 (MQ=255)       cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCCCGCCGCGCTCGCAGGCCAGGCG‑CGCCTGGCCCGGG‑TGGCCGTTTTGCAGGGc                                                                                                >  7:69396/1‑136 (MQ=255)       cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCGGCGCCCCCCGGCCGGGG‑GGGCCTGTTTG‑AGGGc                                                                                                >  1:101195/1‑136 (MQ=255)       cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCGGCG‑CGCCTGGGCGGGC‑TGGCCGTTTTGGCGGGc                                                                                                >  6:214471/1‑136 (MQ=255)       cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑ATGCACGCCGCGCTCGAAGGCCCTGCG‑CGCCGTGCCGGGC‑TGGGCGTTTTgaacggg                                                                                                >  8:273290/1‑130 (MQ=255)       cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCATGCGCAAATA‑CTGCACGCCGCGCTCGCACGCCACGCC‑CGCCCGGCCGGGCGGGGCGTTTTTCCAggc                                                                                                 >  2:295306/1‑135 (MQ=255)            gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCCGCG‑CGCCGGGGCCGGC‑GGGCCGTTTTGGAgggggggtt                                                                                           >  7:206490/1‑130 (MQ=255)                     aaCAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCCGCG‑CGCCTGGCCGGGC‑TGGCCGTTTTTCAGGGCCTGTTGGGCGAGgt                                                                                  >  7:179708/1‑136 (MQ=255)                         gTTTCTCAAGATCCGTCAGCTGCTATTCATACACCTTCGGCGCCGGCTGCTTGCGCAAATC‑CGCCACGCCGCGCGCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                              <  7:238437/136‑1 (MQ=255)                          tttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCGGCG‑CGCCGGGCCGGGC‑TGGCCGTTTTGCAGGGCCGGTTGGCCGCGGgtggtg                                                                             >  7:221149/1‑136 (MQ=255)                              tcAACATACGACCGCTGCTCTTCATCAAACTTCCGCTCCGGATGATTGCCAAAATA‑GTGCACGCCGCGAGCGACGGCCAGGCG‑CGCGTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                         <  8:148442/136‑1 (MQ=255)                                         cAGCTGCTCTCCATACACCGTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCCCGCCTCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGTCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                               <  6:191768/135‑1 (MQ=255)                                           gctgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCCGGCG‑CGCCAGGCCCGGC‑TGGCCTTTTTTCAGGGCCAGGTGGCCGAGGGTGTTGAAGAACTCTATATgcag                                                            >  2:78633/1‑136 (MQ=255)                                                             tgggCTCCGGCTGCGTGCTCAAATA‑GAGCCCGCCGCGCGAGCAGCACCGGGG‑CCCCGGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                           <  8:254095/133‑1 (MQ=255)                                                               ggCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCACGGCCGGTTGGCCCAGGTTGTTGAATAACTCTAAATGCAGCCTCACCAGGATATTGCgg                                         >  2:115352/1‑135 (MQ=255)                                                                 cTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGCAGGCCCGGCG‑CGCCTGGCCGGGC‑TGGCCTGTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCCGCCGCACCACGATATGGCGGAtt                                      >  8:124492/1‑136 (MQ=255)                                                                 cTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCTATATGCCGCAGCACCAGGATATGGCGGAtt                                      >  8:236074/1‑136 (MQ=255)                                                                    cGGCTGCTTGC‑TAAATACCTCCCCGCCGCGCGCGACGGCCTGGCG‑CGCCGGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                     <  5:251912/134‑1 (MQ=255)                                                                          gttGCTCAAATA‑GGGCCCGCCGCGCTCTAAGCCCCGGCG‑GGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc                              <  6:144683/134‑1 (MQ=255)                                                                             gctaaaaTC‑CTCCCCGCCGCGCGCGACCGCCCGGCG‑CGCCGGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                          <  4:208658/132‑1 (MQ=255)                                                                              cgcAACTC‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                          <  2:67740/135‑1 (MQ=255)                                                                                   aTA‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑GGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGCACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTg                    >  6:85812/1‑136 (MQ=255)                                                                                           gcgccgcGCTCTACGGCCCGGCG‑CGCCGGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa             <  1:74384/135‑1 (MQ=255)                                                                                                  gcTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCggg      >  3:150104/1‑136 (MQ=255)                                                                                                      cacggCCGGGCG‑CGCCGGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGGTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  6:272603/133‑1 (MQ=255)                                                                                                                       |                                                                                                                         CCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATA‑CTGCACGCCGCGCTCGAAGGCCAGGCG‑CGCCTGGCCGGGC‑TGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB143/2094701‑2094936

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