Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,689,824 T→G 59.6% R238R (CGT→CGG pyrC' → Dihydroorotase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,689,8240TG59.6% 0.7 / 4.8 27R238R (CGT→CGGpyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base G (2/14);  total (6/21)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCG  >  KT2440/5689698‑5689950
                                                                                                                              |                                                                                                                              
cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCCCGTCAGCCCGCTCCCCACCGCCCGGGGcgcgcg                                                                                                                        >  6:171568/1‑135 (MQ=255)
       aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGGGCGCGGAAGCCACAACACCAGCGCCCCGGGcggccggcgggta                                                                                                                 >  1:143302/1‑130 (MQ=255)
         gttCCTCGGCCGGCCGGGCAGCCCGGAGCGCCCCCAAACGGGGGCCGGGGGGCGCAACGGGCGGCTGGGGCAACCGCGCGGCGTGCGGGCGCACTTCCGCCTCACCACCCGCGCCGGGGGCGCGCGGCTGATCgc                                                                                                               <  4:155236/134‑1 (MQ=255)
            ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGGGCCCTTCAGCCCCCTCACCACCGCCCCTGGCGCGCGGCGGATCgcccc                                                                                                            >  7:257325/1‑134 (MQ=255)
              tCGGCCTGCCGGGAATCCCGCATAGCGCCGAAACGGTGGCCCTGGCGCGAAACCTGCTGCGGTTGTAACCGAGGGGGGTGCGTGCGCACTTCACCCCTACCCCCCGCGCCGGGGGCGCGCGGCTGATCgcccagg                                                                                                          <  6:123148/135‑1 (MQ=255)
              tCGGCCGGCCGGGCATCCCGGAGAGCGCCGAAACCGGGGCCGGGGCGCGCAACCTGCTGCTGGTGGAACAGCGCGGCGGGCGTGCGCACTGCAGCCAGAGCACCCGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                          <  8:257325/135‑1 (MQ=255)
                                 gtagcgcgcGGAAACCGTGGCCCGGGCGCGCAACCTGCTGCGGTTGCAACAGCGCGGCGGGCGGGCTCACTTCAGCCATACCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGcc                                                                                       <  2:250653/130‑1 (MQ=255)
                                       gccctaaaCCGGGGCCCTGGCGGGAAACTTGTGGGGGTTGTAACTGCGGGGGGTGCGGGCGCATTGCCGCCATATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                 <  8:293965/130‑1 (MQ=255)
                                       ccccGAAACCGTGGCCCGGGCGCGAAACGTGCGGCGGGGGTAACAGCGCGGGGGGCGTGCGCACTACAGCCACACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                 <  2:51756/133‑1 (MQ=255)
                                              aCCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTAATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                            <  3:288712/133‑1 (MQ=255)
                                                 ggggCCCGGGCGCGCAACCTGCTGCGGGTGGAACAGAGGGGCGTGCGGGCGCACTCCCGCAAGAGCCCCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                       <  3:157743/133‑1 (MQ=255)
                                                 ggggCCCGGGCGCGCAACCTGCGGCGGGGGCACCCGCGGGGCGGGCGGGCGCACTCCCCCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                       <  1:244788/133‑1 (MQ=255)
                                                 ggggCCCGGGCGCGCAACCGGCTGGTGGTGCAACGGCGGGGCGGGCGTGCGCAGTACCGCCACACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                       <  2:143302/133‑1 (MQ=255)
                                                 gTGGCCGGGGCGCGCAACCTGCTGCTGGTGCAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                       <  3:118119/135‑1 (MQ=255)
                                                                             gggAACAGCGGGGCGGGCGGGCGCACTTCCGCCAGATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACg                                           <  8:25082/134‑1 (MQ=255)
                                                                               gAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATAACCCGCCCCCGTGGCGCGCGGCTGATCGCCCCGGCCCCGGCCGTGGGCCGGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCGCACGGACgcg                                         >  6:169056/1‑133 (MQ=255)
                                                                                     cgcggggGGCGGGCGCACTCCCCCCAGACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg                                   <  2:277482/134‑1 (MQ=255)
                                                                                            ggcgtgccCACTTCCGCCATTACCCCTGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                            <  5:74575/134‑1 (MQ=255)
                                                                                            ggcgggcgCCCTCCAGCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                            <  2:283392/134‑1 (MQ=255)
                                                                                              cgggcgcACTTCCGCCCGCTCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                          <  7:150747/132‑1 (MQ=255)
                                                                                                 gggcAGTTCCGCCACATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                       <  4:332407/133‑1 (MQ=255)
                                                                                                  cccACTTCAGCCATCCCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCa                      <  1:305167/133‑1 (MQ=255)
                                                                                                        gaaGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGt                <  8:285834/133‑1 (MQ=255)
                                                                                                                 gcccccgcgcCCGGGGCGCGCGGCGGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGc       <  1:213437/129‑1 (MQ=255)
                                                                                                                  cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTACCCGCCGCTGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTCCCCCGTGCa      >  2:141010/1‑135 (MQ=255)
                                                                                                                  cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCCCCGATCTAGCGCTGTACCAGCTGATCCTCACGGACGACTCGCTGCGCGAGTGCTCCCGCCTGTCCCACGTGCa      >  6:248962/1‑135 (MQ=255)
                                                                                                                      ggcgcCCGGGGCGCGCGGCAGATCGCCCAGGCCCAGGCCGTGGGCCTGCCGGGCACCGCCGACCTAGCGCTGTCCCAGCTGATCCTCCCGGACGAGGCGCTGCGCGAGCTCTCCAGCCTGTACCACGTGCAgccg  <  3:304008/134‑1 (MQ=255)
                                                                                                                              |                                                                                                                              
CATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCG  >  KT2440/5689698‑5689950

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