Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,245,145	A→C	34.3%	L178R (CTG→CGG)	CJ019_03982 ←	Crotonyl‑CoA hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,245,145	0	A	C	34.3%	9.5 / 4.7	33	L178R (CTG→CGG)	CJ019_03982	Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (7/13);  new base C (4/8);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTG  >  ACB122/4245020‑4245276                                                                                                                              |                                                                                                                                    cgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGGCTTCGCTGCCCAACCCGGCGGCCAGGGCGCGCTc                                                                                                                           >  3:315747/1‑136 (MQ=255)                 ggTCGGCGCATCGGCGAAGGGCGGGCCATTGGCCAGGCGCCCGGCCAGTGCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcg                                                                                                           <  4:100064/136‑1 (MQ=255)                     tGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcgtcg                                                                                                        <  4:355540/135‑1 (MQ=255)                        gcATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGCTCGGCCCCCCCCCACCCCTCGAggg                                                                                                    >  3:239049/1‑135 (MQ=255)                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCACCACTCGTCGAGGCACATCAGgcgc                                                                                         >  2:72109/1‑135 (MQ=255)                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCAACAAGTCGAGGAACAGCAGGCCCt                                                                                        >  7:10697/1‑136 (MQ=255)                                           cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGGGCGCTCGGCGCCCCCCAACTCGTGCAGGAACAGCCGGCGCGTGGccgg                                                                                 >  3:107084/1‑134 (MQ=255)                                           cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGAACCGCCGGCGCGTGGCCTg                                                                                 >  7:245099/1‑136 (MQ=255)                                             ttGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCt                                                                               <  8:131417/136‑1 (MQ=255)                                                    ggCGCCCGGCCAGTGCGGCGGCAACGGCAGGCAACTGCGCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                        <  8:22202/136‑1 (MQ=255)                                                     gcgcCCGGCCAGTTCGGCGGCAACGCCAGGCAACTGCACGTCGGCGCAGACTTCGCGGACCACCCCGCCGCCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                        <  1:72109/135‑1 (MQ=255)                                                             ccAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACCACTCGGCGAGGAACAGCAGGCGCTTGGCCTGCGCGCGGCCGAGccgccg                                                               >  3:308448/1‑136 (MQ=255)                                                              cAGTGCGGCGCCAACGGCCGGCAACTGCGCGTCGGCGCAGACTTCGCGGACCACCCCGGCGCCCAGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGt                                                              <  4:280712/136‑1 (MQ=255)                                                              cAGGTCGGCGGCAACGGCAGGCAACTGCGCGTCGGCGCAGACTTCGCGCCCCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGt                                                              <  8:126308/136‑1 (MQ=255)                                                                         cacgggcaggcaACTGCTCGCCGGGGAAGACTCGCCTCCCCAACCCCGCGGCCCGGGCGGGAGCGCGGCCCCACAACTCGGCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCa                                                   <  8:111774/132‑1 (MQ=255)                                                                         cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCAACTCGTCGAGGCACAGCAGGCGCTTGGCCTGCTCGCTGCCGCTCAGCCGTGGCAGGTGcc                                                    >  6:185670/1‑135 (MQ=255)                                                                                     aCTGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                        <  2:322025/135‑1 (MQ=255)                                                                                     aCTGCTCGGCGGCGCAGACTGCGCTGACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                        <  2:316818/135‑1 (MQ=255)                                                                                      cTGCTCGGCGGCGCAGACTGCGCAGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCGTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                      <  4:77884/136‑1 (MQ=255)                                                                                        gcgcgtcgGCGCAGACTTCGCTGACCACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCAt                                     <  2:268715/132‑1 (MQ=255)                                                                                          tcgtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCCCAACTCGTCGAGGAACAGCCGGGGCTTGGGCTGCTCGCTGCCGCTCAGCCGTGGCAGGTGGCAGCTGGCGCCGGCATCgg                                  >  7:310391/1‑136 (MQ=255)                                                                                           cgtcgGCGCAGACTTCGCTGACAAACCCGGCGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                 <  3:170588/136‑1 (MQ=255)                                                                                                 cgcAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGt                            <  2:32403/135‑1 (MQ=255)                                                                                                    aGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACCACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCGGTAGcc                        >  7:334854/1‑136 (MQ=255)                                                                                                               cccAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  4:132798/135‑1 (MQ=255)                                                                                                               acccccccGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCGGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  4:13878/131‑1 (MQ=255)                                                                                                               acccacccGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  4:352857/136‑1 (MQ=255)                                                                                                               accaccCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  8:74461/136‑1 (MQ=255)                                                                                                               aCCAACCGGGCGCGCAGGTCGCGATCGGCGCCCCACAACTCGTCGAGGACCAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGGCGTGTCAGTTTCGATCTTGCTGCGTATTCGTGGCTGTAGCCCATGCTGgtgt             <  4:316785/136‑1 (MQ=255)                                                                                                                ccACCCGGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  5:78576/135‑1 (MQ=255)                                                                                                                ccACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt             <  1:299446/135‑1 (MQ=255)                                                                                                                     ccGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCg        >  1:347329/1‑135 (MQ=255)                                                                                                                           gcCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTg  <  2:59151/135‑1 (MQ=255)                                                                                                                              |                                                                                                                                    CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTG  >  ACB122/4245020‑4245276

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