Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 108,314 A→C 6.8% E186D (GAA→GAC aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported108,3140AC 6.8% 194.5 / 8.3 74E186D (GAA→GACaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base A (38/30);  new base C (2/3);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTCCTG  >  Exported/108245‑108383
                                                                     |                                                                     
ggATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      >  1:4479318/1‑71 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                                                                     <  1:5081467/71‑1 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                                                                     <  1:2364702/71‑1 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:2058886/70‑1 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:1249502/70‑1 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:3102709/70‑1 (MQ=255)
 gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:1688861/70‑1 (MQ=255)
  aTAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTc                                                                    >  1:4948167/1‑71 (MQ=255)
   tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                                                   >  1:1650229/1‑71 (MQ=255)
   tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                                                   >  1:3270176/1‑71 (MQ=255)
   tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                                                   >  1:1417434/1‑71 (MQ=255)
   gaaCTTCCGCCGCGAAGTCGACGGCAAGGGCCTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCt                                                                   >  1:1308823/2‑71 (MQ=255)
   gaaCTTCCGCCGCGAAGTCGACGGCAAGGGCCTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCt                                                                   >  1:1142796/2‑71 (MQ=255)
            tCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                                             <  1:2446882/68‑1 (MQ=255)
            tCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                                             <  1:2377000/68‑1 (MQ=255)
              aGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTc                                                           >  1:3796696/1‑68 (MQ=255)
                gAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:96439/55‑1 (MQ=255)
                gAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:5028764/55‑1 (MQ=255)
                   gTTCTCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                                                            <  1:803461/62‑1 (MQ=255)
                   gTTCACGGCAATGGCCTCTCTTCCTATCCGCCCCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                                                     >  1:3341770/1‑69 (MQ=255)
                   gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                                                            <  1:1681075/62‑1 (MQ=255)
                   gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                                                            <  1:3827363/62‑1 (MQ=255)
                   gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                                                            <  1:4049966/62‑1 (MQ=255)
                   gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                                                            <  1:4356506/62‑1 (MQ=255)
                     tCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTccc                                                         >  1:2118353/1‑63 (MQ=255)
                     tCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTccc                                                         >  1:544506/1‑63 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCg                                                    >  1:4633758/1‑67 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt                                                   >  1:4345289/1‑68 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt                                                   >  1:944058/1‑68 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                >  1:1731088/1‑71 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                >  1:1188865/1‑71 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                >  1:2543948/1‑71 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                >  1:1973453/1‑71 (MQ=255)
                      cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                >  1:3405465/1‑71 (MQ=255)
                       aCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                <  1:4580337/70‑1 (MQ=255)
                       aCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                                                <  1:4015569/70‑1 (MQ=255)
                         ggCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAt                                             >  1:2389253/1‑71 (MQ=255)
                         ggCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAt                                             >  1:901876/1‑71 (MQ=255)
                            aaTGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtggg                                          <  1:4713268/71‑1 (MQ=255)
                              tGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                                             <  1:2142976/50‑1 (MQ=255)
                                gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTcc                                                          >  1:4759201/1‑51 (MQ=255)
                                gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctg                                      >  1:2655520/1‑71 (MQ=255)
                                gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctg                                      >  1:3543206/1‑71 (MQ=255)
                                gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtct                                       <  1:2425248/70‑1 (MQ=255)
                                gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtct                                       <  1:4802160/70‑1 (MQ=255)
                                 ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg                                     <  1:928690/71‑1 (MQ=255)
                                 ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg                                     <  1:5077737/71‑1 (MQ=255)
                                 ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg                                     <  1:5996/71‑1 (MQ=255)
                                 ccTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                                                     <  1:1216028/55‑1 (MQ=255)
                                   tctctTACTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctgggt                                   >  1:3214/1‑71 (MQ=255)
                                    ctctTCCTATCCGCACCCGAAACTGATGCCGGAAt                                                                      <  1:3545755/35‑1 (MQ=255)
                                       ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctgggt                                   >  1:274237/1‑67 (MQ=255)
                                       ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTcc                                 <  1:684846/69‑1 (MQ=255)
                                       ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTcc                                 <  1:2824243/69‑1 (MQ=255)
                                       ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGa                               >  1:542396/1‑71 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGTCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTgg                           <  1:2490622/69‑1 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                                                     >  1:4302950/1‑43 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                                                     >  1:1997979/1‑43 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                                                     >  1:951107/1‑43 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTg                            >  1:1552038/1‑68 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTg                            >  1:3677007/1‑68 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGtt                         <  1:924481/71‑2 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTg                         <  1:3658707/71‑1 (MQ=255)
                                             tCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTg                         <  1:1874679/71‑1 (MQ=255)
                                               cGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTa                                              >  1:2477214/1‑48 (MQ=255)
                                                     cGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtc                                        >  1:1721947/1‑48 (MQ=255)
                                                     cGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTAc                 >  1:1730308/1‑71 (MQ=255)
                                                     cGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTAc                 >  1:1135992/1‑71 (MQ=255)
                                                              tGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTAtt                    <  1:892675/59‑1 (MQ=255)
                                                                  ggATTTCTGGGAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTgg                           <  1:759181/48‑1 (MQ=255)
                                                                  ggAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTcc    >  1:4338762/1‑71 (MQ=255)
                                                                  ggAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTcc    >  1:3798498/1‑71 (MQ=255)
                                                                  ggAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTcc    >  1:3618528/1‑71 (MQ=255)
                                                                  ggAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTcc    >  1:3165735/1‑71 (MQ=255)
                                                                    aaTTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTAcc                <  1:2197341/57‑1 (MQ=255)
                                                                    aaTTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTAcc                <  1:2844606/57‑1 (MQ=255)
                                                                    aaTTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTCCTg  <  1:1530168/71‑1 (MQ=255)
                                                                     aTTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTCCTg  >  1:516887/1‑70 (MQ=255)
                                                                     |                                                                     
GGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGATTGGTGCTATTTACCAGGCTAAATTCCTG  >  Exported/108245‑108383

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