Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,094,808	A→C	53.6%	L186R (CTG→CGG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,094,808	0	A	C	53.6%	‑3.6 / 12.4	41	L186R (CTG→CGG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/3);  new base C (19/2);  total (34/7)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB122/2094680‑2094936                                                                                                                                 |                                                                                                                                 ccGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCCACAGTTTCGCAAGATACGTCAGCTGCTCTTCATACACCGACGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCCCAGGCCCCgcgca                                                                                                                            >  6:30331/1‑134 (MQ=255)        cATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCCGCTGCTCTGCATACCCCTTCGGCTCCGGCTGCTTGCGCAAATAATGCACGCCGCGCTCGAAGGCCACGCGCGCCCgggc                                                                                                                     >  2:412613/1‑133 (MQ=255)            cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACCCCTTCGGCTCCGGCTGCTTGCGCCAATACTGCCCGCCGCGCGCGCAGGCCCGGCGCCccgggccgggg                                                                                                                 >  1:125570/1‑134 (MQ=255)                    tACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCGTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGACGGCCCCGCCCCCCGGGCCCGGCGGGCCgtt                                                                                                         >  6:401118/1‑135 (MQ=255)                            cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGGAAATACTGCACGACGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGGTGGGCCTTTTGCAggg                                                                                                 >  4:54290/1‑135 (MQ=255)                            cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCACGCGCCCCAGGCCGGGCTGGCCgtttttcccgg                                                                                                 >  2:339567/1‑129 (MQ=255)                            cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGCCGTTTTGGAggg                                                                                                 >  1:129920/1‑135 (MQ=255)                            cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACCCCGCGCACGAAGGCCCGGCGCGCCTGGCCGGGCAGGCGGTTTTGGAggg                                                                                                 >  6:464496/1‑135 (MQ=255)                                 gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGTGCCTTTTTGCCGGGGCTGt                                                                                            >  6:128180/1‑135 (MQ=255)                                 gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCAGTCACCGTAGGCTCCGGCAGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCAGGCGCGCCGGGCCGGGCAGGGCTGTTTTCAGGGCCGGt                                                                                            >  6:102695/1‑135 (MQ=255)                                  cccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCGGTTTGGAggggcgggt                                                                                           >  3:87073/1‑129 (MQ=255)                                    cGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCCCCGGGGCGGGCGGGCCGTTTTTAAGGGCCGGGTgg                                                                                         >  5:326559/1‑135 (MQ=255)                                     gggCCAACTGTTGCTCAAGATTGGGCAGCGGGTCTTCATACCCGTTCGGCGCCGGCGGCTGGCGAAAATCCTGCCCGCCGGGCTCGAAGGCCAGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGc                                                                                        <  8:362949/135‑1 (MQ=255)                                         cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCGGGCGGGGCGGTTTTGAGGGCCGGGTGGGCgcg                                                                                    >  7:409950/1‑133 (MQ=255)                                              gTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGGGGGCCGTTTTGGCGGGGCGGTTGGCCGAgggggg                                                                               >  4:415159/1‑131 (MQ=255)                                              gTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCCCCGCGCCGGCCGGGGCGGTTTGCACGGCCGGGTGGGCGCGgttgg                                                                               >  3:210652/1‑134 (MQ=255)                                                    cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCGGGCGGGCCGTTTTGCAGGGCCGGTTGGCCGCGGGTGTTgaagc                                                                         >  7:362949/1‑134 (MQ=255)                                                    cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCCCCCGCCCGGCCCGGCGGGCCGATTTTGAGGGGCGGGTGGCCCGGGgtgttgttg                                                                          >  1:407304/1‑131 (MQ=255)                                                    cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGGAGGGCCTGTTGGCCGCGGTTTTTtaaga                                                                         >  4:121463/1‑135 (MQ=255)                                                      aGATCCGGCAGCTGCTATTCAAACACCTTCGGCTCCGGCTGCTTGCCAAAATACTCCCCCCCGCGCGCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAAc                                                                       <  8:380267/135‑1 (MQ=255)                                                      aGATACGTCAGAGGATCTTCAAACACTTGCGGCGCGGGCTGCTTGCGCAAATACTCCACGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                         <  8:350361/133‑1 (MQ=255)                                                               aGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCAGGGCGTGCTGGGCGTTATGAACGGCCTGGTGTGGGAGGTTTTTGAATAACTatatttta                                                              >  6:38116/1‑131 (MQ=255)                                                                     tctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGAAAGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                         <  4:83963/134‑1 (MQ=255)                                                                                   cGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCTGTTTGCAGGGGCTGTTGGGCGAGGGTGGTGAAGAACGCTATATTCAGCAGCCCCCGGATATGGCg                                          >  8:415750/1‑135 (MQ=255)                                                                                   cGGCTCCGGCTGCTTGCGCAAATACCGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCGTTTTGCAGGGCCGGGTGGGCGCGGGTGTTGAAGACCTCTATATGCAGCCGCACCCCGATATGGCg                                          >  6:27053/1‑135 (MQ=255)                                                                                    ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGGCGGGGTGGGCGTTTTTCTGGGGGTGGTGGCCGAGGTTGTTGCTGAACACTATATTGAGCCGCACCAGGatattggcg                                         >  5:269860/1‑130 (MQ=255)                                                                                    ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCAGGCCGTTTTGCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCGATAAGCACCAGCACCACGAAATCGccg                                         >  2:283416/1‑133 (MQ=255)                                                                                    ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCCGGCGGGCCGGTTTGCAGGGGCGGTTGGGCGCGGTTGGTGCAGAACGCTATATTCAGCCGCACCACGCTATGGc                                           >  5:189839/1‑133 (MQ=255)                                                                                      cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGGCTTTTTGGAGGGGCTGGTGGCCGTGGTTGTTGAAGAACTCGATATGCAGCAGCACCACGATATTGCGGag                                       >  4:180622/1‑134 (MQ=255)                                                                                      cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAt                                       >  6:202269/1‑135 (MQ=255)                                                                                             tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGGTGGCCTTTTTGCAGGGCCGGGTGGGCGAGGTTGTTGAAGAACTAGATATGGAGCAGCACCAGGATATGGCGGCTTTCCTTg                                >  4:70040/1‑135 (MQ=255)                                                                                             tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCTTTTTGCCGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAtttacgtg                                >  2:276374/1‑130 (MQ=255)                                                                                                                 gccgcGCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGCACGCTATATGCAGCCGCCCCAGGATATGGCGGATTTCCTTGACCCAATGCGTGGCcgcg              >  1:259270/1‑133 (MQ=255)                                                                                                                 gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGGCGTTTTGGCGGGCCGGTTGGGCGAGGTTGGTGAAGAACTCTATATGCAGCAGCACCAGGATATGGGGGATTTACTTTACCCAATGCTTGGCcgcg              >  8:16285/1‑133 (MQ=255)                                                                                                                 gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCGTTTTGCCGGGGCGGTTGGCCGCGGTTGTTGAAGAACTCTATATGCAGCCGCCCCCGGATATGGCGGGTTTCCTTGCCCCAATGCTTGGCCTCgcc            >  2:80878/1‑133 (MQ=255)                                                                                                                 gccgcGCTCGAAGGCCAGGCGCGCCAGGCCGGGCGGGCCGTTTTGCCGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATGCCGCCGCCCCCGGTTATGGCGGATTTCCTTTACCCAATGCTTGGCcgcgcc            >  2:195382/1‑133 (MQ=255)                                                                                                                  ccgcGCTCCCCGGCCGGGCGCGCCGGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa           <  3:54290/135‑1 (MQ=255)                                                                                                                   cgcgCTCGAAGGCCCGGCGCGCCTGGCCGGGCGGGCCTTTTTGCAGGGCCTGTTGGCCGCGGTTGGTGAAGACCTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAg          >  3:210354/1‑135 (MQ=255)                                                                                                                      gcTCGAAGGCCACGCGCGCCAGGCCGGGCTGGCCTTTTTGCAGGGCCTGGTGGCCGCGGTTTTTGTAGAACTCTAATAGCAGCAGCACCAGGATATGGCGGATTTACTTTCACCCATGCTTGGGCTCGACAGCgc       >  2:139860/1‑134 (MQ=255)                                                                                                                           atggCCCGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  2:314897/133‑1 (MQ=255)                                                                                                                           acggCCTGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  8:409950/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                                 CCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB122/2094680‑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: ‑