Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,689,816 A→C 51.5% S236R (AGC→CGC)  pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,689,8160AC51.5% ‑4.2 / 8.9 21S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (6/4);  new base C (2/9);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCA  >  KT2440/5689692‑5689946
                                                                                                                            |                                                                                                                                  
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCGTCAGCCCGCTCCCCCCCGccccggg                                                                                                                          >  4:99685/1‑131 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGGGGCCTTAAGCCCGATCACCACCCCCCCGGgcgcgc                                                                                                                     >  8:269758/1‑135 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGGGCACATGACCCCGATCACCAACGCCCCTGgcgcgc                                                                                                                     >  3:139334/1‑135 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGGCGACATCAGCCCACTCACCACCCCCCCTGgcgcgc                                                                                                                     >  1:50204/1‑135 (MQ=255)
             cgCTGCCGCGGCCGGCCGGGCATCCCGGAGCGCGCCGAAACGGTGGCCCTGGCGCGCAACCTGCGGCTGGGGGAACAGCGCGGCGGGCGTGCGCACTTCACCCAGATGACCAGCGCCCGGGGCGCGCGGCTGATc                                                                                                             <  7:2449/134‑1 (MQ=255)
                   cccGGCCTGCCGGGCACCCCGTACCGCGCCGAAACCGGGCCCCTCGCGCGCAACCTGCGGCTGGGGCACCAGCGGGGCGGGCGGGCGCAGTCCCCCCAGAACACCCGCGCCCGGGGCGCGCGGCTGATCgcccag                                                                                                       <  2:170727/133‑1 (MQ=255)
                                                cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCCCGTCCGCCACCTCACCCCCGCCCCTGGCGCGCCGCGGCTCGCCCCCGCCCCCGCCGTGGGGCGGCCCGGAAccgc                                                                          >  8:33862/1‑135 (MQ=255)
                                                    aCCGGGGCCCGGGCGCGCAACCTGCTGCTGGTGCAACAGAGCGGCGTGCGTGCGCACTCCACCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                        <  7:269758/133‑1 (MQ=255)
                                                       ggggCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                   <  4:73870/133‑1 (MQ=255)
                                                       gTGGCCCTGGCGCGCAACATGCTGCTGGTGGAACAGCGCGGCGGGGGGGCGAACTTCAGCCAGAACCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                   <  2:50204/135‑1 (MQ=255)
                                                       gTGGCCCTGCCGCGCAACCGGCTGCCGGGGCCACAGAGCGGCGGGCGTGCGCACTTAAGCCATACCACCAGCGCCGGTGGCGCGCGTCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                   <  1:20619/135‑1 (MQ=255)
                                                            cgtGGGGCGCACCCTGGGGGTGGTGCACCCGCGCGGGGGGCGGGCGCACTGCCCCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCt                                                              <  1:235883/133‑1 (MQ=255)
                                                                                   tgtaaCAGCGCGGGGGGCTGGCGCACTTCCGCCTGATCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACg                                       <  7:227820/132‑1 (MQ=255)
                                                                                             cggcggGCGTGCGCACTGCAGCCAGAACCCCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                             <  4:259153/135‑1 (MQ=255)
                                                                                                  tgcgtgcgCACTTCACCCATAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                        <  4:46421/135‑1 (MQ=255)
                                                                                                  tgcgggcgCACTTCCCCCAGAACACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                        <  2:11372/135‑1 (MQ=255)
                                                                                                  ggcgtgcgCACTTCAGCCACAACCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                        <  8:272633/134‑1 (MQ=255)
                                                                                                       gcgcACTTCCGCCACACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                   <  4:278152/135‑1 (MQ=255)
                                                                                                             ttCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGGTCTCCAGCCTg             >  1:196058/1‑135 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCa  >  6:288282/1‑135 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCGTGGGCCTGCCGGTAACCGCCGCTGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCa  >  4:252854/1‑135 (MQ=255)
                                                                                                                            |                                                                                                                                  
AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCA  >  KT2440/5689692‑5689946

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: