Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,689,824 T→G 45.1% 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,8240TG45.1% 12.6 / 5.7 32R238R (CGT→CGG) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base G (4/11);  total (11/21)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCT  >  KT2440/5689692‑5689955
                                                                                                                                    |                                                                                                                                   
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGCGCCCGTGCGCCCCATCACCACCCCCCCTgg                                                                                                                                   >  6:318894/1‑135 (MQ=255)
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCACATCACCAGCGCCCGTgg                                                                                                                                   >  4:254717/1‑135 (MQ=255)
  gggCCAGGCCCTGCTTCGGCGGCCTGCGGGGCATCCCGGAGAGCGCCGAAACCTTGGCCCGGGCGCGAAACCGGCTGCTGGGTGAACCGAGGGGCGTGCGTGCGCAATTCAGCCAGTTCACCTGCGCCGGGGgcg                                                                                                                                 <  8:389258/133‑1 (MQ=255)
        tgaccAGCTGCCTCGGCCGGCGGGGAAGCCCGGAGAGCGCCTAAACCGTGGCCCGGGCGCGCAACCTGCTGCTGGTGCAACAGAGCGGCGTGCGGGCGCACTTCACCCAGATCACCAGCGCCCGGGGCGCGCGGc                                                                                                                           <  8:360956/132‑1 (MQ=255)
           ccAGCTTCCGCGGCCTGCCGGGCACCCCGGAGAGCGCCGAACCCGGGGCCCGGCCGCGCAACCTGCTGCTGGTGCAACAGAGCGGCGGGCGTGCGCACTTCCGCCAGAGCCCCCGCGCCCGGGGCGCGCGGCt                                                                                                                          <  1:66459/133‑1 (MQ=255)
                  ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCCCTTCACCCACATCACCAGCCCCCCGGGCGCGCGGCGGCTCgcccc                                                                                                                 >  7:383690/1‑134 (MQ=255)
                  ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGGGGGGGGCCGTCAGCCCGCTCACCCACGCCCGGGGGGCGCGGGGGGTCgcccc                                                                                                                 >  5:126251/1‑134 (MQ=255)
                    gcggccggcCGGGCATCCCGGAGAGCGCCGAACCCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTGCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                               <  8:295123/134‑1 (MQ=255)
                    ccggccggcCGGGCACCCCGGCGCCCCCCTAAACCGGGGCCGTGGCGCCCACCCGGCTGGGGGGGCACCAGCGCGGCGGGCGGGCGAACTCCAGCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                               <  4:279543/134‑1 (MQ=255)
                         cTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGACCACCAGCGCCCGTGGCGCGCGGCTGATCgcccaggcccag                                                                                                          <  3:316518/135‑1 (MQ=255)
                           gCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGGGCCCTTAAGGCCCCACAACACCGCCCGGTGCGCCGGGGGGATCgcccacgcccccgc                                                                                                        >  1:317526/1‑131 (MQ=255)
                                             cgcCGAAAGGGGGGCCGTGGCGGGCAACGGGCGGGTGGGGCAACTCCGCGGGGGGCTGGCGCAGTACCCCCTGAACCCCAGCGCCGGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                      <  5:254794/135‑1 (MQ=255)
                                                  accccGTGGCCCGGCCGCGCACCCTGATCCGGGGGGAACAGAGGGGCGGGCGTGGGATATTCCGCCCGCGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCGCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                                 <  8:326820/132‑1 (MQ=255)
                                                    aCCGTGGCCGTGGCGCGCAACGTGCTGGTGGTGCAACTGCGCGGCGGGCGGGGGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                                 <  1:36502/133‑1 (MQ=255)
                                                       gTGGCCCTGGGGCGCAACCTGGGGCTGGTGTAACTGTGGGGCGGGCGTGGGCAGTCCTCCCTCAACACCCCGGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                            <  4:128558/135‑1 (MQ=255)
                                                       gTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                            <  3:194278/135‑1 (MQ=255)
                                                       gTGGCCCTGGCGCGCAACCTGCTGCGGGTGGACCAGAGCGGCGGGCGTGCGAACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                            <  2:13697/135‑1 (MQ=255)
                                                           cccTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGCCCATACCCCCGATCCCCCGCGCCCCGGGCGCGCGGCGGGTCGCCCCCGCCCCCGCCGTGGGCCGGGCCGGAACCCCCGCTGTAgc                                                                          >  3:354595/1‑133 (MQ=255)
                                                            ccTGGCGCGCAACCTGCTGCTGGTGCAACAGCGCGGCGTGCGTGCGCACTTCAGCCAGACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCt                                                                       <  3:254717/135‑1 (MQ=255)
                                                                     cAACCTGGTGCTGGTGCAACAGAGGGGCGTGCGTGGGAATTTACGCCATAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                              <  8:310938/135‑1 (MQ=255)
                                                                                             cggcggGCGTGCGCACTTCCGCCCTATCCCCAGGGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                      <  8:269116/135‑1 (MQ=255)
                                                                                                  ggggggcgcACTGCCCCCAGACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                 <  4:335666/130‑1 (MQ=255)
                                                                                                    cgtgcgCACTGCAGCCCGAGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                               <  8:363758/135‑1 (MQ=255)
                                                                                                       gcgcACTCCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                            <  4:103457/135‑1 (MQ=255)
                                                                                                        cgcACATCAGCCAGAACACCAGCGCCCCTGGCGCGCCGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCCGTAACCGCCGATGTACCCCTGTACCAGCTGATCCTAACCGAAGAGTCCCTGCGCGCGTTATCCa                           >  8:284350/1‑135 (MQ=255)
                                                                                                            cTTCAGCCAGATCCCCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCt                       <  8:386253/135‑1 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCa           >  4:215746/1‑135 (MQ=255)
                                                                                                                          ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgc         >  2:177765/1‑135 (MQ=255)
                                                                                                                             gcgcCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                 <  6:365811/108‑1 (MQ=255)
                                                                                                                             gcgcCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                 >  5:365811/1‑108 (MQ=255)
                                                                                                                                 ccGGGGCCCGCGGCTGAGCGCCCCGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGGGCAGCCGCCGCt  <  1:341500/135‑1 (MQ=255)
                                                                                                                                   gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                               >  1:246492/1‑104 (MQ=255)
                                                                                                                                    |                                                                                                                                   
AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCT  >  KT2440/5689692‑5689955

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