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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,245,1450AC43.5% 14.9 / 11.6 51L178R (CTG→CGG) CJ019_03982Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (6/17);  new base C (16/7);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 7.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGT  >  ACB122/4245020‑4245265
                                                                                                                             |                                                                                                                        
cgcgTATCCGCCGCTTGGTCTGCGCAACGGCTAAGGTCGGGCCATGGGCCAGGCGCCCGGCCAGTTCGGCGGAAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                 <  2:240001/135‑1 (MQ=255)
cgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                 <  7:419585/135‑1 (MQ=255)
cgcgTAGCAGCCGCTTGGTCGGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCGCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                 <  4:441438/135‑1 (MQ=255)
          ccGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCca                                                                                                       <  7:373889/135‑1 (MQ=255)
                ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAgg                                                                                                                        <  5:19431/112‑1 (MQ=255)
                ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAgg                                                                                                                        >  6:19431/1‑112 (MQ=255)
                  tCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcg                                                                                                <  3:176702/134‑1 (MQ=255)
                      cgcATCGGCGATGGGCGGGCCATTGGCCAGGCGCCCGGCCAGGTGGGGGGCAACGGCAGGAAACTGGTGGTCGGCGCATACTGCGCGCACCAACCGGGCGGCGAGGGCGCGATGGGCGCCCCACAACTCGTCGAg                                                                                           <  4:228067/135‑1 (MQ=255)
                       gcATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCCACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAg                                                                                           >  4:88113/1‑134 (MQ=255)
                         aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGGCGAGGaa                                                                                        >  2:364660/1‑135 (MQ=255)
                           aGGCGATGGTGGGGCCTTTGGCCTGGCGCCCGGCCAGTTGGGCGGCAACGGCTGGCAACTGCTCGTGGGCGCATACTTCGCTGACCAACCCGGGGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAca                                                                                      <  3:217168/135‑1 (MQ=255)
                                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACCACTCCTCGCCGAACACCCGgcgc                                                                              >  7:130940/1‑135 (MQ=255)
                                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCCACACGTCGAGGAACAGCCGgcgc                                                                              >  3:185795/1‑135 (MQ=255)
                                     ggggCCATGGCCCAGGCGCCCGGCCAGTGCGGCGCCAACGCCAGGCAACTGCGCGGCGGCGCAGACTTCGCTCACCACCCGGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCtt                                                                            <  3:139899/134‑1 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGCCCACCCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACGCGTCGCGGGACCGCCGGCGCGTgggcg                                                                       >  7:205758/1‑132 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCCCCCCGGCGGCCCGGGCGCGCGGGGGGGCCCCCCACGCGTCGGGGGACCGCCGGCGCGTgggcg                                                                       >  5:314110/1‑132 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGGGCGCTCGGGGCCCCCCCACCCGGGGGGGGACACCCGGCGCTgggggt                                                                       >  1:211631/1‑132 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGTCGAGGAACAGCAGGCGCTTGGCCt                                                                       >  7:143290/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCCCTCGGGGCCCCCCCACACGTCGAGGCACCCCACGCCCGTGGCCt                                                                       >  7:133483/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCt                                                                       >  8:348564/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTAGGCGGCAACGGCAGGCAACAGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGCCGCCCCACAACTCCTCGAGGAACACCAGGCGCATGGcca                                                                       >  3:134025/1‑134 (MQ=255)
                                             gggCCTGGCGCCCGGCCAGTGCGGCGGCAACGGCTGGCAACTGGTGGTGGGGGCTGACTTGGCTGACCACCCGGGCGGCATGGGCGGGATCGGGGCCCCACAACTCGGCGAGGAACAGCAGGCGCTTGGCCTGCt                                                                    <  7:136865/134‑1 (MQ=255)
                                                ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATAGGCGCCCCCCCACACGGCGGGGCACACCACGCGCCTGGCCGGcgggg                                                                 >  4:139899/1‑131 (MQ=255)
                                                    gcgcCGGGCCTGTGCGGCGGCAACGGCAGGCACCTGCGCGGCGGCGCAGACTTCGCTCACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                             <  4:284050/135‑1 (MQ=255)
                                                     cgcCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGGCCCCCCCCACGTCGCGGGACACCCGGCCCGTGGGCGGCTCGCGGCCg                                                            >  6:278816/1‑135 (MQ=255)
                                                            ccAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGGCGAGGCACAGCACGCGCGTTGCCTGCTCGCTGCCGCTGAGcc                                                     >  8:137439/1‑135 (MQ=255)
                                                              tgTTGGGCGGCAACGGCAGGCACCGGCGCGGCGGCGCAGACTTCGCTCACCACCCCGGCGCCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGt                                                   <  3:457201/134‑1 (MQ=255)
                                                                        cAACGGCAGGCACCTGCGCGGCGGCGCAGACTTCGCTCACCAACCCGGCGGCCAGGGCGCGATGGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGcc                                         <  7:290092/135‑1 (MQ=255)
                                                                             gcaggcaACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCCGGGCCCCCCCACACGTCGGGGGACAGCCCGCGCTTGGCCGGCTAGCTGCCGATCAGCCCGGGCCGGTGGCAGcgg                                    >  2:243948/1‑133 (MQ=255)
                                                                                    acggcgcgtcgGCGCAGCCTGCGCGGACCCTCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGCGGAACAGCAGGCGCTTGTCCTGCTCGCTGCCGAGCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                             <  8:351151/129‑1 (MQ=255)
                                                                                    aCTGCGCGTCGGCGCAGACTTCGCCGACCACCCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                             <  8:296343/135‑1 (MQ=255)
                                                                                      tGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                           <  2:10060/135‑1 (MQ=255)
                                                                                      ggcgcgtcgGCGCAGACTTCGCTGACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                           <  2:341540/131‑1 (MQ=255)
                                                                                      ggCTCGGCGGCGCAGCCTGCGCGCCCCACCCCGGCGCCCCGGGCGCGAGCGCCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                           <  5:113319/134‑1 (MQ=255)
                                                                                           gtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACACGGCGAGGCACAGCAGGCGCGTGGCCTGCTCGCTGCCGCTCCGCCGGGGCAGGTGCCCGCTGGCGCCGGCCTCggg                      >  3:134774/1‑135 (MQ=255)
                                                                                           gtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACACGGCGCGGCACAGCAGGCGCGTGGGCGGCTCGGTGCCGCACACCCGGGGCAGGTGCCCGCTGGCGCCGGCCTCggg                      >  6:436231/1‑135 (MQ=255)
                                                                                            tcgGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                     <  8:133483/135‑1 (MQ=255)
                                                                                            gcgGCGCAGACTTCGCCCACCACCCGGGCGGCCTGGGCTGGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCGCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                     <  3:37744/134‑1 (MQ=255)
                                                                                            gcgGCGAAGACTTCGCTGACCGACCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                     <  4:27086/134‑1 (MQ=255)
                                                                                            ccgGGGCAGACTTCGCGAACAACCCCGGCGGCCAGGGCGGGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                      <  1:53617/133‑1 (MQ=255)
                                                                                            caggCGCACACTTCGCTGACCACCCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                     <  8:136828/133‑1 (MQ=255)
                                                                                             cgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGGTCGGCGCCCCCCAACTCGTGGAGGTACAGCAGGGGCGTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCCGCTGGCGCCGGCCTCGGGGc                    >  5:85025/1‑135 (MQ=255)
                                                                                                 gcAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTa                <  7:361861/135‑1 (MQ=255)
                                                                                                          gctccccACCCCGGCGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAGTAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGc        <  3:207728/129‑1 (MQ=255)
                                                                                                            tGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACCACTCGTCGAGGAACAGCAGGCGCGTGGCCTGCTCGCGGCCGAGAAGCCGTGGCAGGTGCCCGCGGGCGCCGGCATCCGGGCTGTAGCCCATGCTgg     >  4:200779/1‑135 (MQ=255)
                                                                                                             gACCAACCCGGCGGCCCGGGCGCGATAGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGCTCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgt    >  8:116117/1‑135 (MQ=255)
                                                                                                             gACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCtctg    >  7:136828/1‑132 (MQ=255)
                                                                                                               ccACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt  <  1:254260/135‑1 (MQ=255)
                                                                                                               ccAACCGGGCGGCCAGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt  <  3:177055/135‑1 (MQ=255)
                                                                                                               ccAACCCGGCGGCCCGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt  <  7:123549/135‑1 (MQ=255)
                                                                                                               cacccccGGCGGCCAGGCCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt  <  6:172473/131‑1 (MQ=255)
                                                                                                                             |                                                                                                                        
CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGT  >  ACB122/4245020‑4245265

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