Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,773,653 C→T 16.3% intergenic (+42/‑141) aceA → / → aceK isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,773,6530CT16.3% 136.2 / 16.1 61intergenic (+42/‑141)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base T (4/6);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCC  >  Exported/2773583‑2773721
                                                                      |                                                                      
ccGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTc                                                                        <  1:4093058/71‑1 (MQ=255)
 cGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTc                                                                        <  1:4405078/70‑1 (MQ=255)
  ggCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAc                                                                      >  1:2033487/1‑71 (MQ=255)
      ccACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgccg                                                                  >  1:2716499/1‑71 (MQ=255)
        aCTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                <  1:3086690/71‑1 (MQ=255)
         cTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgcc                                                                   >  1:4812974/1‑67 (MQ=255)
           gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCa                                                                       >  1:861432/1‑61 (MQ=255)
           gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgc                                                                    >  1:4796055/1‑64 (MQ=255)
           gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCAcgc                                                                    >  1:3336804/1‑64 (MQ=255)
           gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                              >  1:1321008/1‑70 (MQ=255)
           gaagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTACAcgc                                                                    >  1:5198874/1‑64 (MQ=255)
            aagaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                            >  1:2625098/1‑71 (MQ=255)
                     aGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCgg                                                    <  1:2431421/70‑1 (MQ=255)
                         cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCa                                                                       >  1:1625255/1‑47 (MQ=255)
                         cTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCa                                                                       >  1:2425950/1‑47 (MQ=255)
                         cTAAGCAACAACAACAGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCAc                                               >  1:2642731/1‑71 (MQ=255)
                             gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                          >  1:2082521/1‑71 (MQ=255)
                             gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                          >  1:2750278/1‑71 (MQ=255)
                             gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                          >  1:1488604/1‑71 (MQ=255)
                             gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                          >  1:4513233/1‑71 (MQ=255)
                             gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                          >  1:4473628/1‑71 (MQ=255)
                               aacaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                <  1:1965872/48‑1 (MQ=255)
                                acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                <  1:1562522/47‑1 (MQ=255)
                                acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCAc                                               <  1:651087/64‑1 (MQ=255)
                                acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCt                                       <  1:551396/71‑1 (MQ=255)
                                 caacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCt                                       >  1:5141301/1‑70 (MQ=255)
                                 caacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTg                                      <  1:1795518/71‑1 (MQ=255)
                                  aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGa                                     >  1:2894040/1‑71 (MQ=255)
                                  aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGa                                     >  1:3714733/1‑71 (MQ=255)
                                   acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGAt                                    >  1:499976/1‑71 (MQ=255)
                                             cTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                       <  1:2311645/43‑1 (MQ=255)
                                               gACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTgcg                       >  1:5337156/1‑71 (MQ=255)
                                                 cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGcgcg                     >  1:4337777/1‑71 (MQ=255)
                                                 cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGcgcg                     <  1:310363/71‑1 (MQ=255)
                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCt                  >  1:4861730/1‑70 (MQ=255)
                                                      ggCCGGATAAGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTa                <  1:2980922/70‑1 (MQ=11)
                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGc                                  <  1:1754404/54‑1 (MQ=255)
                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCtt                 <  1:702953/70‑1 (MQ=255)
                                                      gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCtt                 <  1:3613819/68‑1 (MQ=1)
                                                      gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTa                <  1:895741/69‑1 (MQ=11)
                                                      gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTa                <  1:3007464/69‑1 (MQ=11)
                                                      gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTa                <  1:3307294/69‑1 (MQ=11)
                                                      agCCGGATAAGGCGCTCGCGCCGCATCCGGCAA‑CGGTGCCAGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt               <  1:1428683/70‑1 (MQ=14)
                                                       gCCGGATAAGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTAt               <  1:4359086/70‑1 (MQ=11)
                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt               <  1:1913176/71‑1 (MQ=255)
                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt               <  1:2817935/71‑1 (MQ=255)
                                                       accGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTAt               <  1:1166876/70‑1 (MQ=11)
                                                        ccGGATAAGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCa             <  1:3793475/71‑1 (MQ=11)
                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTATc              <  1:4271156/71‑1 (MQ=11)
                                                         cGGATAAGGCGCTCGCGCCGCATCCGGCAA‑CGGTGCCAGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCag            >  1:1015364/1‑70 (MQ=14)
                                                          ggATAAGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTAt               <  1:5187437/67‑1 (MQ=11)
                                                          ggATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTATCag            >  1:455277/1‑70 (MQ=1)
                                                            aTAAGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGcc         >  1:3814546/1‑71 (MQ=11)
                                                              aaGGCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTa       >  1:234268/1‑71 (MQ=11)
                                                               aGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCt        <  1:3737784/70‑1 (MQ=255)
                                                                ggCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTa       >  1:3284257/1‑69 (MQ=11)
                                                                ggCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTACa     <  1:3713134/71‑1 (MQ=11)
                                                                ggCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTACa     <  1:169435/71‑1 (MQ=11)
                                                                ggCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTAc      <  1:423264/71‑1 (MQ=255)
                                                                  cGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTAt               >  1:4935700/1‑59 (MQ=11)
                                                                    ttCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGc                                  <  1:214201/40‑1 (MQ=40)
                                                                    ttCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGcc  <  1:5093041/71‑1 (MQ=255)
                                                                    ttCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGcc  >  1:4729389/1‑71 (MQ=255)
                                                                      |                                                                      
CCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCC  >  Exported/2773583‑2773721

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