Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 382,476 A→C 43.9% E11A (GAG→GCG)  soxD → Sarcosine oxidase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440382,4760AC43.9% 7.9 / 2.4 26E11A (GAG→GCG) soxDSarcosine oxidase subunit delta
Reads supporting (aligned to +/- strand):  ref base A (4/9);  new base C (7/5);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACC  >  KT2440/382345‑382607
                                                                                                                                   |                                                                                                                                    
acccACGGGCCGCGCCGTTCCCCCTGGACCGTTTCAACAAGGGCGCACTGAGCCACGACCCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                   <  5:193356/133‑1 (MQ=255)
aTCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTAAAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                   <  4:18695/135‑1 (MQ=255)
                gTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCCCTAACCGGAGACCCCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTTCGCTCCGAAGAagag                                                                                                                   >  8:240054/1‑135 (MQ=255)
                  gccccATGCAGCTTTTCAACACGGGCGCACTGTGCGACAACCCGGGCGCCGCCGCCGTCGCCCACTAGCGGGACACACCGTCTTGTTGCATATTTTCTGGCCCCACGGCGCCGCGCTGCGCTCCTAAGAAGAGtc                                                                                                                 <  6:243027/132‑2 (MQ=255)
                            ccGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCAGCCGTCGCCCACTAACCGGAGACCCCGTCATGTTTCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGCATACTTCCCCCCct                                                                                                        >  6:92414/1‑133 (MQ=255)
                              gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGCAGAAGAGTTACACGCCGCTgg                                                                                                    >  8:56104/1‑135 (MQ=255)
                              gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCCCTGCGCTCCCAAGAAAAGGTCCAccccacagg                                                                                                    >  8:290393/1‑130 (MQ=255)
                              gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTgg                                                                                                    >  7:97165/1‑135 (MQ=255)
                              gTTTCTACAACGGCGAACAGATCGACGAACACGGCGCCGCCGCCGTCGCCCCCTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGTGCTGCGCTCCCAAGAACAGTTACCCCCctctct                                                                                                    >  6:203250/1‑133 (MQ=255)
                                   tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCAGCGCGCCCAAGCAGAGTTACCACCCACTGGcccgc                                                                                               >  3:9458/1‑132 (MQ=255)
                                              agtGATCGCGGACCCCGGCGCCGCCGCCGTCGCCCACTACCCGGCCACACCGTCAGGTTCAATATTTTCTGCCCCCACTGCGGCGCGCTGCGCTCCGAACAAGCGTTCCACGCCTCTGGCCAGGCGCACATTgcc                                                                                    <  2:234/133‑1 (MQ=255)
                                                   tCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCGGCGCTCCGCAGAAGAGTTCCACGCCTCTGGCCCGGCGCACATTgcccgccc                                                                               >  7:27830/1‑135 (MQ=255)
                                                          acacGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                        <  1:13743/135‑1 (MQ=255)
                                                                    gccgccgTCGCCCACTAACCGGAGACACCTTCTTGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                              <  4:183509/135‑1 (MQ=255)
                                                                         cgccGCCCACTAACCGGAGACACCGTCAGGTTGGATTTTTTCTGTCCCGACTGCGCGGAGCTGCGCTCCGAAGACGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTcc                                                         <  4:274679/132‑1 (MQ=255)
                                                                          ggcGCCAACAACCCGGACACCCCGTCTTGTTGCATTTTTGCGGCCCCCCCGGCGCCGAGCGGCGCCCCGAAGAAGGGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                        <  5:249521/133‑1 (MQ=255)
                                                                            cGCCCACTACCCGGAGACACCGGCAGGTTGCATATTTTCTGGCCCCACCGCGGCGCGCGGCGCGCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcgac                                                      <  6:86/135‑1 (MQ=255)
                                                                                    acccGGAGACACCGTCATGTTGCATATTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTg                                              <  4:9458/133‑1 (MQ=255)
                                                                                       cGGATACACCGTGATTTGGCATATTTTCTGGCCCCCCTGCGGCTCGCTGCGCGCCGAAAACGCGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTgggg                                           <  6:232152/135‑1 (MQ=255)
                                                                                             cacGGTCATGTTGCATATTTGCTGTCCCCAGTCCGTCGCGCTGCGCACCGAAGACGAGCTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                       <  6:146635/133‑1 (MQ=255)
                                                                                                    aTGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCCGGCGCACATTGCCCGCCCGCTGGCCCCTAACGCCTGCTCCGACGAGGCGTGGGGTACCTAAATGTTc                              >  6:48031/1‑135 (MQ=255)
                                                                                                             atTTTCTGTCCCCACTGCGGCGCGCTGCGCGCCGAAGAAGAGTTCCCCCCCTCTGGCCCCGCGCACATAGCCCCCCCCCGGGCCCCTAACGCCGGCTCCGACCCGGAGTGGGGTACCGACATGTTCTACCGTgc                      >  8:25439/1‑133 (MQ=255)
                                                                                                                  cTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                >  7:40463/1‑135 (MQ=255)
                                                                                                                   tGTCCCCACTGCGGCGAGCTGCGCTCCGACGACGCGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                <  5:141261/134‑1 (MQ=255)
                                                                                                                       cccACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGCCGAGGAGTGGGGTACCTACATGTTCTACCGGGCCAACCcgcgcg           >  1:160162/1‑135 (MQ=255)
                                                                                                                                gcgAGCTGCGCTCCGAAGAAGAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCAcc  <  6:98515/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
ATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACC  >  KT2440/382345‑382607

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