Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,902,196	0	T	G	57.1%	7.5 / 18.5	21	V167G (GTG→GGG)	CJ019_03691	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base G (0/12);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATG  >  ACB122/3902064‑3902324                                                                                                                                      |                                                                                                                                 gcgcCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTgcg                                                                                                                               >  7:72403/1‑136 (MQ=255)                       ccaGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAg                                                                                                         >  4:36455/1‑136 (MQ=255)                                           tgCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGCCGATGCAGCGCTGGTGCGGCGCCTGAGCAACGCCCTGCAGCGCCCCTGCCGTGAAACgg                                                                                      >  5:9632/1‑135 (MQ=255)                                            gCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTc                                                                                    >  2:150163/1‑136 (MQ=255)                                                ccGCTACCGCCTGAACACCGCCACCAACCCGAGACGCGCCAGCAAGTTGGGCACCGGGCCCG‑CGATGGCGACGATGCCGGGCGGGGGCGGCGCCGGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGAc                                                                                 <  6:178744/134‑1 (MQ=255)                                                   cTACCGCCTGACACCCTCCACCAACCCCAGCACCGCCAGAAACTTGGCCACCGGCCCCG‑CGAAGGCTACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGc                                                                             <  4:179790/136‑1 (MQ=255)                                                      ccGCCTCATCAACTCCACCAACCCGACCCGCGCCAGCAACTTTGCCAACTGGCCCG‑CGAAGGCGACGAGGCAGCGCTGGGGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCt                                                                            <  8:190031/134‑1 (MQ=255)                                                                aaCTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAg                                                                 >  5:225246/1‑135 (MQ=255)                                                                      agcaACCCGCGCCGCCCCTGCTAGTTTGGCAACGGGCCCG‑CGACGGCGACGATGCCGGGCGGGGGCGGCGCCGTACCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGa                                                          <  4:230278/136‑1 (MQ=255)                                                                       ccaACCCTCGCAGCCCCAGAAAGTTTCGCAACGGGCCCG‑CCCCGGCCAGTTTGCGGGGGGGGGGCGGGGCCTACCCCACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGa                                                          <  6:229536/134‑1 (MQ=255)                                                                                gcagcGCCAGAAACTTTGCAACCTGCCCCG‑CGACGGCGCCGAGGCAGCGCGGGGGCGGCGCCGGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                <  3:243462/136‑1 (MQ=255)                                                                                cccgccccaccaAC‑TTGGCAACCGGGCCGCCTAGGGCGCCTAGGCCGCGCGGGGGGGGCGCGTCACCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                <  8:1140/130‑1 (MQ=255)                                                                                cccccgcCAGCAACTTTGCCAACGGCCCCG‑CGAAGGCGACGATGCAGCGCTGGGGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                <  4:62945/132‑1 (MQ=255)                                                                                       cAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCCGCGCCCCGGCCTTGAACCGGTCGCCTGGCTGTCACATCAAGACCAGGACCTGGCCCGGGAACAGt                                         >  4:181/1‑136 (MQ=255)                                                                                            tgtttGGCCACTGGCCGG‑CGATGGCTACTATGCAGGGGGGGTGGGGCGCGGTAGCAACGAGCTGCTGCGCCCCTGCCTTGAAAGGGTCGACTGGCTGTCACATCTAGACCGGGACCTGGCCCGGGATCAGTTCat                                     <  6:120633/133‑2 (MQ=255)                                                                                                                  gggCGCCGAGGCCGCGCGGGGGCGGCCCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAc               <  4:91907/135‑1 (MQ=255)                                                                                                                  cggCGACGTGGCAGCGCGGGGGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAc               <  2:12093/135‑1 (MQ=255)                                                                                                                  cggCGACGATGCAGCGCTGGGGCGGCGCCGGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAc               <  7:32274/135‑1 (MQ=255)                                                                                                                  cggCGACGAGGCAGCGCGGGGGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGGCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAc               <  7:62911/135‑1 (MQ=255)                                                                                                                  aGGCGACTATGCAGCGCGGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAc               <  7:157729/136‑1 (MQ=255)                                                                                                                                cgcgGGGGCGGCGCCTGCGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATg  <  6:225246/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 GCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCG‑CGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATG  >  ACB122/3902064‑3902324

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