Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,245,154 T→G 49.8% D175A (GAT→GCT)  CJ019_03982 ← Crotonyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,245,1540TG49.8% 4.1 / 19.0 37D175A (GAT→GCT) CJ019_03982Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (9/7);  new base G (3/16);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 2.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAAGCGCGCCGAG  >  KT2440/4245020‑4245288
                                                                                                                                       |                                                                                                                                       
cgcgTATCAGCCGCTGGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                                         <  4:56020/135‑1 (MQ=255)
                                aaGGGCGGGCCTTTGGCCAGGCGCCCGGCCTGTGCGGCGGCAA‑CGGAAGGCAACTGCTCGGCGGCGCATAGTGCGCTGCCCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                                         <  8:289255/135‑1 (MQ=255)
                                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGCTCGGCGCCCCACCACTCGTCGAGGACCACCAGgcgc                                                                                                      >  4:16545/1‑135 (MQ=255)
                                   gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCCCGCGGCGCCCCCCACCACCTCGGGGGACAGCCGgcgc                                                                                                      >  3:234694/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGCCCAACCCGGCGGCCCGGGCGCGCGCGGCCCCCCCCACCCCGTCGGGGGACCCCAGGCGC‑GTGGGCt                                                                                              >  1:78225/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGGGCGCTGGGCGCGCCCCCACTCGGGGCGGGACCGCCGGCGCGTTGGcc                                                                                               >  1:144033/1‑135 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCGCGGGGCCCCCCCACTCGGGGCGGGACCCCACGCGC‑GTGGccg                                                                                              >  7:299693/1‑134 (MQ=255)
                                          cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCCACTCGTCGAGGAACACCAGGCGC‑GTGGCCt                                                                                              >  3:101806/1‑135 (MQ=255)
                                               ccctggCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCGCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCg                                                                                         <  6:214961/131‑1 (MQ=255)
                                               ccctggCCCCCGGCCCGGTGGGCGCCAA‑CGGCAGGCAACGGCGCGGCGGCCCAGACTTCGCGGACCAACCCGGCGCCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCg                                                                                         <  6:166537/131‑1 (MQ=255)
                                                ccAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCCACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGGCTGCTc                                                                                          >  8:233011/1‑133 (MQ=255)
                                                    gcgcCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTGCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGcc                                                                                    <  8:173311/135‑1 (MQ=255)
                                                    gcgcCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACGGCTCGGCGGCGCAGACGTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGcc                                                                                    <  7:190784/135‑1 (MQ=255)
                                                    gcgcCCGGCCAGTGGGGCGGCAC‑CGCCAGGCAACTGCGCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCTAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGcc                                                                                    <  1:173834/135‑1 (MQ=255)
                                                     cgcCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACGCGGGGGGGGACCGCAGGCGC‑GTGGCCGGCGCGCGGc                                                                                     >  5:282124/1‑133 (MQ=255)
                                                       cccGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGTCGCGGAACCGCCGGCGC‑GTGGCCGGCGCGCTGCCgct                                                                                 >  4:162940/1‑133 (MQ=255)
                                                       cccGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACGCCTCGCGGAACCCCCGGCGC‑GTGGCCTGCTCGCGGCCgc                                                                                  >  4:58541/1‑133 (MQ=255)
                                                         cGGCCCGTGGGGCGGCAC‑CGGCCGGCACCTGCGCGGCGGCGCAGACTGCGCTCACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCa                                                                               <  5:49392/135‑1 (MQ=255)
                                                          ggCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGCTCAg                                                                              >  8:250663/1‑135 (MQ=255)
                                                              tgTTCGGCGGCAA‑CGGCAGGCAACTGCTCGGCGGCGCCGACTGCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGt                                                                          <  6:51737/134‑1 (MQ=255)
                                                              cgtggggcgccAA‑GGGCCGGCAACTGCGCGCGGGGGCACACTTCGCTCCCCAGCCGGGGGGCCGGGCGGGGAGGGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGt                                                                          <  8:258004/130‑1 (MQ=255)
                                                              agctgggcggcAC‑CGGCAGGCAACTGCTCGTCGGCGCAGACTGCGCTGACCCCCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGt                                                                          <  2:239229/130‑1 (MQ=255)
                                                               tgtgggcgggAACCGGCCGGCACCTGCGCGGCGGCGCTGCGGTGGCAGCCCAACCGGCCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGt                                                                          <  2:65052/131‑1 (MQ=255)
                                                                        cAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGCTCAGCCGTGGCAGGTg                                                                  >  8:222217/1‑133 (MQ=255)
                                                                              cccggaaaCGGCTCGGCGCCGCAGACTTCGCCGACCAACCCGGCCCCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTg                                                           <  7:35990/129‑1 (MQ=255)
                                                                                            gtcgGCTCTGAGTTCGCTTACCACCCCGGCGGCCCGGGGGCGAGGGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                             <  2:242198/135‑1 (MQ=255)
                                                                                            ggcgGCGCAGACTCCGCTCACCAACCCGCCGGCCCGGGCGCGAGCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                             <  7:292378/133‑1 (MQ=255)
                                                                                            ggcgGCCCAGACGTCGCTGCCCCACCCGGCGGCGCGGGCGCGAGCGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                                             <  1:105959/133‑1 (MQ=255)
                                                                                                 cgcAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGt                                        <  6:215069/135‑1 (MQ=255)
                                                                                                 cgaaGACTACCCTGACCACCCCGGCGGCCAGGGCGGTATCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGt                                        <  2:268876/132‑1 (MQ=255)
                                                                                                    aGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCACGGCGCCCCACAACACGTCGCCGCACCGCCCGCGC‑ATTGCCCGCTAGCAGCCGATACGCCGTCGGACGTTGCCGCTGGCGCCGGCCACGGGGcgttgtg                                     >  6:98485/1‑129 (MQ=255)
                                                                                                       gtgcGCTGACACCCCGGGCGCCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                                   <  2:262016/131‑1 (MQ=255)
                                                                                                             tGACCAACCCGGCGGCCACGGCGCGGAACGCGCCCCCCAACTCGTGGAGGCACACCCGGCGC‑TTGGCCAGCGCGCGGGCCCAAAGCCCTGGCAGGTGCCCCCGGGCGCCCGCCGCGGGGCGGTATCCCagtcggg                            >  6:284955/1‑129 (MQ=255)
                                                                                                                cccccccGCCGGCGGGGGCGTGTTGGGCGCCCCACAAGTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt                         <  1:91618/131‑1 (MQ=255)
                                                                                                                ccACCCCGGCGGCCCGGGCGCTAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt                         <  5:82714/135‑1 (MQ=255)
                                                                                                                ccACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑GTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt                         <  3:286647/135‑1 (MQ=255)
                                                                                                                ccAACCCGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt                         <  1:183423/135‑1 (MQ=255)
                                                                                                                  accccGGCGGCCCGGGGCGGAGCGGCGCCCCACAACTCGTCGAGGAACCGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAg                       <  6:34893/133‑1 (MQ=255)
                                                                                                                                       tCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTAGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAAGCGCGCCgcg  >  8:219509/1‑133 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAA‑CGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGC‑TTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTTGAAGCGCGCCGAG  >  KT2440/4245020‑4245288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑