Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 759,433 A→C 48.4% intergenic (+42/‑126) mmuM_1 → / → CJ019_00729 Homocysteine S‑methyltransferase/L‑methionine gamma‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440759,4330AC48.4% 4.0 / 18.0 26intergenic (+42/‑126)mmuM_1/CJ019_00729Homocysteine S‑methyltransferase/L‑methionine gamma‑lyase
Reads supporting (aligned to +/- strand):  ref base A (3/9);  new base C (10/3);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ACTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTGCAGCGCCCGGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGA  >  KT2440/759303‑759551
                                                                                                                                        |                                                                                                                      
aCTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGGGTTTGCAGCGCCCCGCCCGCCGCCa                                                                                                                    >  7:26901/1‑135 (MQ=255)
 cTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCCGGGCGGGTTTTGCGCGCCCCGCCCgccgccc                                                                                                                    >  3:127549/1‑133 (MQ=255)
     cGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCGGTTGGAGCGCCCGCCCCCCCGCCATTat                                                                                                                >  7:43969/1‑133 (MQ=255)
      gTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGGGGCTTTTTCGGCGCCCCGGCCGCAGCCGATaaa                                                                                                               >  7:84374/1‑134 (MQ=255)
      gTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTTGAGCGCCCGGCCCGCCGCCAATaaa                                                                                                               >  2:22065/1‑134 (MQ=255)
      gTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTGCCGCGCCCCGCCCGCCGCCATTaaa                                                                                                               >  4:51800/1‑134 (MQ=255)
      gTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCGTTTGCAGCGCCCGGCCCGCAGCCATTaaa                                                                                                               >  4:260545/1‑134 (MQ=255)
                 gCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCACCCTGGCCTAAGCGCCCACGTGGCCGGGCGGCGTTTGCCGCGCCCGGGCCGGAGCCATTAAAAAAAACGcccgg                                                                                                   >  2:47447/1‑132 (MQ=255)
                             ggcggTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGACACGCACCCTGGCCTAAGCGCCCACGTGGCCGGGCGGGTGTTGCCGCGCCCCGCCCGCCGCCAGTTTAAAAAACGCCAGGCCGCCCCggcgt                                                                                       >  7:387135/1‑132 (MQ=255)
                             ggcggTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTGCAGCGCCCGGCCCGCCGCCATTAAAAAAAACGCCAGGCAGCCCCGGCtt                                                                                       >  7:95164/1‑135 (MQ=255)
                                ggTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCCGGGGGCGTTTTCAGCGCCCCGCCCGCCGCCCGTAAAAAAAAACCCCGTCGGACCCggccttcc                                                                                    >  4:371781/1‑129 (MQ=255)
                                                  ccGGTGCAAACCGCCGTGCTGAAACGCAACCTGGCCAACGCGCCCACGTGGCCGGGCGGTTTTTGCAGCGCCGGCCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCg                                                                  <  3:387312/135‑1 (MQ=255)
                                                   cGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCGTTTTGAGCCCCCCGCCCGCCGGCCTGAAAAAAAAACGCACGCAGACCCGGGTGTCCTCGTGTTTGAGGccggcct                                                                 >  8:163137/1‑130 (MQ=255)
                                                                  aGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTTGAGCGCCCCGCCCGCCGCCATTAAAAAAAACGCCAGTCAGCCCCGGATTACCTCGTTTTTGAGGACGCCCTCATCTCACTATAGta                                                  >  6:361249/1‑134 (MQ=255)
                                                                    tgtaaaCGCACCTTGCCATAGGCGCCCAGGGGCTGTGGGGGCTTTTGGCCCGCCGGGCCCGCAGCCATAAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCg                                                <  4:104919/132‑1 (MQ=255)
                                                                    cTGAAACGCAACCTGGCCTAAGGGCCCACTTGGCGGGGGTGCTTTGGCCGCCCCGGGCCCGCCGCCTTTAAAAACAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCg                                                <  6:171031/135‑1 (MQ=255)
                                                                             caCCCGGCCTAACGCCCCCAGGGGCGGGGGGGGTTTTGCCCCCGCCGGCCCCGAACCAATTAAAAAACACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGc                                        <  7:179103/134‑1 (MQ=255)
                                                                                   gCCTAAGCGCCCACGGGGCCGGGCGGCTTTTGCCGCGCCCGGCCCGCAGCCATAAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGAtgct                                 <  8:95164/135‑1 (MQ=255)
                                                                                   cccTACGCGCCCCGGGGGCGGGGCGGCTTTTGCCGCGCCCGGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGAtgct                                 <  4:100192/134‑1 (MQ=255)
                                                                                   ccatacgcgcCCACGGGGCGGGGCGGTTTTTCCCGCGCCGGGCCCGCCGCCATAAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGAtgct                                 <  8:84374/129‑1 (MQ=255)
                                                                                    ccaatgcgcCCCCGGGCCGGGGGGCTTTTACCAGCCCCCGGCCCCCACCATTTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGAtgctg                                <  2:320268/130‑1 (MQ=255)
                                                                                      tacgcgcCACCGGGGCCGGGCGGCTTTTCCAGCGCCCGGCCCGCAGCCATTAAAAAACACGCCAGTCAGACCCTGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCa                              <  4:41646/132‑1 (MQ=255)
                                                                                        aGCCCCCACGGGGCCGGGGGGATTTGGCCGCGCCCGGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATc                            <  8:326553/135‑1 (MQ=255)
                                                                                                                      gCAGCGCCCCGCCCGCCGCCATTTAAAAAAAAGCCAGTCAGACCCGGATTCCATCGTTGTTAAGGCCTCCGTGAACTCCCTATAGTACCCTGAAGGCCGATGCTGCATCGGTGTACCCAACTCACCTGGGCgggc  >  4:267312/1‑134 (MQ=255)
                                                                                                                       gatgggccgggcccgCTGCCATAAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCggg   <  4:381056/129‑1 (MQ=255)
                                                                                                                       cccccgcCGGCCCCGCAGCCAATAAAAAACACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCggg   <  3:51800/130‑1 (MQ=255)
                                                                                                                        gccgccccggccccGCCCCCATAAAAAAACACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCggg   <  5:361249/129‑1 (MQ=255)
                                                                                                                                        |                                                                                                                      
ACTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAAGCGCCCACGTGGCCGGGCGGCTTTTGCAGCGCCCGGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGA  >  KT2440/759303‑759551

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: