Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,245,145	0	A	C	53.3%	12.0 / 34.4	53	L178R (CTG→CGG)	CJ019_03982	Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (8/13);  new base C (15/9);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAA  >  ACB122/4245020‑4245278                                                                                                                              |                                                                                                                                      cgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                              <  5:25443/135‑1 (MQ=255)           ccGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTAGCTGCCCCACCCGGCGGCCCGGGGGCGGGCCGCGccccc                                                                                                                    >  6:444238/1‑134 (MQ=255)                   tCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCGCGGCGCCCCCCCACtcgt                                                                                                            >  7:278372/1‑135 (MQ=255)                      gcgcATCGGCTATGGGCGGGCCTTTGGCGCGGCGCCGGGCCTGGGCGGGGCCAACGGCTGGCAACTGCGCGGCGGCGCAGACTGGGCTCACCACCCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACtcgtcg                                                                                                          <  6:534896/134‑1 (MQ=255)                        gcATCGGCGACGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGAAAGTGCGCGTCGGCGCAGACTTCGCTCACCACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGa                                                                                                         <  4:443549/133‑1 (MQ=255)                          aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGATCGGCGCCCCCCACCGCGTCGAggc                                                                                                      >  3:377856/1‑133 (MQ=255)                            aGGCGAAGGTCGGTCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGACGCCCCACAACTCGTCGAGGAAca                                                                                                   <  6:181631/135‑1 (MQ=255)                               cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGAAcagca                                                                                                >  7:441462/1‑135 (MQ=255)                                gAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCCACGCGGCGAGGGAccgcag                                                                                               >  3:219531/1‑135 (MQ=255)                                 acgggcGGGCCTTTGCCCCGGCGCCCGGCCAGTTGGGCGGCAACGGCAGGCAACTTCGCGTCGGCGCAGACTTCGCTGACCAACCGGGCGGCCCGGGCGCTAGCGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                              <  4:199345/130‑1 (MQ=255)                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGCCCAACCCGGCGGCCCCGGCGCGCGAGGCGCCCCCCCACCCGGCCAGGCAAAACAGgcgc                                                                                           >  5:118106/1‑135 (MQ=255)                                           cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATGGGCGCCCCACAACTCGGCGAGGCACACCAGGCGCTTGGcca                                                                                    >  3:414473/1‑134 (MQ=255)                                           cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCCACTCGTCGAGGAACCGCAGGCGCGTGGCCt                                                                                    >  3:222464/1‑135 (MQ=255)                                                gcccggcgcccggcCAGTGCGCCGGCAACGGCAGCCAACTGCTCGGCGGCGCAGACTTCGCTGACACACCGGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGCACAGCAGGCGCTTGGCCTGCTCg                                                                               <  8:105485/135‑1 (MQ=255)                                                gCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCCCGCGTCGAGGAACAGCACGCCCTTGGCCGGcgcg                                                                               >  8:188995/1‑135 (MQ=255)                                                 ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACGCGGCGCGGGACACCCGGCGCGTGGCCGGCTCGc                                                                              >  3:195253/1‑135 (MQ=255)                                                     gcgcCCGGCCTGGGGGGCGCCACCGGCCGGCAACGGCTGGGCGGCCCAGATTGCGGTGACCAACCCGGGGGCCCGGGCGCGAGGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                          <  7:246221/135‑1 (MQ=255)                                                      cgcCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCGCGGCGCCCCCCCCCGCGGGGGGGGACAGCCGGCGGGTGGCCGGCGCGCGGcgg                                                                         >  7:499447/1‑133 (MQ=255)                                                        cccGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTAGGGGGCCCCCCACGCGGGGACGAACCCCCCGCGCTTGGGCGTCTCGCAGCCgct                                                                       >  3:190780/1‑133 (MQ=255)                                                         ccGGCCTGTGCGGCGGCAACGGCAGGAAACTGCTCGGCGGCGCAGACTTCGCTGACCAACCGGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATc                                                                      <  6:516822/135‑1 (MQ=255)                                                                    ggcggcAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCcaca                                                                                                                  <  8:384724/80‑1 (MQ=255)                                                                    ggcggcAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCcaca                                                                                                                  >  7:384724/1‑80 (MQ=255)                                                                          aaCGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGCTCGGCGCCCCCCAACTCGTCGCGGAACAGCAGGCCCTTGGCCGGCTCCCTGCCGATCAGCCGGGGCCCGTGCCa                                                     >  5:280196/1‑135 (MQ=255)                                                                                   caACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACTCGTCGAGGAACAGCCGGCGCTTGGCCGGCTCGCTGCCGCGCAGCCGGGGCCGGTGGCCGCGGGCGcc                                            >  3:417115/1‑135 (MQ=255)                                                                                     ccTGCGCGGCGGCGCAGACTTCGCGGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                          <  6:191817/134‑1 (MQ=255)                                                                                       tGCGCGGCGGCGCAGACTGCGCTGCCCACCCCGGCGCCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                        <  8:111618/135‑1 (MQ=255)                                                                                          gcggcgGCGCAGACTGCGCTCCCCACCCCGCCGGCCTGGGCGGTCGGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCg                                     <  5:464049/131‑1 (MQ=255)                                                                                            gtcgGCGAAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAAATCGTCGAGGAACATCAGGCGCTTGTCCTGCTCGCTGCCGTAGAGCCGTGGCCGGTTCCCGCTGGCGCCGGCCTCggg                                   >  6:466709/1‑135 (MQ=255)                                                                                            ggcgGCGCACACTGCGCGCACCACCCGGGCGCCCTGGGGGCGAGGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                   <  3:30913/133‑1 (MQ=255)                                                                                             tcgGCGCAGACTGCGCTGACCACCCCGGCGGCCCGGGCGCTAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                                  <  8:280820/135‑1 (MQ=255)                                                                                             ggggCGAAGACTTCGGCGACCACCCCGGCGGCCCGGGCGCTACCGGGGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                                  <  8:85385/133‑1 (MQ=255)                                                                                                 cgcAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACACCTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCCTCGGGGCt                                >  7:436865/1‑133 (MQ=255)                                                                                                 cgcAGACTGCGCTGCCCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGt                              <  4:146859/135‑1 (MQ=255)                                                                                                      aCTTGGCTGACCACCCCGGGGCCCTGGGGGCGATGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                         <  8:35571/135‑1 (MQ=255)                                                                                                      aCTGCGCTCACAAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                         <  4:331069/135‑1 (MQ=255)                                                                                                        ttCGCTGACCACCCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCAt                       <  4:417115/135‑1 (MQ=255)                                                                                                           gCTCACCACCCGGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCt                    <  6:467966/135‑1 (MQ=255)                                                                                                             tGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGCACAGCCGGCGCGTGGCCTGCTCGCTGCCGCTCAGCCGTGGCCGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCCGGCTgg                  >  8:148268/1‑135 (MQ=255)                                                                                                              gACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCACAACACGTCGAGGAACACCAGGCGCGTTGCCTGCTCGCGGCCGCTAAGCCGTGGCAGGTGCCCCCAGGCGCCCGCCAACGGGCGGTAGCCCAAGCTGgt                 >  4:217422/1‑135 (MQ=255)                                                                                                              gACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCGTCTCGCTGCCGATCAGCCCTGGCAGGTGCCCGCTGGCGCCGGCATCGGGGCTGTAGCCCAAGCTGgt                 >  6:238849/1‑135 (MQ=255)                                                                                                              gACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACACCCCGCGCTTGGCCTGCTCGCTGCCGATCAGCCCTGGCAGGTGCCAGCTGGCGCCGGCAACGGGGCTGTTGCCCATGCTGgt                 >  7:356256/1‑135 (MQ=255)                                                                                                                cccccccGGGGGCCGGGGGGGGTGCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  7:466309/131‑1 (MQ=255)                                                                                                                ccACCCGGGGGGCCCGGGCGCTAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  5:530193/135‑1 (MQ=255)                                                                                                                ccACCCGGGCGCCCCGGGGGGGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  3:217422/135‑1 (MQ=255)                                                                                                                ccACCCCGGGGGCCTGGGCGGGATGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  3:162002/135‑1 (MQ=255)                                                                                                                ccACCCCGGCGGCCTGGGGTCTATCGGCGCCCCACAACTTGTCTAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCTGGTGCCAGCTGGCTCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  5:114368/135‑1 (MQ=255)                                                                                                                ccACCCCGGCGGCCAGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  2:280981/135‑1 (MQ=255)                                                                                                                ccACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  4:219531/135‑1 (MQ=255)                                                                                                                caaccccGGCGGCGTGGGGGCTCGCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt               <  5:345826/131‑1 (MQ=255)                                                                                                                            ccAGGGCGCGATCGGCGCCCCCAAACACCGCGAGGCACATCCGGGGCTTGGCCGGCTCGCTGCCCCGCCGCCGTGGCCGGTGCCCCCGTGCCCCCGCATCGGGGGGGTAGCCCCTGCTGGGGGTGCCgggcgtgg   >  2:346988/1‑129 (MQ=255)                                                                                                                            ccAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                        >  7:449091/1‑98 (MQ=255)                                                                                                                            ccAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                        <  8:449091/98‑1 (MQ=255)                                                                                                                             cAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGaa  <  8:397017/135‑1 (MQ=255)                                                                                                                              |                                                                                                                                      CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAA  >  ACB122/4245020‑4245278

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