Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,586,078 A→C 34.9% E230D (GAA→GAC)  accA1 → Acetyl‑/propionyl‑coenzyme A carboxylase alphachain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,586,0780AC34.9% 26.2 / 6.2 35E230D (GAA→GAC) accA1Acetyl‑/propionyl‑coenzyme A carboxylase alphachain
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base C (3/10);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 8.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTGGGTG  >  KT2440/4585946‑4586211
                                                                                                                                     |                                                                                                                                     
cTCGGCCCTGCGGGCGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGGGGAAAAGTCCGTGCTCACGCCCCGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGAAACTGCCTGTCCCT‑Caacgaccg                                                                                                                                     <  6:63789/135‑1 (MQ=255)
cTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑Caacgaacg                                                                                                                                     <  6:384587/135‑1 (MQ=255)
           gTGAGGCACAGTCGTGGTTCGGCGATGCGCGCATGCTGGTGGAAAAGGACGTGCTCCAGCCCCGCCCCGTGGAGCTTCAGGTGTTAGCCGAGCCGCCCGGCGACTGCCTGTACCT‑CGACGACCGCGACTGTTCCa                                                                                                                          >  5:255821/1‑135 (MQ=255)
               ggCACAGGGGTCGTTCGGCGAGGCGCGCATGGGGGTGGAAAAGTCCGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCC‑CACCGACCGCGACTGTTCCATTCa                                                                                                                      <  2:323789/135‑1 (MQ=255)
               gccactGTCGTGGTGCGGCGTGGCGCGCCGGCTGGTGGAAAAGTACGTGATCAAGCCACGCCAGGTGGAGATTCCGGGGTGCGCCGAGCACCACGCCAACTGCCGGTACCT‑CAACGACCGCGACTGTTCCATTCa                                                                                                                      <  4:125566/133‑1 (MQ=255)
                  actggcgtcgtGCGGCGAGGCGCGCATGCGGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATGCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑AAACGAACGCGACTGTTCCATTCAgcg                                                                                                                   <  5:199245/130‑1 (MQ=255)
                  acaGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAgcg                                                                                                                   >  7:38819/1‑135 (MQ=255)
                  aaaggcgtcgtGGGGCGATGGGCGCAGGCTGGTGTAAACGTCAGTGCTCCCCCCAGGCAACGTGGATATTCTGGTGGTCGCCTAGACCCACGGAACGTGCGGGTACAT‑CAACCAACGCGACTGTTCCATTCAgcg                                                                                                                   <  6:188903/130‑1 (MQ=255)
                       cgtcgtTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGACCGCGACTGTTCCATTCAGCgccgcc                                                                                                              >  3:54823/1‑135 (MQ=255)
                       cgtcgtTCGGCGATGCGCGCATGCGGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CCCCGCCCGCGCCTGTTGCATTCAGcgccg                                                                                                                >  3:333058/1‑133 (MQ=255)
                           gtggggCGATGCGCGCATGCGGTTGTAAAAGTACGGGCTCACGCCACGCCACGTGTAGATTCGGGTGTGCGCCGATCAGCACGGCAACTGCCTGAACCT‑CACCCACCGCGACTGTTCCATTCAGCGCCGCCATCa                                                                                                          <  8:210100/131‑1 (MQ=255)
                             tCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCa                                                                                                          >  8:208020/1‑133 (MQ=255)
                                gCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAgg                                                                                                     <  3:326215/135‑1 (MQ=255)
                                              cTGGTGTAAAAGAACGTGCTCAAGCCACGCCACGTGGATATTCAGGGGTGCGCCGATCACCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGc                                                                                        <  3:88020/134‑1 (MQ=255)
                                                    gAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGcc                                                                                 <  8:38819/135‑1 (MQ=255)
                                                         gTACGTGCTCAAGCCACCCCCCGTGGAGATTCAGGTTTTCGCCTATCAGCACGGCAACGGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCCGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCgg                                                                            <  1:237176/135‑1 (MQ=255)
                                                         gTACGTGCCCACGCCAGGCCACGTGTAGATTCAGGTGTGCGCCGATCAGCACGCAAATTGCGTGAACCA‑CACCGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCgg                                                                            <  6:103210/135‑1 (MQ=255)
                                                            cGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTCGAAGAGGCCCCTGCCCCCGGCCt                                                                         >  6:360721/1‑135 (MQ=255)
                                                              tGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGt                                                                       >  7:322594/1‑135 (MQ=255)
                                                                  cAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGcc                                                                   >  1:167022/1‑135 (MQ=255)
                                                                   acgccacgccagGGGGAGATTCAGGGGTGCGCCTACCCGCACGCCACCCGCCTGTACCT‑CACCCACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGccc                                                                  <  4:324418/133‑1 (MQ=255)
                                                                   accccacgccaggTGGAGATTCAGGGGTGGGCCGATCAGCAGGGCAACTGCGGGTACCT‑CACCGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGccc                                                                  <  6:98566/132‑1 (MQ=255)
                                                                       aacgccacgTGTAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGCAc                                                              >  2:60096/2‑135 (MQ=255)
                                                                             acgTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGCACGCGACTGTTCGATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGACCTGCGcc                                                        >  4:290282/1‑135 (MQ=255)
                                                                                                          cccGGCACCGGCCGGTACCT‑CACCCACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGGGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTgcc                           <  2:393995/133‑1 (MQ=255)
                                                                                                           agggCAACTGCCTGTACCT‑CAACGACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTgccg                          <  4:54823/133‑1 (MQ=255)
                                                                                                                 aatGCCTGTAACTACAACAACCGCTACTTTCCAATTCTGCGCCGCATTCATACGGTGTTGCAAAAGGCCCCT‑CCCCCGGCCTGTCGCACGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAgg                    <  2:296234/133‑1 (MQ=255)
                                                                                                                  cTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGc                   >  7:174057/1‑135 (MQ=255)
                                                                                                                      cTGAACCT‑CAACCACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATc               <  6:151988/135‑1 (MQ=255)
                                                                                                                       tGTACCT‑CAACGAAGGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCg              >  3:147494/1‑135 (MQ=255)
                                                                                                                       tGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCg              >  2:134490/1‑135 (MQ=255)
                                                                                                                           ccT‑CAGCGAACGCGAGTGTTCCATTCAGCGCCTCAATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGTCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTa          <  8:330635/135‑1 (MQ=255)
                                                                                                                           ccT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTgccgc                         >  5:294501/1‑120 (MQ=255)
                                                                                                                           ccT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTgccgc                         <  6:294501/120‑1 (MQ=255)
                                                                                                                            cT‑CAAGGACCGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTAc         <  1:19471/135‑1 (MQ=255)
                                                                                                                            cT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTAc         >  4:82971/1‑135 (MQ=255)
                                                                                                                            cT‑CAACCAACGCGACTGTTCCTTTCAGCGCCTCCATCAGAAGGTGGTGGAAGATGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTa          <  3:268730/134‑1 (MQ=255)
                                                                                                                             cacaccgaccgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTAc         <  7:314510/130‑1 (MQ=255)
                                                                                                                                  cgaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTggg    >  4:285735/1‑135 (MQ=255)
                                                                                                                                   gaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCgg                                                      >  3:380219/1‑84 (MQ=255)
                                                                                                                                   gaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCgg                                                      <  4:380219/84‑1 (MQ=255)
                                                                                                                                   gaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGCAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTGGGt   >  3:201157/1‑135 (MQ=255)
                                                                                                                                   gaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGTCTATCGGTTACGTGGGt   >  7:330635/1‑135 (MQ=255)
                                                                                                                                   gaacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTGGGt   >  2:140091/1‑135 (MQ=255)
                                                                                                                                    aacgCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTGGGTg  <  7:208020/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCT‑CAACGAACGCGACTGTTCCATTCAGCGCCGCCATCAGAAGGTGGTGGAAGAGGCCCCTGCCCCCGGCCTGTCGCCCGAACTGCGCCGGGCCATGGGCGAAGCGGCGGTGCGTGCCGCCCAGGCTATCGGTTACGTGGGTG  >  KT2440/4585946‑4586211

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