Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,607 A→C 47.4% 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,6070AC47.4% 3.8 / 15.7 53Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base C (18/7);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 5.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135475‑4135741
                                                                                                                                    |                                                                                                                                      
gCCCGCACCCTCACCGAATCCGCCAAGGGGATGGTGGCATGCCCCTCGTGCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAgg                                                                                                                                      <  1:75863/135‑1 (MQ=255)
   cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGccg                                                                                                                                   >  1:4614/1‑134 (MQ=255)
   agCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATTGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCCACCCGGCCACCCCCCCGCCCt                                                                                                                                   >  4:208855/2‑135 (MQ=255)
         cTCCCCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGgcgc                                                                                                                             >  1:113673/1‑135 (MQ=255)
               gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCTGGGCGCCACGGc                                                                                                                       >  3:202923/1‑135 (MQ=255)
               gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGACGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCTGGGCGCAACTGc                                                                                                                       >  8:118891/1‑135 (MQ=255)
                     gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCCCCCAGGCCTGGGCGCACCTGCCCGAcc                                                                                                                 >  2:16626/1‑134 (MQ=255)
                     gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGCCCTGGGCGCAACTGCCCCACg                                                                                                                 >  1:56871/1‑135 (MQ=255)
                           gggATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCGTGGGCGCAACTGCCCGACGCCGATg                                                                                                           <  1:54872/135‑1 (MQ=255)
                            gttttgttgCATGACCCTCGTTCGTTTCCATATGGGCATCGCCGGCCACGGGCGCGATGCCCAGAGCGGCAACACCTGCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                          <  4:132115/132‑1 (MQ=255)
                            ggtttgttgCTTGACCCTCGTTGGTTGCCATATGTGCATCGACGGCCAATGGCGCGAGGCCCAGAGCGCCAAGACCTACGACCGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                          <  2:134656/132‑1 (MQ=255)
                            ggATTGTTGCATGACCCTCGTTCTTTTCCTGATGGGCATCGCCGGCCAATGGCGCGATGCCCAGAGCGCCAATACCTTCGACAGCCTCAACCCGGCCACCGCCCTGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                          <  4:94575/135‑1 (MQ=255)
                            ggATTGTTGCATGACCCGGGTTCGTTTCCAGTTGTGCTTCGACGGCCAAGGGCGCGAGGCCCAGAGCGGCAAGACGTTCGACAGCCTCACCCGGGCCACGGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                          <  4:84917/135‑1 (MQ=255)
                                  ttgCATGACCCTCGTGCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCg                                                                                                     <  8:250846/134‑1 (MQ=255)
                                   tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGAt                                                                                                   >  6:240007/1‑135 (MQ=255)
                                   tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCACGGCGCCACAGCCCCACGCCGCAGACGCCGaa                                                                                                   >  1:194854/1‑134 (MQ=255)
                                                 tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACAGCCCGCCGCCCAAGAAGCCGATGTGGGGCTGGCCGt                                                                                     >  6:104943/1‑135 (MQ=255)
                                                 tcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCAGGGCGCAACTGCCCGACGCCGCAGAAGCCGATGTGGAGCTGGCCg                                                                                      >  8:132314/1‑134 (MQ=255)
                                                  cgttCCCAGCTGCGCCGCGCCGCCCAAGGGCGCGAGGCCCCGCGCGCCAACCCCTGCGACCCCCTCACCCCGGCCACCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGtt                                                                                    <  7:165472/135‑1 (MQ=255)
                                                        cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGGGCCACTGCCCGCCGCCGCAGACGCCGCGGGTGAGCTGGCCGTGCCGgcc                                                                              >  5:35874/1‑135 (MQ=255)
                                                        cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCCGGGCGCCACAGCCCCACGCCCCGGCACCCGATGTGGCGCGTGCCGGTCACgcc                                                                              >  2:16213/1‑135 (MQ=255)
                                                         aGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCCCCCCGGCGCCACAGCCCGCCCCCCATGACGCCCAGGTGGAGCGGGCCGATCAGtccg                                                                             >  7:180290/1‑135 (MQ=255)
                                                                 aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACTGCCCGCCGCCGCTGAAGCCGCTGTGGCGCGGGCCGGTCAGGCCGccccccgg                                                                     >  6:287829/1‑130 (MQ=255)
                                                                      cGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGACTCAACCCGGCCACCGCCCAGGCCTGGGCGCGACTTCCCGCCGCAGATGAAGCCGATGTGGCGCTGGCCGTTCAGGCCGCCCCGCGTGccgt                                                                >  7:297011/1‑133 (MQ=255)
                                                                      cGGCCAAGGGCGCGAGGCCCAGCGCGGCTAGCCCTGCGCCAGCCGCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCtt                                                                <  7:32580/135‑1 (MQ=255)
                                                                           aagggCGCGATGCCCAGCGCGCCAAGACCTTCGACAGCCTCACCCCGGCCCCCCCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  8:162995/132‑1 (MQ=255)
                                                                           aaggccgcgAGGCCCAGCCCGCCACCACCTTCGACCCCCGCACCCCGCCCCCCGCCCCGGCCGGGGCGGAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  2:60265/130‑1 (MQ=255)
                                                                           aaTGGCGCGATGCCCAGAGCGGCAACACCTGCGACCCCCCCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  1:178985/135‑1 (MQ=255)
                                                                           aaTGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  8:292662/135‑1 (MQ=255)
                                                                                   gATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGCTGTGGCGCTGGCCGTTCAGGCCGCCCCCCGTGCCGTCGCCCCCAAGGCg                                                   >  2:281720/1‑135 (MQ=255)
                                                                                     tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGCCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  8:3838/135‑1 (MQ=255)
                                                                                     ggCCCAGAGCGGCAAGCCCTGCGACAGCCTCACCCCGGCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  2:169040/134‑1 (MQ=255)
                                                                                        ccAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCCCAGCCCGCCGCCGCGGAAGCCGCTGTGGGGCTGGCCGGTCAGGCCGCCCCGCGGGCCGTCGACAGCCACGCGTGgcg                                              >  1:264018/1‑135 (MQ=255)
                                                                                        ccAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGCAGCCGATGTGGCGCTGGCCGTTCAGGCCGCCCCCCGTGCCTTCGCCCGCAAGGCGGGgcg                                              >  6:235756/1‑135 (MQ=255)
                                                                                           gcgcgcCAAGCCCTTCGACCCCCCCACCCCGCCCCCCCCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                           <  8:44759/133‑1 (MQ=255)
                                                                                           gagCGGCACGCCCTGCCACCGCCTCAACCCCGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                           <  4:128892/135‑1 (MQ=255)
                                                                                           gagCGGCAAGACCTTCGACAGCCGCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                           <  2:42986/135‑1 (MQ=255)
                                                                                              cGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCAGGGCGCAACAGCCCGCCGCCGCAGAAGCCGCTGTGGAGCTGGCCGTTCAGGCCGCCCCGCGGGGCTTGGCCCGCAAGGCGTGGCGCAGCAt                                        >  2:217936/1‑135 (MQ=255)
                                                                                               ggCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCCACAGCCCGCCGCCGCTGAAGCCGATGTGGAGCT‑GCCGTTCAGGCCGCCCCGCGTGCCGTCGACAGCAAGGCGTGGCGCAGCATCa                                      >  1:49947/1‑135 (MQ=255)
                                                                                                    gACCTTCGACAGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGGTGCGCCGCCTCCCCGc                                  >  6:3224/1‑135 (MQ=255)
                                                                                                    gACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCACTGCCCCCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCCCCCACCCTGCCGTAGCCCGCAAGGCGTGGCGCAGCATCACCGc                                  >  2:293950/1‑135 (MQ=255)
                                                                                                    gACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCAGGGCGCCACAGCCCGCCGCCGCTGAAGCCGCTGTGGAGCTGGCCGTTCAGGCCGCCCACCGGGCCGTCGCCAGCCAGGCGTGGGGCAGCATCAccc                                   >  3:103043/1‑133 (MQ=255)
                                                                                                    gACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCTGGGCGCAACTGCCCGCCGCCGCAGGAGCCGATGTGGAGCTGGCCGTTCAGGCCCCCCAGCGGGCCGTAGTCAGCAAGGCGGGGCGCCGCAACAACGc                                  >  8:190187/1‑135 (MQ=255)
                                                                                                       cttccccAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAc                               <  4:202923/129‑1 (MQ=255)
                                                                                                        ttCGACAGCCTCAACCCGCCCACCGCCCAGCCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  4:152933/135‑1 (MQ=255)
                                                                                                               gCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  7:300433/135‑1 (MQ=255)
                                                                                                               gCCGCAACCCGGCCACCGCCCAGCCCTGGGCGCACCTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  6:249964/135‑1 (MQ=255)
                                                                                                                   cAACCCGGCCACCGCCCCCGCCGGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGGTGAGGCCGCCCCGCGTGGCGTCGCCAGCCAGGCGTGGCGCCGCCTCACCCCTACCCCCCGCGGCaa                   >  2:282135/1‑135 (MQ=255)
                                                                                                                    accccGCCCACCGCCCGGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAg                  <  3:36811/133‑1 (MQ=255)
                                                                                                                    accccGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAg                  <  8:180290/133‑1 (MQ=255)
                                                                                                                       ccGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg               <  4:30123/135‑1 (MQ=255)
                                                                                                                         ggCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct             <  6:35874/135‑1 (MQ=255)
                                                                                                                              ccGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgtc        <  5:235756/135‑1 (MQ=255)
                                                                                                                                    aGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGGGCTGGCCGTTAAGGCCGCCCCGCGTGCCTTCGACCGCCAGGCGTGGCGCCGCCTCACACCTACCGCCCGCGGCAAGCGGCTGcgtcgggagg  >  6:197645/1‑130 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  KT2440/4135475‑4135741

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