Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,514,652 Δ1 bp 10.8% coding (731/1257 nt) icd_1 ← Isocitrate dehydrogenase [NADP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,514,6520T.10.8% 112.2 / 9.0 37coding (731/1257 nt)icd_1Isocitrate dehydrogenase [NADP]
Reads supporting (aligned to +/- strand):  ref base T (13/20);  new base . (2/2);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGA  >  KT2440/4514518‑4514784
                                                                                                                                      |                                                                                                                                    
gcagcaTGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCt                                                                                                                                      >  6:124913/1‑135 (MQ=255)
gcagcaTGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCt                                                                                                                                      >  2:30743/1‑135 (MQ=255)
 cagcaTGGCGGCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGACATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCtt                                                                                                                                     <  4:18840/135‑1 (MQ=255)
      gggcggcggcGATGGCGTCCTGGACGATGACCGCGCGACCGGTCTTGGGGGTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAgg                                                                                                                                <  4:19746/134‑1 (MQ=255)
             ggCGATGGCTTCCTGCACGATGACCACGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGACGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCg                                                                                                                        <  7:277495/136‑1 (MQ=255)
                 aTGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTGAAGGCACCCTCGGTGa                                                                                                                    >  3:215114/1‑135 (MQ=255)
                 aTGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTGAAGGCACCCTCGGTGa                                                                                                                    >  3:218888/1‑135 (MQ=255)
                               gATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTg                                                                                                      <  8:137592/136‑1 (MQ=255)
                                  gACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTAGcc                                                                                                    <  6:109754/135‑1 (MQ=255)
                                        gcgACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc                                                                                             <  7:236318/136‑1 (MQ=255)
                                        gcgACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc                                                                                             <  4:215114/135‑1 (MQ=255)
                                         cgACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc                                                                                             <  4:218888/134‑1 (MQ=255)
                                                   tGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCACGATGTTGCCCTTGTGCACCAGGGACa                                                                                   >  6:57565/1‑135 (MQ=255)
                                                             gAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTcgcg                                                                         >  4:273805/1‑135 (MQ=255)
                                                             gAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCgg                                                                        >  7:237410/1‑136 (MQ=255)
                                                             gAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGATGCCCTGGTGCACCAGGGTCAGTGATTcgag                                                                         >  2:116264/1‑133 (MQ=255)
                                                                      ccAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATcc                                                               <  7:207801/136‑1 (MQ=255)
                                                                          gggCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATccacc                                                            <  2:19859/135‑1 (MQ=255)
                                                                          gggCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATccacc                                                            <  1:155382/135‑1 (MQ=255)
                                                                          gggCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATccacc                                                            <  5:247967/135‑1 (MQ=255)
                                                                                ccATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTAt                                                      >  4:125300/1‑135 (MQ=255)
                                                                                 cATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTAtt                                                     <  5:124913/135‑1 (MQ=255)
                                                                                 cATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTAtt                                                     <  1:30743/135‑1 (MQ=255)
                                                                                      aGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAgg                                                <  3:116222/135‑1 (MQ=255)
                                                                                                cgcCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCTCCCTCGGTGAACTTCACGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCa                                      <  5:57565/135‑1 (MQ=255)
                                                                                                cgcCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCAc                                     <  8:61376/136‑1 (MQ=255)
                                                                                                   cGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCa                                   <  1:98124/135‑1 (MQ=255)
                                                                                                   cGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAg                                  <  8:237410/136‑1 (MQ=255)
                                                                                                        tCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGc                                                                                                                               <  5:247443/38‑1 (MQ=255)
                                                                                                        tCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGc                                                                                                                               >  6:247443/1‑38 (MQ=255)
                                                                                                        tCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGc                              >  4:29537/1‑135 (MQ=255)
                                                                                                        tCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCt                             >  8:214383/1‑136 (MQ=255)
                                                                                                           gcACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCTTgg                          <  8:102788/135‑1 (MQ=255)
                                                                                                                  ccACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTGATCCACCACGGATTGCAGGGCCTTGCGCACCACCCGCTTGGTACcct                    >  6:34814/1‑135 (MQ=255)
                                                                                                                            cccccAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACAcc         <  7:214383/136‑1 (MQ=255)
                                                                                                                                cAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTCCCCTCCTTCGACACCGGc      >  2:77841/1‑135 (MQ=255)
                                                                                                                                  gTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc                                                                                             >  2:241408/1‑46 (MQ=255)
                                                                                                                                  gTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc                                                                                             <  1:241408/46‑1 (MQ=255)
                                                                                                                                   tCCTTGAAGGCACCCTCGGTGAACTTCATGATCTTGCCCAAGTGCACAAGGCTCAGTGATTCGCGGTCGTTAAACACCACGTCTTGCAGGGCCATGCGCACCAGACGCTTCCTACCCTCCATCCACACCGGCTTg   >  6:70219/1‑135 (MQ=255)
                                                                                                                                    ccTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATACACCACGTATTGCAGGGCCTTGCGCACCAGCCGCTTGGTACCCTCCTTCGACACCGGCTTGa  >  6:107429/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGCACCAGGGTCAGTGATTCGCGGTCGTTATCCACCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGA  >  KT2440/4514518‑4514784

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