Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,607 A→C 45.7% 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,6070AC45.7% 8.5 / 8.4 38Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (8/11);  new base C (14/4);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTC  >  KT2440/4135473‑4135738
                                                                                                                                      |                                                                                                                                   
aaGCCCGCACCCTCCCCGATTCCGCAAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                     <  5:280334/135‑1 (MQ=255)
aaGCCCGCACCCTCACCGAATCCGCCACGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                     <  2:431745/135‑1 (MQ=255)
          ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCACCCCGGGgccc                                                                                                                          >  8:307441/1‑134 (MQ=255)
           cTCGCCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGATCGTTTCCAGATGTGCATCGACGGCCAAAGGCGCGATGCCCAGCGCGTCTAGACCTTCGACCGCCGCCACCCGGCCACCCCCCCGGCCAggggcc                                                                                                                          >  2:180975/1‑132 (MQ=255)
               ccGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACt                                                                                                                      >  4:95317/1‑135 (MQ=255)
                     ccGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCCACCCGGCCACCGCCCCCCCCCGGGCGCCCCAGCCCgc                                                                                                                >  1:170712/1‑134 (MQ=255)
                           aGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCCGGGCGCACCTGCCCGCCGCCGAt                                                                                                         >  7:208184/1‑136 (MQ=255)
                              ggATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  2:292369/135‑1 (MQ=255)
                              ggATGGTTGCTTGACCCTCGTGCGTTCCCATATGTGCAACGACGGCCAATGGCGCGTGGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  6:1943/135‑1 (MQ=255)
                                   gttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCCACTGCCCGCCCCCGATGCAGcc                                                                                                   >  4:438101/1‑134 (MQ=255)
                                     tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACAGCCCGCCGCCGATGAATcccaag                                                                                               >  7:147755/1‑132 (MQ=255)
                                          gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCCCCGGGGCGCCCCTCCCCCCCGCCCACCAAGCCCAAGTGGAg                                                                                          >  7:80024/1‑136 (MQ=255)
                                                 gttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGCCAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACGGCCCGCCGCCGAGGAAGCCGCTGTGGGGCTGGcc                                                                                    >  1:91552/1‑135 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGCCGCCGCTGCAGCCGATGTGGCGCTGGCCGGTGAGgcc                                                                           >  8:446469/1‑136 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACAGCCCGCCGCCCCTGAAGCCCCTGTGGAGCTGGCCGTTCAGgc                                                                            >  2:6971/1‑135 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACAGCCCGACGCCGATGAAGCCGCTGTGGAGCTGGCCGTTCAGgc                                                                            >  6:339277/1‑135 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCCACTGCCCGCCCCCCAAGGAGCCGATGTGGAGCTGGCCGATAAGgcc                                                                           >  8:295165/1‑136 (MQ=255)
                                                         ccAGATGTGCATCGACGGCCAATGGCGCGATGCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCACCTGCCCCACGCCGATGCAGCCGCTGTGGAGCTGGCCGTTCAGgcc                                                                           >  4:425510/1‑135 (MQ=255)
                                                             aggcgCATCCCCGGCCACTGGCGCGCGGCCCAGAGCGGCAAGCCCTACGACCCCCTCACCCCGGCCCCCGCCCCGGCCGGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                                       <  7:107368/132‑1 (MQ=255)
                                                                    tCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCGCCACGGCGCCCCTGCCCCCCGCCGCTGAAGCCGATGTGGGGCGGGCGGTTCTGGCCGCCCCGCGAGc                                                                >  8:21796/1‑136 (MQ=255)
                                                                             aaTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGGGCCCCTGCCCGCCGCCGCTGAAGCCGCTGTGGAGCTGGGCGGTGAGGCCGCCCCGCGGGCCGTCGCCa                                                        >  1:247214/1‑135 (MQ=255)
                                                                                 ccgcgAGGCCCCGCGCGGCAAGACCTTCGCCAGCCTCACCCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                   <  8:233975/135‑1 (MQ=255)
                                                                                  cgcgCTGCCCACCGCGGCACGACCAGCGACCGCCTCACCCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGCCGCCGATGAAGCCGAGGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                   <  6:275076/135‑1 (MQ=255)
                                                                                     gATGCCCCGAGCGGCAAGCCCTTCGACCGCCGCCACCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCg                                                <  5:354613/135‑1 (MQ=255)
                                                                                       tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGCTGTGGCGCTGGCCGTTCAGGCCGCCCCGCGGGCCGTCGCCAGCAAGGCGTg                                              >  5:223691/1‑135 (MQ=255)
                                                                                       tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGCCGCCGATGACTCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCGTGGCCAGCAAGGCGTgg                                             >  7:322178/1‑136 (MQ=255)
                                                                                           cAGAGCGCCAAGACCTTCGACAGCCTCACCCCCGCCACCCCCCAGGCCGGCGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                          <  1:157485/135‑1 (MQ=255)
                                                                                             gcgCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                        <  2:247214/133‑1 (MQ=255)
                                                                                                    aaGCCCTTCGACAGCCTCAACCCGGCCACACCCCCGGCCGGGGCCCAACAGCCCGCCGCCGAGGACGGCGCTGTGGAGCGGGCCGTTCAGGCCCCCCAGCGGTCCTGCGCCCGCCAGGCGTGGCGCGCCGTCAcc                                 >  6:112932/1‑135 (MQ=255)
                                                                                                              aCAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgccc                       <  4:278921/135‑1 (MQ=255)
                                                                                                                cgCCTCACCCCGCCCCCCGCCCAGCCCGGGGCGCAACTGCCCGACGCCGATGCAGCCGAGGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgcccgc                     <  5:131529/134‑1 (MQ=255)
                                                                                                                 gCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                    <  2:250779/135‑1 (MQ=255)
                                                                                                                     cAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGGGCTGGCCGTTCAGGCCGCCCAGCGTGGCGTCGACAGCAAGGCGGGGCGCAGCATAACCCCTACCCCCCGCGGCCAg               >  8:230744/1‑136 (MQ=255)
                                                                                                                         gcGGCCCCCGCCCAGGCGGGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg            <  5:304998/134‑1 (MQ=255)
                                                                                                                           tgCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCGTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  2:170712/134‑1 (MQ=255)
                                                                                                                           ggCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  6:340284/135‑1 (MQ=255)
                                                                                                                              ccccGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgt      <  8:80024/134‑1 (MQ=255)
                                                                                                                               aCCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGTGGCGCCGCATCACCGCGACCGCCCGCGGCCAGCGGCTGcgt      >  2:59710/1‑135 (MQ=255)
                                                                                                                               aCCGCCCCGGCCTGGGCGCAACTGCCCCCCGCCGCGGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGGCGGGGACAGCAAGGCGTGGCGCCGCCTCACCGCTACCGCCCCCGGCCAgcggcggcgg      >  2:15706/1‑134 (MQ=255)
                                                                                                                                   cccTGGCCTGGGCGCAACTGCCCGACGCCGATGAATCCGTTGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCTCTACCGCCCGCGGCAAGCTGCTGCgtcgtc  <  4:382722/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 < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: