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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,689,7940TG46.4% 3.6 / 5.4 25R228R (CGT→CGGpyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base G (4/8);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  KT2440/5689661‑5689926
                                                                                                                                      |                                                                                                                                    
gcgACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGcgtgcgtg                                                                                                                                     <  5:362921/135‑1 (MQ=255)
    ccTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCCACCTGCTGCTGGTGGAACCGAGCGGcgtgggtggggt                                                                                                                                 >  2:370029/1‑131 (MQ=255)
     cTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTgcgcccg                                                                                                                               >  8:326956/1‑133 (MQ=255)
                 ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGGCCGTCAGcccccaac                                                                                                                   >  8:189036/1‑130 (MQ=255)
                     tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAg                                                                                                               >  4:88127/1‑135 (MQ=255)
                               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)
                                                                                                                                   gcgtgcgCACTTCAGCCAGAGCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCATGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt   <  2:62820/135‑1 (MQ=255)
                                                                                                                                   gcgggcgCACCACAGCCATACCCCCCGCGCCGGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc  <  8:205594/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  KT2440/5689661‑5689926

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: