Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 382,484 T→G 39.9% S14A (TCC→GCC)  soxD → Sarcosine oxidase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440382,4840TG39.9% 18.8 / 5.0 40S14A (TCC→GCC) soxDSarcosine oxidase subunit delta
Reads supporting (aligned to +/- strand):  ref base T (8/14);  new base G (8/8);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGAA  >  KT2440/382354‑382612
                                                                                                                                   |                                                                                                                                
ccGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGCTCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGCCACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCGcccg                                                                                                                              >  2:102797/1‑133 (MQ=255)
                     gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCCCTGCGGCGCGCGGCGCGCCGGAGAAGAGTTCCACGCCGCTGGc                                                                                                        >  3:320235/1‑136 (MQ=255)
                     gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGc                                                                                                        >  4:258097/1‑136 (MQ=255)
                         caaCAACGGCGCCCTGATCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGGCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCCTCTGGCCa                                                                                                      <  8:269290/132‑1 (MQ=255)
                          tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGccccgc                                                                                                   >  4:228597/1‑132 (MQ=255)
                                     aCTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCCCTGCGCGCCGAAGAAGACTTACACGCCACAGGCCAGGCGCACATTgcc                                                                                         >  6:181519/1‑135 (MQ=255)
                                     aCTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTgcc                                                                                         >  8:272799/1‑135 (MQ=255)
                                           cgcccaaCCCGGCGCCGCCGCCGGCGCCCAATAACCG‑GAGACACCGTCTTGTTGCATATTTTCTGCCCCCCCTCCGGCGAGCGGCGCGCCGAAGACGAGGTCCACGCCTCTGGCCAGGCGCACATTGcccgcccg                                                                                   <  7:17811/130‑1 (MQ=255)
                                           cgacgaACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGcccgcccg                                                                                   <  2:244083/135‑1 (MQ=255)
                                            gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCccgcccgc                                                                                  >  8:223682/1‑135 (MQ=255)
                                            gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCCCAGAAGAGATCCACCCCCAATGCCACGCGCACAATgcccgccccc                                                                                  >  1:284597/1‑133 (MQ=255)
                                            gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCCCAGAAGAGTTACACGCCTCTGGCCCCGCGCACCATGCCCGCCCCCt                                                                                 >  4:23839/1‑136 (MQ=255)
                                                          cgccgccgTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGGGCGGGGCGGCGCAGAAGAGTTGCCCCCCGCTGGCCCCGCGCCCAGTGCCCCCCCCCGGGGCCCGAACGcc                                                                    >  1:314997/1‑135 (MQ=255)
                                                           gccgccgGCGCCAAC‑AACCGCGCCACCCCGTCTTTTTTCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCGAAGAACAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  1:349153/135‑1 (MQ=255)
                                                                cgTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAACAAGAGTTACACGCCTCTGGCCAGGCGCACATTTCCCGCCCGCTGGCCCCTAAAGCCTGCTCcg                                                             >  4:158805/1‑136 (MQ=255)
                                                                  gcGCCCACTCCCCG‑GCCACACCGTCTTGTAGCATATTTTCTGCCCCCACGGCGGCGCGCTGCGCGCCGACGAAGCGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcgac                                                           <  3:310920/135‑1 (MQ=255)
                                                                       cACTAACCG‑CAGACACCGTCATGTTGCATATTTTCTCCCCCCACGGCGGCGAGCTGCGCTCCTAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgagga                                                      <  3:23839/136‑1 (MQ=255)
                                                                         caaaCCG‑GACACACCGTCATGTTGCATTTTTTCTGACCCCACTGCGGCGAGCGGCGCGCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGt                                                    <  3:258097/134‑1 (MQ=255)
                                                                          gaaCCG‑GAGACACCGTCATGTTGCATATTTTCTGGCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGt                                                    <  8:22944/134‑1 (MQ=255)
                                                                               cgttcccacCGTCATGTTGCATATTTTCTGCCCCCCCTGCGGCGCGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTgggg                                                <  5:181519/129‑1 (MQ=255)
                                                                                 gagaCACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCGCCGAAGAAGAGTTCCACGCCGATGGCCAGGCGCACATTTCCCCCCCCCTGGCACCTCACGCCTGCTGCGACGAGGAGGGGGGTaa                                             >  4:130445/1‑135 (MQ=255)
                                                                                 ccgaCACCGTCAGGTTGCATATTTTCTGTCCCCACTGCGGCGAGGTGCGCTCCTAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTa                                              <  8:279017/133‑1 (MQ=255)
                                                                                   gaCCCCGTCAGGTGGCATATTTGCTGGCCCCACTGCGGCGAGCTGCGCCCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  6:299670/135‑1 (MQ=255)
                                                                                   gaCACCTTCATTTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGACGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  7:3065/135‑1 (MQ=255)
                                                                                   gaCACCGTCATGTTGAATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  5:152695/135‑1 (MQ=255)
                                                                                   gaCACCGTCATGTTCCATATTTTCTGTCCCCCCTCCGGCGCGCTGCGCTCCGAAGAAGTGGCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  5:152701/135‑1 (MQ=255)
                                                                                           cATGTTGAATATTTTCTGCCCCCCCGGCGGCGCGCGGCGCTCCGAAGAAGAGGCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGtt                                    <  2:314997/135‑1 (MQ=255)
                                                                                            aTGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCCGGCGCACATTGCCCGCCCGCGGGACCCTAACGCCTGCTCCGCCGAGGAGTGGGGTACCTACATGTTc                                   >  1:318047/1‑135 (MQ=255)
                                                                                            aTGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTc                                   >  4:35213/1‑135 (MQ=255)
                                                                                                tGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTAc                                <  7:223682/134‑1 (MQ=255)
                                                                                                     atTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTACACGCCTCTGGCCCGGCGCACATTTCCCGCCCGCGGGACCCGAACGCCTGCGCCGCCGAGGAGTGGGGTAACTACATGTTCTACCCTGa                           >  8:9418/1‑134 (MQ=255)
                                                                                                         tCTGGCCCCACTGCGCCGAGCTGCGCGCCGAAGAAGACTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGT‑GGGTACCTACATGTTCTACCGTGACAAccc                     <  3:111130/135‑1 (MQ=255)
                                                                                                          cTGTCCCCACTGCGGCGCGCTGCGCGCCCAAGAAGAGTTCCACGCCTCTGGCCACGCGCACATTGCCCGCCCGCGGGCCCCGAACGCCTGCGCCGACGAGGAGTGGGGGACCGACAAGTTCTACCGGGACCAccc                     >  7:81813/1‑135 (MQ=255)
                                                                                                          cTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                     <  4:138016/135‑1 (MQ=255)
                                                                                                          cTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                     <  3:35213/135‑1 (MQ=255)
                                                                                                           tGTCCCCTTTCCGGCGAGCGGCGCTCTGAATAAGAGTGCCTCGACTCTGGGCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                     <  1:336629/134‑1 (MQ=255)
                                                                                                               cccACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCgg               <  4:145591/136‑1 (MQ=255)
                                                                                                                 cACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTa             <  4:127618/136‑1 (MQ=255)
                                                                                                                      cggcgAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTc          <  4:42435/134‑1 (MQ=255)
                                                                                                                             cggcgcACCGATGCTGAGGTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGaa  <  6:129876/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
CCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCG‑GAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGAA  >  KT2440/382354‑382612

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