Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,245,145 A→C 49.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
*ACB1224,245,1450AC49.3% ‑6.5 / 9.2 34L178R (CTG→CGG) CJ019_03982Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base C (10/5);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCC  >  ACB122/4245036‑4245269
                                                                                                             |                                                                                                                            
ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACACCtc                                                                                                     >  4:145334/1‑135 (MQ=255)
    tGCGCTTCGGCGACGGGCGGGCCATGGGCCTGGCGCCCGGCATGTGCGGCGGCACCGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACtcgtcg                                                                                                 <  2:46350/135‑1 (MQ=255)
         aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCACCCCGGCCAGGa                                                                                             >  4:259831/1‑134 (MQ=255)
         aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGATCGGCGCCCCCCCACCCGTCGAggg                                                                                             >  2:256211/1‑133 (MQ=255)
         aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACCACTCGTCGAGGa                                                                                             >  2:67389/1‑134 (MQ=255)
         aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGaa                                                                                            >  3:196236/1‑135 (MQ=255)
           gggCGAGGGGGGGGCCTTTGGCCTGGCGCCCGGCCTGTGCGGCGGCAACGGCTGGCAAGGGCTCGGGGGCGAAGACTGGGCTGACCACCCCGGCGGCCGGGGCGGGAGGGGCGCCCCACAACTCGTCGAGGAAca                                                                                          <  6:214670/134‑1 (MQ=255)
           aGGCGATGGGCGGGCAATTGGCCTGGCGCCCGGCCTGTGCGGCGGCAACGGCAGGAAACTGCTCGTGGGCGCATACTTCGCTGCCCAACCCGGCGGCCTGGGCGCTTTCGGCGCCCCACAACTCGTCGAGGAAca                                                                                          <  6:72323/135‑1 (MQ=255)
           aGGCGAAGGGCGGGCCATGGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACGGCGCGTCGGCGCAGACTGCGCTGCCCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAAca                                                                                          <  6:302202/135‑1 (MQ=255)
                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATAGGCGGCCCCCCACGCGGCGGCGAACACCCGgcgc                                                                                  >  3:46327/1‑135 (MQ=255)
                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCCCCACCCGGCGCGGACcacccgcc                                                                                    >  4:304614/1‑127 (MQ=255)
                     cGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCACGGCGCGCTCGGCGCCCCCCCACTCGTCGACGGACAGCACGCGCTTg                                                                               >  7:117982/1‑136 (MQ=255)
                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAGGCACAGCACGCGCGTGGccgg                                                                          >  7:44550/1‑134 (MQ=255)
                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGAACCGCAGGCGCTTGGCCt                                                                           >  4:315015/1‑135 (MQ=255)
                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACACCTCGCGGGACACCACGCGCtgtggca                                                                           >  2:140660/1‑129 (MQ=255)
                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCCGTCGGCGGCCCCCCACACGTAGGGGGACAGCACGGGCTTgggca                                                                           >  6:339214/1‑132 (MQ=255)
                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCCACTCGTCGAGGCACAGCAGGCGCTTGGCCt                                                                           >  2:254264/1‑135 (MQ=255)
                                       cccGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGCCTTCGCTGCCCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCCCGCGTCGCGGCACAGCCGGCGCGTGGCCGGCGCGCGGCCgct                                                              >  3:110149/1‑133 (MQ=255)
                                         cGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACGCGTCGCGGCACAGCCGGCGCGTGGGCGGCTCGCTGCCGCTCac                                                           >  7:225937/1‑135 (MQ=255)
                                              tgtggggcggcACGGGCCGCCACCGGCGCGGCGGCGCAGACTGCGCTCCCCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGt                                                       <  3:204580/130‑1 (MQ=255)
                                                                    aCTGCGCGGCGGCGCAGACTTCGCTCCCCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGGGCCAGCTGGCGCCgg                                 <  4:110149/135‑1 (MQ=255)
                                                                         tcgtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGTCGCGGAACAGCCGGCGCGTGGCCTGCTCGCGGCCGATCCGCCGGGGCAGGTGGCCGCGGGCGCCGGCCTCg                            >  1:236174/1‑135 (MQ=255)
                                                                           ggcgGCGCAGACTGCCCTCACCCCCCCGGCGGCCGGGCCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                          <  3:228465/133‑1 (MQ=255)
                                                                           ggcgGCGCACACTGCCCGCACCCCCCCGGCGCCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                          <  5:88507/133‑1 (MQ=255)
                                                                           ggcgGCCCAGAGTGCGCCCACCCCCCCGGCGCCCCGGGCGCGAGCGCCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                          <  1:256211/133‑1 (MQ=255)
                                                                            ccgGCGCTTACTGCGGTGCCCACCCGGGGGCCCTGGGCGGGTGGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                         <  3:127006/134‑1 (MQ=255)
                                                                               gcgcAGACTTCGCTCACAAACCGGGCGCCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGt                     <  8:42496/136‑1 (MQ=255)
                                                                                   aGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATAGGCGCCCCACCACTCGGCGAGGAACAGCCGGCGCTTGGCCTGCTAGCTGGCCATCCGCCGGGGCAGGTGCCCGCTGGCGCCGGCCTCGGGGCGGTAGc                  >  4:173201/1‑135 (MQ=255)
                                                                                     aCTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCa               <  8:117982/136‑1 (MQ=255)
                                                                                          gCTGACCACCCCGGCGGCCAGGGCGGTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGc            <  6:153919/134‑1 (MQ=255)
                                                                                              accaccCCGGGGGCCCGGGGGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt      <  8:148668/136‑1 (MQ=255)
                                                                                               ccACCCCGGGGCCCTGGGCGCGTTCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt      <  6:287190/135‑1 (MQ=255)
                                                                                               ccAACCCGGCGGCCAGGGGGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt      <  3:112446/135‑1 (MQ=255)
                                                                                                  aCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGcc  >  8:127317/1‑136 (MQ=255)
                                                                                                             |                                                                                                                            
GGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCC  >  ACB122/4245036‑4245269

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: