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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,245,1450AC45.6% 6.7 / 16.2 51L178R (CTG→CGG) CJ019_03982Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base C (17/5);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTT  >  KT2440/4245022‑4245275
                                                                                                                           |                                                                                                                                  
cgTATCAGCCGCTTGGTCTGCGCCTAGGCGAAGGTCGGGCCATTGGCCCGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGGGCAGCCTTCGCTGCCCAACCCGGCGGCCCCGGcgcgcgcgg                                                                                                                        >  7:356650/1‑131 (MQ=255)
        cgtgttGGTCTGCGAATTGGCGATGGTCGGGCCTTTGGCCTGGCGCCCGGCCTGTTCGGGGGAAACGGAAGGAAAGTGATCGTCGGCGCAGACTTCGCTGACAAACCCGGGGGCCCGGGGGCTATCGGCGCCCca                                                                                                                 <  6:134814/131‑1 (MQ=255)
        ccGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCACGGCGCGCTCGGCGccccc                                                                                                                 >  2:305066/1‑134 (MQ=255)
          gCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCCACTGCTCGTCGGCGCAGACTTCGCTGCCCAACCCGGCGGCCCGGGCGCGGTCGGCGCCCcaca                                                                                                               >  3:184989/1‑135 (MQ=255)
               accggcgcATAGGCGAAGGGCGGGCCATGGGCCAGGCGCCCGGCCAGTGCGGCGGCAACGGCAGGCAACGGCTGGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcgt                                                                                                         <  7:261362/132‑1 (MQ=255)
                 cTGCGCATAGGCGAAGGGCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACtcgtcg                                                                                                       <  8:74250/136‑1 (MQ=255)
                       aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACACCGCGAGGCAc                                                                                                 >  7:67247/1‑136 (MQ=255)
                         aGGCGACGGGCGGGCCATTGGCCAGGCGCCCGGCCTGTGCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCAACCGGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAca                                                                                                <  3:117037/135‑1 (MQ=255)
                            cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAGGAAccgcc                                                                                             >  3:326644/1‑134 (MQ=255)
                            cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAcagca                                                                                             >  4:360671/1‑135 (MQ=255)
                             gAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCCCTGACCAACCCGGCGGCCCCGGCGCCCACGGCGCCCCCCACCACGTCGCGGCAccccacg                                                                                           >  7:143212/1‑130 (MQ=255)
                                     ggCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCGCGGCGGCGCAGACTTCGCTGCCCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTgg                                                                                    <  6:343349/135‑1 (MQ=255)
                                        cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACACGTCGAGGGACCCCCGGCGCGTgggcg                                                                                 >  4:381057/1‑132 (MQ=255)
                                        cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCGCGGCGGCGCAGACTGCGCGGACCAACCCGGCGCCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCt                                                                                 <  3:337536/135‑1 (MQ=255)
                                        cATTGGCCAGGCGCCCGGCCAGTGCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCt                                                                                 <  1:89889/135‑1 (MQ=255)
                                        cATTGCCCCGGCGCCCGGCCTGTCCGCCGCCACCGGCAGGCAACTGCGCGGCGGCGCAGACTTCGCTCCCCACCCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCt                                                                                 <  3:63346/135‑1 (MQ=255)
                                              ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGATGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGCACAGCAGGCGCTTGGCCTGCTCGCt                                                                          >  8:355212/1‑136 (MQ=255)
                                                aGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGc                                                                        <  8:441855/136‑1 (MQ=255)
                                                  gcgcCCGGCCTGTTCGGCGGCAACGGCAGGCAACGGCTCGTCGGCGCAGACTGCGCTGACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                       <  2:287281/135‑1 (MQ=255)
                                                     cccGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACGCGTCGAGGCACCGCCGGCGCGTGGCCTGCTCGCTGCCgct                                                                    >  1:44454/1‑133 (MQ=255)
                                                     cccGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGTCGAGGGACAGCCGGCGCGTGGCCGGCTCGCGGCCgcg                                                                    >  4:337536/1‑133 (MQ=255)
                                                          ccAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCACCTCGTCGAGGGACCCCCGGCGCGTGGCCGGCTCGCTGGCGCTCAGcc                                                               >  3:359615/1‑135 (MQ=255)
                                                                      cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACTCGTCGGGGCACACCAGGCGCGTGGCCGGCTCGCTGCCGCTCCGCCGTGGCCGGGGcc                                                   >  1:288342/1‑135 (MQ=255)
                                                                      cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACTCGTCGAGGAACAGCAGGCGCTTGGCCAGCTCGCGGCCGCTCAGCCCGGGCCGGTggc                                                   >  4:75444/1‑133 (MQ=255)
                                                                      cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGCGGCACACCAGGCGCGTTGCCAGCACGCTGCCCCTCCGCCCGGGCAGGTGCCa                                                  >  7:411131/1‑136 (MQ=255)
                                                                           gcaggcaACTGCTCGTCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGCTCGGGCCCCCCCCACGCGTCGACGGCCACCCCGCGCCGTGCCGGCGCCCTGCCGCTCCGCCGTGGGCGGTGCCAGcgg                                              >  4:267148/1‑133 (MQ=255)
                                                                            caggcaACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGgc                                            <  8:238191/136‑1 (MQ=255)
                                                                                  acggcgcgtcgGCGCAGACTTCGCTCCCCACCCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                       <  1:36711/129‑1 (MQ=255)
                                                                                  aCTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCCATAGGCGCCCCACAACTCGTCGAGGCACAGCACGCGCTTGGCCTTCGCGCTGCCGCTCAGCCGTGGCAGGTGCCAGCTGGCGCCgg                                       >  6:310344/1‑135 (MQ=255)
                                                                                  aCTGCTCGTCGGAGCAGACTTCGCTGACCAACCCGGAGGCCAGGGCGCGAACGGCGCCCCACAAGTCGTCGAGGAACAGCAGGCCCTTGCCCTGCTCGCTGCCGATCAGCCCTGGCACGTGCCCGCAGGCGCCgg                                       >  6:230934/1‑135 (MQ=255)
                                                                                   cTGCGCGGCGGCGCCGACTTCGCTGACCACCCCGCCGGCCAGGGCTCGAGCGGCGCCCCACAACTCGTCGCGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                     <  7:423818/136‑1 (MQ=255)
                                                                                    tGCTCTTCGGCGCAGACTTCGCTGCCCAACCCGGGGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                     <  6:92831/135‑1 (MQ=255)
                                                                                    tGCGCGGCGGCGCAGACTTCGCCGACCACCCCGGCGGCCAGGGCGCGTTCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                     <  2:326536/135‑1 (MQ=255)
                                                                                    tGCGCGGCGGCGCAGACGGCGCTCACCACCCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                     <  2:139063/135‑1 (MQ=255)
                                                                                      cTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCTCGGCGCCCCCCAACTCGTCGAGGCACAGCAGGCGCGTGGCCGGCTCGCGGCCGCTCAGCCGTGGCAGGTGGCCGCTGGCGGCGGCaac                                   >  2:138092/1‑133 (MQ=255)
                                                                                         gtcgGCGCAGACTGCGCGGACCAACCCGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                               <  8:48556/136‑1 (MQ=255)
                                                                                         ggcgGCCCATACTTCGCTGACCACCCCGGCGGCCAGGGCGCGATCGGGGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                               <  7:51772/134‑1 (MQ=255)
                                                                                         gccgGCGCAGACTGCGCGGCCCACCCCGGCGGCCAGGGCGCTATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                               <  8:55675/134‑1 (MQ=255)
                                                                                          gcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                               <  6:354424/134‑1 (MQ=255)
                                                                                          gcgGCGCAGACTTCGCGCCCCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                               <  3:381057/134‑1 (MQ=255)
                                                                                           cgGCGCAGAGTTCGCGTACCAACCGGGCGGCCAGGGGGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCt                             <  8:143212/136‑1 (MQ=255)
                                                                                            ggCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGGGCCCCCCAACGCGGGGCGGCCCAGCAGGCGCGTGGCCTGCTCGCGGGCGAGCAGCCGGGGGAGGTGCCAGCGGGCGCCGGCCTCgggggg                             >  4:344798/1‑133 (MQ=255)
                                                                                            ggCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAGGAACAGCAGGCGCTTGGGCTGCTCGCGGGCGATCAGCCGGGGCAGGTGCCCGCGGGCGCCGGCATCgggggg                             >  5:357287/1‑133 (MQ=255)
                                                                                                 aGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGCACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGcc                       >  7:254186/1‑136 (MQ=255)
                                                                                                   aCTTCGCTGACCACCACGGCGCCAAGTGGGCTAGGGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCTATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                      <  6:110080/135‑1 (MQ=255)
                                                                                                   aCTGCGCTGCCCAACCCGGCGGCCGGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCCGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGccc                      <  1:87088/135‑1 (MQ=255)
                                                                                                           gACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCAACTCGTCGACGCACAGCACGCGCGTTGCCTGCTCGCTGCCGATAAGCCCTGGCACGTGCCAGCAGGCCCCCGCATAGGGGCGGTAGCCCAAGCGGgt              >  4:57940/1‑135 (MQ=255)
                                                                                                            cccACCCGGGCGGCCGGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt            <  7:248269/135‑1 (MQ=255)
                                                                                                             ccAACCCGGGGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt            <  2:155353/135‑1 (MQ=255)
                                                                                                               aaCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGc         >  8:51568/1‑136 (MQ=255)
                                                                                                                         ccAGGGCGCGATAGGCGCCCCACAACTCGGCGAGGAACAGCAGGCGCTTTGCCTGCTCGCTTCCGATACGCCGTGGCAGGTGCCAGCGTGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCtt  >  6:12269/1‑133 (MQ=255)
                                                                                                                           |                                                                                                                                  
CGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGCCTT  >  KT2440/4245022‑4245275

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