Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,245,145 A→C 45.8% 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,1450AC45.8% 8.6 / 24.3 32L178R (CTG→CGG) CJ019_03982Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (4/9);  new base C (10/4);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 5.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

GCGCCGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTT  >  ACB122/4245013‑4245275
                                                                                                                                    |                                                                                                                                  
gcgcCGTCGCGGAGCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAgg                                                                                                                                  <  3:72078/135‑1 (MQ=255)
   ccGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGGCAGTTCGGCGGCCACGGCAGGCCACTGCTCGTCGGCGCAGACGTCGCTGCCCCACCCGGCGGCCCGGgcg                                                                                                                               >  7:2613/1‑135 (MQ=255)
       cgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCGGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCCGACTTCGCTGCCCCCCCCGGCGGCCCCGGcgcgcg                                                                                                                           >  8:279623/1‑133 (MQ=255)
       cgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                           <  3:137189/135‑1 (MQ=255)
       cgcgTATCAGCCGCTTGGGCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                           <  7:29887/135‑1 (MQ=255)
         cgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCg                                                                                                                         >  8:371030/1‑135 (MQ=255)
                   gCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCGGCCCAACCCGGCGGCCCGGGCGCGCTCGGCGccccccc                                                                                                               >  5:261836/1‑132 (MQ=255)
                         tcggcgcATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGACAGTTCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACAAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcg                                                                                                          <  2:57953/131‑1 (MQ=255)
                             cgcATAGGCGAAGGGGGGGCCATTGCCCAGGCGCCCGGCCTGTTCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAg                                                                                                     <  2:306977/135‑1 (MQ=255)
                              gcATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTTGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACACGTCGac                                                                                                     >  7:198253/1‑133 (MQ=255)
                                aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGCTCGCCGGCCCCCACCACGACGCGGac                                                                                                  >  6:358685/1‑134 (MQ=255)
                                  aGGCGAAGGGCGGGCAATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGAAACTGCGCGGCGGCGAAGACTTCGCTGACCAACCCGGCGGCCTGGGCGCTATCGGCTCCCCACAACTCGTCGAGGAAca                                                                                                <  8:306872/135‑1 (MQ=255)
                                   ggCGAAGGGCGGGCAATTGGCCAGGCGCCCGGCAAGTGCGGCGGCAACGGCAGGCAACGGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAAca                                                                                                <  8:258086/134‑1 (MQ=255)
                                     cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGGCGAGGAAcagca                                                                                             >  1:38391/1‑135 (MQ=255)
                                     cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACGCGTCGAGGCAcagca                                                                                             >  8:263459/1‑135 (MQ=255)
                                          gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCCCCCCCCCACACGTCGCCGGACACCAGgggc                                                                                        >  3:137974/1‑135 (MQ=255)
                                          gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCCCTCGGCGCCCCACCCCACGTGGAGGCACCACCGgcgc                                                                                        >  8:364799/1‑135 (MQ=255)
                                             gggCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTAGGCGCCCCCCCACTCGGCGAGGCACAGCAGGCGCTTg                                                                                     >  4:187856/1‑135 (MQ=255)
                                                      gCCCGGCGCCGGGCCAGTTCGGCGGCACCGGCTGGCAAGTGCTCGCCGGCGCAGACTGCGGTGACAACCCCGGCGGCCAGGGCGCTAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCg                                                                            <  2:217094/135‑1 (MQ=255)
                                                           gcgcCCGGCCTGTGCGGCGGAAAGGGACGGCACGGGCTCGCCGGCGCAGAATCCCCTCACCACCCCGCCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                       <  2:210575/135‑1 (MQ=255)
                                                           gcgcCCGGCCAGTGCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                       <  4:154554/135‑1 (MQ=255)
                                                              cccGGCCAGTTCGGCTGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                           <  5:243090/80‑1 (MQ=255)
                                                              cccGGCCAGTTCGGCTGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                           >  6:243090/1‑80 (MQ=255)
                                                                                    gcaggcaACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGATCGGCGCCCCCCAACACGTGCAGGAACCCCCGGCGCTTGGCCTGCGCGCTGCCCAGACGCCGGGGCAGGTGCCCGcgg                                              >  1:268799/1‑133 (MQ=255)
                                                                                             tccTCGCCGGCGCAGTCTGCGCTCACCACCCGGGCGGCCCGGGCGGGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                     <  2:280114/133‑1 (MQ=255)
                                                                                              gggcggcgGCGCAGACTGCGCTGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCAt                                    <  4:137974/129‑1 (MQ=255)
                                                                                                tcggcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCg                                  <  6:261836/135‑1 (MQ=255)
                                                                                                gcggcgGCGCCGACTTCGCTCACCACCCCGGCGGCCCGGGCGCGAGCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCg                                  <  7:279623/131‑1 (MQ=255)
                                                                                                     ggCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGGCCCCCAACTAGGCGAGGAACCGCAGGCGCTTTGGCGGCGCCCTGGCGCAGAGCCGGGGCAGGTGGCACCTGGGGCCGGCAAAGGGGcg                             >  2:272600/1‑134 (MQ=255)
                                                                                                                      ccACCCCGGCGGCCGGGGCGCGAGCGGCGCCCCACACCTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGAGCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt            <  7:216163/135‑1 (MQ=255)
                                                                                                                       caccccGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTa           <  3:179342/132‑1 (MQ=255)
                                                                                                                                ggcCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCtt  <  2:130154/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GCGCCGTCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTT  >  ACB122/4245013‑4245275

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