Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,662,000 T→A 8.5% H265Q (CAT→CAA)  mmgE → Citrate/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,662,0000TA 8.5% 136.2 / 4.2 47H265Q (CAT→CAA) mmgECitrate/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (17/26);  new base A (2/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AAGTCCTGGGCTGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTCATGCGCAGACCGC  >  KT2440/2661869‑2662131
                                                                                                                                   |                                                                                                                                   
aaGTCCTGGGCTGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATAc                                                                                                                                   >  1:228005/1‑134 (MQ=255)
         gCTGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAg                                                                                                                           <  4:262891/133‑1 (MQ=255)
           tGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCCCGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCt                                                                                                                       >  7:281093/1‑135 (MQ=255)
           tGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCt                                                                                                                       >  1:70804/1‑135 (MQ=255)
              gggTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTg                                                                                                                      <  1:44995/133‑1 (MQ=255)
                gTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGTGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGgcgc                                                                                                                  >  8:170025/1‑135 (MQ=255)
                gTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCAATACCAACAAAGACCTGgcgc                                                                                                                  >  7:73015/1‑135 (MQ=255)
                gTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGCGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCCTTCAGCCATACCAACAAAGACCTGgcg                                                                                                                   >  5:36960/1‑134 (MQ=255)
                 tGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCt                                                                                                                 <  8:287346/135‑1 (MQ=255)
                 tGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCt                                                                                                                 <  7:83349/135‑1 (MQ=255)
                      ccTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCGTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGc                                                                                                            <  2:185684/135‑1 (MQ=255)
                       cTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGcc                                                                                                           >  6:35007/1‑135 (MQ=255)
                           aGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCt                                                                                                       <  1:315499/135‑1 (MQ=255)
                                   tgtgCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTa                                                                                               >  3:135871/1‑135 (MQ=255)
                                   tgtgCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTCGCGCTCAAGCCCGCTGGGCAGTa                                                                                               >  7:192701/1‑135 (MQ=255)
                                    gtgCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTa                                                                                               >  4:122272/1‑134 (MQ=255)
                                          cTGGCCGACATCGCCGTGCGGGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGGCGCAGTGGGGGGTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTg                                                                                        <  7:287240/135‑1 (MQ=255)
                                           tGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTgc                                                                                       <  2:100652/135‑1 (MQ=255)
                                               cGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCt                                                                                   >  1:39826/1‑135 (MQ=255)
                                                           gagTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCAAACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCggt                                                                       >  1:336755/3‑133 (MQ=255)
                                                           gCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCAAACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       >  1:336757/1‑135 (MQ=255)
                                                           gCGTGGGGAAATGGGCGGGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       <  7:161599/135‑1 (MQ=255)
                                                            cGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCAACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       <  3:269566/134‑1 (MQ=255)
                                                            cGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGGGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       <  2:32081/134‑1 (MQ=255)
                                                              tGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTggg                                                                    <  8:21350/135‑1 (MQ=255)
                                                                gggAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCa                                                                  <  6:154146/135‑1 (MQ=255)
                                                                gggAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCCTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCa                                                                  <  4:259591/135‑1 (MQ=255)
                                                                   aaaTGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGAc                                                                                                                         >  3:266570/1‑77 (MQ=255)
                                                                   aaaTGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGAc                                                                                                                         <  4:266570/77‑1 (MQ=255)
                                                                       gggCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTAtgt                                                           <  5:35007/135‑1 (MQ=255)
                                                                              ccGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGga                                                     >  1:3427/1‑134 (MQ=255)
                                                                                 ggCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAAc                                                 <  1:51523/135‑1 (MQ=255)
                                                                                    gTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGGTGGGCAGCTATGTGATGGAGAACGTg                                              <  8:120099/135‑1 (MQ=255)
                                                                                               gTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGt                                   <  2:228005/135‑1 (MQ=255)
                                                                                                   cAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAg                                <  5:153989/134‑1 (MQ=255)
                                                                                                    aGTGGGGGTTCTACGATGTTTCGTTCAGCCAAACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGtt                              <  2:336757/135‑1 (MQ=255)
                                                                                                    aGTGGGGGTTCTACGATGTTTCGTTCAGCCAAACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGtt                              <  2:336755/135‑1 (MQ=255)
                                                                                                         gggTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAAc                                                 <  7:188653/111‑1 (MQ=255)
                                                                                                         gggTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAAc                                                 >  8:188653/1‑111 (MQ=255)
                                                                                                         gggTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCt                         <  5:209878/135‑1 (MQ=255)
                                                                                                               tACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAg                   <  1:333921/135‑1 (MQ=255)
                                                                                                                   aTGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTc               >  1:223550/1‑135 (MQ=255)
                                                                                                                       ttCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       <  2:151687/75‑1 (MQ=255)
                                                                                                                       ttCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGtt                                                                       >  1:151687/1‑75 (MQ=255)
                                                                                                                           ttCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTCATGCGCAg       >  1:28606/1‑135 (MQ=255)
                                                                                                                              aGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTCATGCGCAGAcc    <  2:9335/135‑1 (MQ=255)
                                                                                                                                ccATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTCATGCGCAGACCGc  <  8:17561/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
AAGTCCTGGGCTGCGGGTGACGCCTCCAGCCGAGGTGTGCGCCTGGCCGACATCGCCTTGCGTGGGGAAATGGGCGTGCCGGGCGTGCTCACGGCGTCGCAGTGGGGGTTCTACGATGTTTCGTTCAGCCATACCAACAAAGACCTGGCGCTCAAGCCCGCTGGGCAGTACGAGCTGCGCCTGCCCCAGGCGTTGGGCAGCTATGTGATGGAGAACGTGCTGTTCAAGGTCAGTTTCCCTGCCGAGTTTCATGCGCAGACCGC  >  KT2440/2661869‑2662131

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