Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,689,824 T→G 50.2% 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,8240TG50.2% 8.0 / 24.8 35R238R (CGT→CGGpyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base G (3/15);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 3.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGCG  >  KT2440/5689692‑5689958
                                                                                                                                     |                                                                                                                                      
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACGTCAGCCACCTAACCAGCGCCCCTgg                                                                                                                                      >  2:261917/1‑135 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGGCCGTCAGCCAGCTAACCACCGCCCGGGGcgcgcg                                                                                                                                >  7:217826/1‑136 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACATCAGCCACATCACCACCGCCCCTGGcgcgcg                                                                                                                                >  8:73297/1‑136 (MQ=255)
     ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGGGCCCTTACGCCAGATCACCAGCCCCCGTGgcgcgc                                                                                                                                 >  6:171795/1‑135 (MQ=255)
        tGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCTCCCGTGGCGCGCGGc                                                                                                                              <  5:333042/135‑1 (MQ=255)
             aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTCCGCCCGCTCCCCAGCCCCCGGGGCGCGCCGCTGCTc                                                                                                                         >  2:121087/1‑135 (MQ=255)
             aGCGTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACGGTGGCCGTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCg                                                                                                                        <  7:326956/136‑1 (MQ=255)
                    tgggCGTGCCGGTCATCCCGGAGAGCGCCTAAACCGTGGCCCTGGCCCGAACCCTGCTGCGGGGGGACCTGAGCGGCGGGCGTGCGAACTTCAGCCAGACCACCCGCGCCGGTGGCGCGCGGCTGATCgcccagg                                                                                                                  <  6:325774/133‑1 (MQ=255)
                     cggccggcGGGGCATCCCGTAGAGCCCCGAAACCGTGGCCCTGGCGCGAAACCTGGTGCTGGTGTACCTGTGGGGGGTGCGTGCGAACTGCTCCCACATCACCCGCGCCGGGGGCGCGCGGCTGATCgcccaggc                                                                                                                 <  6:429458/135‑1 (MQ=255)
                         cTGCCGGGCATCCGGGAGCGCGCCTAACCCGGGCCCGTGGCGCGCAACCTGCGGGTGGGGGAACAGCGCGGCGTGCGGGCCCACTCCAGCCAGATCACCCGGGCCGGGGGCGCGCGGCTGATCgcccaggcccag                                                                                                             <  2:33762/135‑1 (MQ=255)
                           gCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGCGCCATTCCCCCCGATCACCACCGCCCGGGGCGCGCGGCTGTTAGCCcccgcccacgcc                                                                                                          >  8:19019/1‑136 (MQ=255)
                                            gcgcCGAACCCGTGGCCCGGGCGCGCAACAGGCGGCGGGGGGAACAGCGCGGCGTGCGGGCGCACTGCAGCCACATCACCCGGGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                         <  7:342818/136‑1 (MQ=255)
                                             cgcCGAAACCGGGGCCAGGGGGCGCAACCTGCTGCGGGGGTAACACAGCGGCTGGGGGGCGAACTTAACCCTAATCACCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                         <  6:414281/135‑1 (MQ=255)
                                                      ggTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                               <  7:236673/135‑1 (MQ=255)
                                                       gTGGCCCTGGCGCGCAACCTGCTGCTGGGGTACCAGCGGGGCGTGCGGGCGCACTGCCGCCAGATCCCCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                               <  4:201994/135‑1 (MQ=255)
                                                        gggCCCCGGCCGCGCACCGGGCGGCGGGGGGACCCGAGGGGGGGGCGTGCGCACTTCCGCCAGACCCCCCCCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                               <  4:328696/134‑1 (MQ=255)
                                                                   gcgcACCCTGGGGCTGGGGGAACAGCGCGGGGGGGGGGCGCACTTCCCCCATACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAg                                                                   <  7:61725/136‑1 (MQ=255)
                                                                    cgcACCCTGCGGGGGGGGGAACAGAGGGGCGGGCGTGCGCACTTCAGCCAGAGCCCCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAg                                                                   <  4:57641/135‑1 (MQ=255)
                                                                                      gAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCACCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCCGCCTTGGGCCTTCCGGTAACCGCCGCTGTAGCGCTGTGACAGCTGATCCTCCCGGCCGAGt                                                >  7:288959/1‑136 (MQ=255)
                                                                                              cggggggcgggcgCCCTTCCGCCCGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                         <  2:378721/129‑1 (MQ=255)
                                                                                              cggcggGCGGGCGCACTCCCCCCAGACCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                         <  4:14850/135‑1 (MQ=255)
                                                                                              cggcgTGCGGGCGCACTTCAGCCAGAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                         <  4:267889/135‑1 (MQ=255)
                                                                                               ggcgTGGGTGCCCACTTCCGCCAGCTCCCCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcgc                                        <  4:161674/135‑1 (MQ=255)
                                                                                                    gggggcgcACTGCCCCCACACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                                   <  6:437551/131‑1 (MQ=255)
                                                                                                    gcgtgcgCACTTCAGCCAGAGCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCATGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                                   <  2:62820/135‑1 (MQ=255)
                                                                                                    gcgggcgCACCACAGCCATACCCCCCGCGCCGGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                                  <  8:205594/136‑1 (MQ=255)
                                                                                                        gcgcAGTGCAGCCAGAGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                               <  6:200718/135‑1 (MQ=255)
                                                                                                                 ccccAGATCCCCCGCGCCCGGGCCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGt                        <  6:82152/131‑1 (MQ=255)
                                                                                                                   ccAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTGCTCCAGCCTGTACCa                    >  5:254788/1‑135 (MQ=255)
                                                                                                                          ccccgggcCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg             <  2:246378/129‑1 (MQ=255)
                                                                                                                             cgcgcCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgccg          <  5:171795/134‑1 (MQ=255)
                                                                                                                                gcCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGccgccg       >  5:202499/1‑135 (MQ=255)
                                                                                                                                 cccGTGGCGCGCGGCTAATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCcgccgc      >  3:243752/1‑135 (MQ=255)
                                                                                                                                    gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACGAGCTGATCCTCACGGACGAGTCGCTGCGCGCGGTCTCCAGCCTGTACCACGTGCCGCCGCCGCTgcg  >  8:359969/1‑136 (MQ=255)
                                                                                                                                    gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGCGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTgc   >  2:180207/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                      
AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGCG  >  KT2440/5689692‑5689958

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: