Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 109,695 C→T 14.7% H647Y (CAC→TAC)  aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported109,6950CT14.7% 150.3 / 16.4 62H647Y (CAC→TAC) aceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (19/33);  new base T (7/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAT  >  Exported/109626‑109765
                                                                     |                                                                      
cTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCca                                                                       >  1:4869140/1‑71 (MQ=255)
cTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCca                                                                       >  1:4986175/1‑71 (MQ=255)
cTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCca                                                                       >  1:3748247/1‑71 (MQ=255)
  gATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                                     >  1:4488963/1‑71 (MQ=255)
  gATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                                     >  1:3778329/1‑71 (MQ=255)
    tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGGCACAtt                                                                   <  1:1969520/71‑1 (MQ=255)
    tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                   <  1:1128471/71‑1 (MQ=255)
    tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                   <  1:3435891/71‑1 (MQ=255)
    tCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                   <  1:1721085/71‑1 (MQ=255)
      ggcggTACTTCCGGTCGTACCACCCTGAACGTCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                                                  <  1:3333204/70‑1 (MQ=255)
      ggcggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                                                  <  1:106912/70‑1 (MQ=255)
      ggcggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                                                  <  1:3839412/70‑1 (MQ=255)
      ggcggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                                 <  1:4323920/71‑1 (MQ=255)
      ggcggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCAACATTCa                                                                 <  1:1727404/71‑1 (MQ=255)
             cTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                                           <  1:4076172/70‑1 (MQ=255)
             cTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                                           >  1:2036212/1‑70 (MQ=255)
             cTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                                           >  1:3036260/1‑70 (MQ=255)
             cTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                                           >  1:2544627/1‑70 (MQ=255)
                ccGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:2328172/70‑1 (MQ=255)
                     cGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATccc                                                  >  1:2794630/1‑71 (MQ=255)
                         ccaccCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                                 >  1:945725/1‑52 (MQ=255)
                         ccaccCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                                 >  1:698482/1‑52 (MQ=255)
                           accCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                                            <  1:2165794/71‑1 (MQ=255)
                            cccTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                                            <  1:4539368/70‑1 (MQ=255)
                            cccTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                                            <  1:2759283/70‑1 (MQ=255)
                                gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:2641005/54‑1 (MQ=255)
                                gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:2593176/54‑1 (MQ=255)
                                gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:250673/54‑1 (MQ=255)
                                gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                                        <  1:377055/54‑1 (MQ=255)
                                    ggCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                   <  1:4217739/39‑1 (MQ=255)
                                    ggCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                   <  1:687377/39‑1 (MQ=255)
                                    ggCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTa                                   <  1:4637397/71‑1 (MQ=255)
                                    ggCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTa                                   <  1:4667080/71‑1 (MQ=255)
                                     gCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                                  >  1:2618749/1‑71 (MQ=255)
                                      cGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                                 <  1:4663100/71‑1 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:1289223/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:3597082/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:3536865/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:2033612/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:168373/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:244504/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:2468025/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:5083203/1‑71 (MQ=255)
                                         aGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc                              >  1:3103786/1‑71 (MQ=255)
                                            tCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           <  1:1251081/71‑1 (MQ=255)
                                            tCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg                           <  1:3448589/71‑1 (MQ=255)
                                            tCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCg                            >  1:2030991/1‑70 (MQ=255)
                                              tGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATctct                                     <  1:2665599/59‑1 (MQ=255)
                                              tGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATctct                                     <  1:744967/59‑1 (MQ=255)
                                              tGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATctct                                     <  1:2889208/59‑1 (MQ=255)
                                               gcagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         >  1:3360042/1‑70 (MQ=255)
                                               gcagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         >  1:4530336/1‑70 (MQ=255)
                                               gcagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct                         >  1:970608/1‑70 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta                       <  1:1209154/71‑1 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCTACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta                       <  1:1657897/71‑1 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                                 <  1:4939755/61‑1 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                                 <  1:1101214/61‑1 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                                 <  1:375717/61‑1 (MQ=255)
                                                cagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACg                                 <  1:2789815/61‑1 (MQ=255)
                                                     cGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgct                   >  1:4720876/1‑70 (MQ=255)
                                                     cGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgct                   >  1:772327/1‑70 (MQ=255)
                                                       aaGATGGTCACAGACACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTa                                   <  1:4731453/52‑1 (MQ=255)
                                                                acaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:4494903/1‑71 (MQ=255)
                                                                acaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:4220570/1‑71 (MQ=255)
                                                                acaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:3644930/1‑71 (MQ=255)
                                                                acaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCt       >  1:2295975/1‑71 (MQ=255)
                                                                   gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc    <  1:99324/71‑1 (MQ=255)
                                                                    ccACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtca   >  1:1722807/1‑71 (MQ=255)
                                                                     cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtcat  <  1:1349106/71‑1 (MQ=255)
                                                                     |                                                                      
CTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAT  >  Exported/109626‑109765

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