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→CGG) pyrC'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

GCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  KT2440/5689661‑5689926
                                                                                                                                      |                                                                                                                                    
gcgACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGcgtgcgtg                                                                                                                                     <  5:362921/135‑1 (MQ=255)
    ccTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCCACCTGCTGCTGGTGGAACCGAGCGGcgtgggtggggt                                                                                                                                 >  2:370029/1‑131 (MQ=255)
     cTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTgcgcccg                                                                                                                               >  8:326956/1‑133 (MQ=255)
                 ggcCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGGCCGTCAGcccccaac                                                                                                                   >  8:189036/1‑130 (MQ=255)
                     tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCG‑CGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAg                                                                                                               >  4:88127/1‑135 (MQ=255)
                               aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACGTCAGCCACCTAACCAGCGCCCCTgg                                                                                                      >  2:261917/1‑135 (MQ=255)
                                    ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGGCCGTCAGCCAGCTAACCACCGCCCGGGGcgcgcg                                                                                                >  7:217826/1‑136 (MQ=255)
                                    ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACATCAGCCACATCACCACCGCCCCTGGcgcgcg                                                                                                >  8:73297/1‑136 (MQ=255)
                                    ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGGGCCCTTACGCCAGATCACCAGCCCCCGTGgcgcgc                                                                                                 >  6:171795/1‑135 (MQ=255)
                                       tGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCTCCCGTGGCGCGCGGc                                                                                              <  5:333042/135‑1 (MQ=255)
                                            aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTCCGCCCGCTCCCCAGCCCCCGGGGCGCGCCGCTGCTc                                                                                         >  2:121087/1‑135 (MQ=255)
                                            aGCGTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACGGTGGCCGTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCg                                                                                        <  7:326956/136‑1 (MQ=255)
                                                   tgggCGTGCCGGTCATCCCGGAGAGCGCCTAAACCGTGGCCCTGG‑CCCGAACCCTGCTGCGGGGGGACCTGAGCGGCGGGCGTGCGAACTTCAGCCAGACCACCCGCGCCGGTGGCGCGCGGCTGATCgcccagg                                                                                  <  6:325774/133‑1 (MQ=255)
                                                    cggccggcGGGGCATCCCGTAGAGCCCCGAAACCGTGGCCCTGG‑CGCGAAACCTGGTGCTGGTGTACCTGTGGGGGGTGCGTGCGAACTGCTCCCACATCACCCGCGCCGGGGGCGCGCGGCTGATCgcccaggc                                                                                 <  6:429458/135‑1 (MQ=255)
                                                        cTGCCGGGCATCCGGGAGCGCGCCTAACCCGGGCCCGTGG‑CGCGCAACCTGCGGGTGGGGGAACAGCGCGGCGTGCGGGCCCACTCCAGCCAGATCACCCGGGCCGGGGGCGCGCGGCTGATCgcccaggcccag                                                                             <  2:33762/135‑1 (MQ=255)
                                                          gCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGCGCCATTCCCCCCGATCACCACCGCCCGGGGCGCGCGGCTGTTAGCCcccgcccacgcc                                                                          >  8:19019/1‑136 (MQ=255)
                                                                           gcgcCGAACCCGTGGCCCGGG‑CGCGCAACAGGCGGCGGGGGGAACAGCGCGGCGTGCGGGCGCACTGCAGCCACATCACCCGGGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                         <  7:342818/136‑1 (MQ=255)
                                                                            cgcCGAAACCGGGGCCAGGG‑GGCGCAACCTGCTGCGGGGGTAACACAGCGGCTGGGGGGCGAACTTAACCCTAATCACCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                         <  6:414281/135‑1 (MQ=255)
                                                                                     ggTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                               <  7:236673/135‑1 (MQ=255)
                                                                                      gTGGCCCTGG‑CGCGCAACCTGCTGCTGGGGTACCAGCGGGGCGTGCGGGCGCACTGCCGCCAGATCCCCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                               <  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)
                                                                                                                                      |                                                                                                                                    
GCGACCTGGCCCAGGGCGGCCTGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGG‑CGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTC  >  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: ‑