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

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

CATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAT  >  ACB122/2094687‑2094946                                                                                                                                   |                                                                                                                                  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)                                                                                                                              gcgccgcgccGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAt  <  5:30331/135‑1 (MQ=255)                                                                                                                              gcgcCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCATGACCCAATGCTTGGCCTCGACAGCGGGCAATCAGCCCTCaa  >  2:459708/1‑134 (MQ=255)                                                                                                                                   |                                                                                                                                  CATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAT  >  ACB122/2094687‑2094946

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