Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,880,363	A→C	45.9%	V96G (GTG→GGG)	CJ019_04574 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,880,363	0	A	C	45.9%	18.5 / 5.5	33	V96G (GTG→GGG)	CJ019_04574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/12);  new base C (13/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCG  >  ACB122/4880231‑4880496                                                                                                                                     |                                                                                                                                        gCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCCCCGGcacc                                                                                                                                       >  8:368412/1‑135 (MQ=255)  cccGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGcacca                                                                                                                                      <  5:351037/135‑1 (MQ=255)   ccGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGcaccac                                                                                                                                     <  6:77578/135‑1 (MQ=255)   ccGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGcacc                                                                                                                                       <  2:154278/133‑1 (MQ=255)    cGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCCCCCCCCCCCGCCCCCGCcccccc                                                                                                                                    >  7:136105/1‑135 (MQ=255)      cAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCgg                                                                                                                                  <  8:174259/135‑1 (MQ=255)        gCCCCAACTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGcc                                                                                                                                <  8:47722/135‑1 (MQ=255)               ctcctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCCGCCCCACCCGCCCCGCA‑cc                                                                                                                        >  1:263375/1‑135 (MQ=255)                         gCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCCGCCGCTCC‑CCCATCCC‑ggt                                                                                                              >  3:279632/1‑132 (MQ=255)                           aCTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCaa                                                                                                           <  4:70272/135‑1 (MQ=255)                                          cTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGCATCAGCACCCCACCCGGCGCCGGCGCCCCCGGCCCCGCC‑CCCCACCC‑GGGCAACGCGGCCAAAAAAgc                                                                                            >  8:355309/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCCCCCCCCCCCGCGCCGGCCCCCCCGGCCGCGCC‑CCCATCCC‑CGTCAACGCGGCCAAAAAAgccg                                                                                          >  8:17054/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCGGCCCCACCCGCCGCGAA‑CCCCTCCC‑GGTCAACGCCGCATAAAAAgccg                                                                                          >  2:160427/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCCGCACCCCCCGCCCCGCA‑ACCATCCC‑GGTCAACGCGGCAAAAAAAcccg                                                                                          >  6:272043/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCCCCGGCCCCCCCCGCCCCGGACCCCACCCC‑GTGCACGCCCGCaaaaaaaccg                                                                                           >  2:170722/1‑132 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCGGCCGCGCA‑CCCCTCCC‑GGTCAACGCGGCAAAAAAAgccg                                                                                          >  1:64414/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCCGCGCCCGCCCCCCCGGCCGCGCC‑CCCCTCCCGGTGCAAC‑CCGCAAAAAAAcccc                                                                                          >  2:115765/1‑134 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCAACCGGCGCCGGCCCCACCGGCCGCGCA‑CACCTCCC‑GGTCAACGCGGCAAAAAAAgccg                                                                                          >  1:70177/1‑135 (MQ=255)                                            cgtcgATATCCGACGTATCGCCGCTGATCACCACGGCACCCACCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCACCACCGGCCCCGAA‑CCCCACCC‑CTTCCACTCGGCAAAAAAAgccg                                                                                          >  6:135403/1‑135 (MQ=255)                                                     ccGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCCGCCCCACCCGCCGCACA‑CCCATCAC‑CTTCAACCCGGCAAAAAAAGCCCGCCgctgca                                                                                 >  4:173166/1‑134 (MQ=255)                                                       gACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCGGCCCCACCGGCCGCGCC‑CCCCTCCC‑GGTCAACGCGGCCAAAACAGCCCGCCCCtgcggc                                                                               >  3:372523/1‑135 (MQ=255)                                                         cGTATCGCCGCTCATCCCCACGCCACCCAGCACCTGGCCGTCCTGATCACGAATCACCACCCCACCCGGCGCCGGCACCCCCGGCCGCTCA‑CCCAGCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTgcgc                                                                             <  4:209378/135‑1 (MQ=255)                                                          gTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCCCCACCGGCCGCGCC‑CCCCTCCC‑CGTCCACGCGGCAAAAAAAGCCGGCCGCTGCTgcgcc                                                                            >  8:243958/1‑134 (MQ=255)                                                                               gCACCCAGCACCTTGCCGCCCTCATCACGAAGCAGCCCCCCCCCCGGCGCCGGCCCCACGGGCCGCACA‑CCCAGCCC‑TTTCCACGCGGAAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAg                                                       <  7:158206/135‑1 (MQ=255)                                                                                    cAGCACCTTGCCGTCCTGAGCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCtt                                                  <  6:173147/135‑1 (MQ=255)                                                                                        aCCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCCCCGGCCGCACC‑CCCATCCC‑GGTCAACGCCGCAAAAAAAGCCGGCCCCTGCTGCGCCTCCACCGCCAGCCCCCGCGCGCCCATGccc                                              >  3:67591/1‑135 (MQ=255)                                                                                          cTTGCCGTCCTGATCACCAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCAACCC‑TTTCAACGCGGCAAAAAAACCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAg                                            <  5:272043/135‑1 (MQ=255)                                                                                               cGTCCTGATCACTAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTAAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                       <  3:173166/135‑1 (MQ=255)                                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCGGCCGCTCC‑CCCATCCC‑GTTCAACGCCGCAAAAAAAGCCGCCCGCTGCTGCGCCTCCCCCGCCACCCAGCCCGCCCCCATGCCCAgcgcca                                       >  8:305567/1‑133 (MQ=255)                                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCCCCCGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCAGCCAGCGCCCGCAACCGCGAGCCCATGGCCAGGGcc                                        >  8:187871/1‑133 (MQ=255)                                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCC‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGCGCCCTTGCCCAGGGCaa                                       >  1:36931/1‑135 (MQ=255)                                                                                               cGTCCTGATCACGAAGCAGCACCCCCCCCGGCGCCGGCACCCCCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                       <  5:149642/135‑1 (MQ=255)                                                                                                                    ccccacccGGCTCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCTCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAg                  <  2:75837/135‑1 (MQ=255)                                                                                                                     cccacccGGCGCCTGCACCACCGGCCGATCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCTCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTTCCAGt                 <  2:70177/135‑1 (MQ=255)                                                                                                                     cccacccGGCGCCGGCCCCACCGGCCGCTCA‑CCCATCCC‑GTTCCACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                 <  7:266556/135‑1 (MQ=255)                                                                                                                     cccacccGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                 <  8:136105/135‑1 (MQ=255)                                                                                                                     cccacccGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                 <  2:64414/135‑1 (MQ=255)                                                                                                                          ccGGCGCCGGCACCACCGGCCGCTGC‑CCCATCCC‑CGTCAACCCGGCAAAAAAACCCCGCCCCGGCTGCGCATACCCCGCCCCCCACCCCCCGCCCCTGCCCCCGGCCAAAACACCCCACGCGTTGCCAGGTGCCa            >  4:66075/1‑135 (MQ=255)                                                                                                                                    caccacCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGCCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  6:292987/135‑1 (MQ=255)                                                                                                                                    caccacCGCCCGCGCA‑CCCAACCC‑GTGCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  1:213557/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCG  >  ACB122/4880231‑4880496

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