Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,607 A→C 49.9% Q32P (CAG→CCG)  tgnC → (Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,135,6070AC49.9% 0.5 / 15.5 30Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base C (11/4);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 6.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTC  >  KT2440/4135473‑4135738
                                                                                                                                      |                                                                                                                                   
aaGCCCGCACCCTCCCCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                     <  3:192504/135‑1 (MQ=255)
aaGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                     <  3:200244/135‑1 (MQ=255)
   cccGCACCCTCACCGAATCCGCCAAGGGGATTCTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAAGGCCAATGGCGCGATGCCCAGAGAGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCAgc                                                                                                                                   >  4:263760/1‑133 (MQ=255)
     cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGCCCt                                                                                                                                >  1:25553/1‑135 (MQ=255)
                     ccGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGCCAGCCTCACCCCGGCCACCGCCCCGGCCAGGGCGCCCCAGCCCGAc                                                                                                               >  7:127867/1‑136 (MQ=255)
                     ccGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGCCAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGAc                                                                                                               >  8:74471/1‑136 (MQ=255)
                       gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGCCCGGGGCGCACCAGCCcgccgc                                                                                                             >  7:30317/1‑136 (MQ=255)
                         cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGCCGGCCAATGGCGCGATGCCCAGAGCGGCGAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCgccgcc                                                                                                            >  3:179160/1‑135 (MQ=255)
                           cggggATTGTTGCATCACCCTCGGTCGTTGCCCGATGTGCATCGACGGCCAATGGCGCGATGCCCAGCGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGa                                                                                                          <  6:182437/134‑1 (MQ=255)
                              ggtttgttgCATGACCCTGGTTCGTTTGCAGATGTGCTTCGACGGCCATTGGCGCTAGGCCCAGCGCGGCAAGACCGGCGACTGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  4:54031/132‑1 (MQ=255)
                              ggATTGTTCCATGACCATCGTGCGTTTCCATTTGGCAATCGACGGCCATTGGCGCGATGCCAAGCGCGCAAAGACCTTCGACCGCCTCAACCCGCCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  6:3582/135‑1 (MQ=255)
                               gATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCCCTGCCCGCCGCCGAGGAAg                                                                                                     >  7:195442/1‑136 (MQ=255)
                                    ttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCCCCTGCCCGCCGCCGCTGCAGCCgc                                                                                                 >  3:216176/1‑134 (MQ=255)
                                            cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCCCGGCGGCGCTGCCCGCCCCCGCGGAAGCCGAAGTggt                                                                                           >  1:272136/1‑132 (MQ=255)
                                            cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCCCCCGGGCCCAACTGCCCGACGCCGCGGAAGCGGCTGTGGAGc                                                                                         >  5:223644/1‑135 (MQ=255)
                                                          cAGATGTCCATCCACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGgcc                                                                           <  5:77921/135‑1 (MQ=255)
                                                                  cagcgccGGCCAAGGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGt                                                                  <  8:195442/130‑1 (MQ=255)
                                                                      gACGGCCAAGGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGcc                                                               <  6:46632/135‑1 (MQ=255)
                                                                              aTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGAccg                                                       >  8:267634/1‑133 (MQ=255)
                                                                                 gcgcGAGGCCCAGCCCGCCAAGACCTGCCACCGCCTCAACCCGCCCCCCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                   <  7:25154/136‑1 (MQ=255)
                                                                                            agagCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGGGCAACTGCCCGCCGCCGAGGAAGCCGATGTGGCGCTGGCCGTTCCGGCCGCCCCGCCGGCCATCGCCCGCCAGGCCTGgcgcc                                         >  6:258185/1‑134 (MQ=255)
                                                                                             gagCGCCAAGACCTGCGACAGCCGCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                        <  6:65010/135‑1 (MQ=255)
                                                                                                     aGACCTTCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCACCTGCCCGACGCCGCTGGAGCCGATGTGGCGCTGGCCGTTCAGGCCGCCCAGCGTGCCGTCGACAGCAAGGCGTGGCGCAGCATCACCg                                >  4:210053/1‑135 (MQ=255)
                                                                                                      gACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCACAGCCCGCCCCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGAGCCTTCGACAGCAAGTCCTGGCGCAGCATCAccc                                >  3:28689/1‑133 (MQ=255)
                                                                                                                 gCCTCACCCCGGCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGGGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                    <  4:235291/135‑1 (MQ=255)
                                                                                                                 cccTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                    <  3:38952/134‑1 (MQ=255)
                                                                                                                  ccTCACCCCGGCCCCCGCCCCGGCGGGGGCGCAACTGCCCGACGCCGAGGCAGCCGATGTGGAGCTGGCCGTTCCGGCCGGCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                    <  7:262016/134‑1 (MQ=255)
                                                                                                                     cAACCCGGCCACCGCCCACGCCGGGGCGCCACTGCCCGACGCCGATGAAGCCGCTGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTGGACACCCAGGCGGGGCGCCGCATCACCGCGACCGCCCGCGGcca                >  2:57145/1‑133 (MQ=255)
                                                                                                                          cGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  8:127867/136‑1 (MQ=255)
                                                                                                                           ggCCACCGCCCGGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  4:236288/135‑1 (MQ=255)
                                                                                                                              ccccGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgt      <  7:236422/134‑1 (MQ=255)
                                                                                                                               aCCGCCCCGGCCGGGGCGGCCCTGCCCGCCCCCGCGGGAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCCCCGTGGCGTCGACAGCAACGCCTGGCGCCGCCTCCCCCCGACCGCCCCCGGCCAGCTGCTGcgtc     >  8:25279/1‑136 (MQ=255)
                                                                                                                                   ccccGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCgtcgtc  <  5:251980/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTC  >  KT2440/4135473‑4135738

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: