Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,247,863	T→G	30.9%	R255R (AGG→CGG)	CJ019_04943 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,247,863	0	T	G	30.9%	44.9 / 3.7	47	R255R (AGG→CGG)	CJ019_04943	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base G (12/4);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCT  >  ACB122/5247728‑5247992                                                                                                                                        |                                                                                                                                  gcTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCt                                                                                                                                   <  8:26660/136‑1 (MQ=255)      cgcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGc                                                                                                                                <  2:33510/134‑1 (MQ=255)      cgcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGa                                                                                                                              <  7:132816/136‑1 (MQ=255)                  cccTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTGCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCGGGCGAGTTCTGCtgtg                                                                                                                   <  3:151045/135‑1 (MQ=255)                     tGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGACGTGCCTGGCGAGTTCTGCTGTGGGc                                                                                                                >  4:137543/1‑135 (MQ=255)                        gcgcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACCCCACGCGCGTGAAGGGCCGGGCGCGTTCTGCGGGGGGCgcc                                                                                                             >  2:169278/1‑133 (MQ=255)                          gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTg                                                                                                           >  4:138226/1‑135 (MQ=255)                          gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGGCGGGCGGGGTCGGCGGTGGGCGAc                                                                                                             >  1:158935/1‑133 (MQ=255)                               ggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCATGTAGTGGCGGGCGGGGTGTGGTGTGGGCGGCTGATg                                                                                                        >  6:336455/1‑133 (MQ=255)                               ggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGa                                                                                                     <  7:205532/136‑1 (MQ=255)                               ggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTGTGCTGTGGGGGACGGATgg                                                                                                       >  6:106625/1‑134 (MQ=255)                                gacagcTACAGGGATGCCAGCCTGGCTGTCCTGTACCGCTTCCGCTGCTTCTCCAGCCAGTAAAAAGGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCGTGGCGAGTTCTGCTGTGGGCGACTGATGGGa                                                                                                     <  5:108919/131‑1 (MQ=255)                                              tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCCCCAGCCAGTAACCACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGGGGCGGGCGCGTGCGGCTGTGGGCGCCGGCTGGGGACACCGCCCCAggg                                                                                       >  5:262715/1‑135 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCCCCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGTGAAGGGGCGGGCGAGTGCTGCGGGGGGCGGCGGCTGGGAAAACAGCACCAGGGCGt                                                                                    >  1:135315/1‑135 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGGGTGAAGTGGCGGGGGGGGTGTGGGGTGGGCGGCTGCTGGGGAAAACGCCCCCGGGCgg                                                                                    >  5:6582/1‑134 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGGGAGGTCTGCTGTGGGCGCATGCTGGGGAAACAGCACCAGGGctt                                                                                    >  2:314895/1‑133 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGGTGTGGGCGCCTGATGGGAAAACCGCACCAGGGCGt                                                                                    >  3:296517/1‑135 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGCAGTGGCGGGCGCGGTCTGCTGTGGGGGCCTGGTGGGGAAACAGCCCCAGGGCGt                                                                                    >  3:274003/1‑135 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGCAGTGCCTGGCGAGTTGTGGTGTGGGGGGCGGAGGGGAAAACAGCACCAGGGCGt                                                                                    >  5:154171/1‑135 (MQ=255)                                                 cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCGGTGGGCGCCGGATGGGGAAACAGCCCCAGGGCGTc                                                                                   >  8:76883/1‑136 (MQ=255)                                                   gccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCg                                                                                 >  8:220465/1‑136 (MQ=255)                                                   gccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCACGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACAGATGGGCAAACAGCACCAGcgcggta                                                                                  >  6:108919/1‑132 (MQ=255)                                                     cagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCACGCGCGTGAAGTGCCGGGCGCGTTCTGGGGTGGGCGCCGGATGGGCAAACAGCACCCGGGCGTCCgg                                                                                >  6:386569/1‑135 (MQ=255)                                                         cggcccTGTACCGCTTCCGCTGCTTCACCAGCCATTCCCAACGGGCATTCCAGCAGCGACGCCAGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                                           <  8:67505/132‑1 (MQ=255)                                                          tGTCCTGTACCGCTTCCGCTGCTTCACCAGCCACTAACAACGGGCATTCCAGCATCGACGCCAGGGCGGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                                           <  5:95696/135‑1 (MQ=255)                                                              cTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTccgccg                                                                       >  1:322545/1‑135 (MQ=255)                                                               tGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGTAGTGCCTGGCGACTTCTGCTGTGGGCGCCGGATGGGAAAACCGCACCACGGCGTCCGGCTTACGCCGtt                                                                     >  7:371086/1‑136 (MQ=255)                                                                         ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGGGTGAAGTGCCGGGGGAGTTCTGGTGTGGGCGACTGATGGGGAAACAGCCCCAGGGCGTCCGGCGTCCGCCGGTCGCAGcccc                                                            >  2:15021/1‑134 (MQ=255)                                                                         ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGGTTGAAGTGCCTGGCGAGTTCTGCTGTGGGGGACTGATGGGCAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGAccc                                                            >  5:208229/1‑134 (MQ=255)                                                                         ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGCGGTCTGCTGTGGGCGCCGGAGGGGAAAACAGCCCCCCGGCGTCCCGCGTCAGCCCTTCGCaccccc                                                            >  6:329290/1‑130 (MQ=255)                                                                         ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGGGGGCGACTGAAGGGAAAACAGCACCCCGGCGTCCGGCGTGAGCCCTTCGCAGCCCa                                                            >  3:94179/1‑135 (MQ=255)                                                                         ccGCTGCGTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGGTTCAACTGGCGGGAGTGTTGTGCTGTTGGCGACAGCTGGGTAACCAGCACCAGGGCGACCTGCGTCCCCCGTTCGGACACCAg                                                           >  8:186281/1‑136 (MQ=255)                                                                          cGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                                          <  7:64659/136‑1 (MQ=255)                                                                           gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                                          <  3:138226/135‑1 (MQ=255)                                                                                   aCCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTc                                                  <  4:218805/135‑1 (MQ=255)                                                                                      agccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGc                                              >  8:103934/1‑136 (MQ=255)                                                                                               caacGGGCATTCCAGCATCGACCCCCGCGGCGTGAAGTGCCTGGCGGGTTATTGTGTGGGGGGCTGATGGGAAAACAGCCCCCCGGCCTGCGGCTTCAGCCGTTGGCAGACCAGCGCAAGCTGCGGCAGCCGcag                                      >  4:167768/1‑133 (MQ=255)                                                                                                     gCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCTGGGGCGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCCGCGTCAGCCGTTCGCAGCCCAGCGCCAGCGCCGGCCGCGGGGGCCCg                                  >  6:326463/1‑133 (MQ=255)                                                                                                             gCATCGACGCCAGGCGCTTCAAGTGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACACCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGGTCGCCAGCgg                       <  7:186281/136‑1 (MQ=255)                                                                                                              cATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCgg                       <  5:106625/135‑1 (MQ=255)                                                                                                                   aCGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATc                   <  2:336191/134‑1 (MQ=255)                                                                                                                      ccAGGCGCTTGAAGTGCCTGGCGCGTTATGCTGTGGGCGACTGATGGGGAAACAGCACCAGGGCGTCCGGCTTCCGCCGTTCGCCGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCCCCGGGCTCCcgcg              >  7:328221/1‑136 (MQ=255)                                                                                                                      ccAGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGTGGGCGACGGATGGGTCAACAGCCCCACGGCGTCCCGCGTGAGCCGTTCGCAGACCAGCGCAACCGCCTGCACCGGCTGCCCGGTCGCCAGCCGGCTCAcgc               >  6:168450/1‑135 (MQ=255)                                                                                                                              ttGACGGGCCTGGGGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcg       <  2:395644/135‑1 (MQ=255)                                                                                                                              ttGAAGGGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcg       <  5:326159/135‑1 (MQ=255)                                                                                                                                  tgTGCGTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCt  <  8:339455/135‑1 (MQ=255)                                                                                                                                  aGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCt  <  7:220465/136‑1 (MQ=255)                                                                                                                                   gTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACTCGTGTTTCGTCGCt  <  3:179292/135‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCT  >  ACB122/5247728‑5247992

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