Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,773,680:1 +C 15.9% intergenic (+69/‑114) 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,6801.C15.9% 200.7 / 51.0 69intergenic (+69/‑114)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base . (32/26);  new base C (4/7);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTT  >  Exported/2773612‑2773750
                                                                     |                                                                       
gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                                                       >  1:2750278/1‑71 (MQ=255)
gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                                                       >  1:4473628/1‑71 (MQ=255)
gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                                                       >  1:4513233/1‑71 (MQ=255)
gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                                                       >  1:1488604/1‑71 (MQ=255)
gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑Tg                                                                       >  1:2082521/1‑71 (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:3714733/1‑71 (MQ=255)
     aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGa                                                                  >  1:2894040/1‑71 (MQ=255)
      acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGAt                                                                 >  1:499976/1‑71 (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:3007464/69‑1 (MQ=11)
                         gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTa                                             <  1:3307294/69‑1 (MQ=11)
                         gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTa                                             <  1:895741/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:2817935/71‑1 (MQ=255)
                          gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt                                            <  1:1913176/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:169435/71‑1 (MQ=11)
                                   ggCGTTTACGCCGCATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTACa                                  <  1:3713134/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)
                                          aCGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGtt                            >  1:27977/1‑71 (MQ=255)
                                          aCGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGtt                            >  1:4200124/1‑71 (MQ=255)
                                           cgccgcATCCGGCAA‑CGGTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTc                                                <  1:4656924/49‑1 (MQ=11)
                                                    cGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGt                   <  1:2852359/70‑1 (MQ=255)
                                                      gCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAg                 <  1:1052830/70‑1 (MQ=255)
                                                           cGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt                                            <  1:2871819/38‑1 (MQ=255)
                                                           cGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:1168418/1‑71 (MQ=255)
                                                           cGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:2656109/1‑71 (MQ=255)
                                                           cGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:423878/1‑71 (MQ=255)
                                                            ggTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:1673876/1‑70 (MQ=255)
                                                            ggTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          <  1:604840/71‑1 (MQ=255)
                                                            ggTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          <  1:2272111/71‑1 (MQ=255)
                                                            ggTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          <  1:1384998/71‑1 (MQ=255)
                                                              tGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          >  1:5650249/1‑69 (MQ=255)
                                                               gCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          >  1:1162796/1‑68 (MQ=255)
                                                               gCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          >  1:4762110/1‑68 (MQ=255)
                                                               gCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt          >  1:2142608/1‑68 (MQ=255)
                                                                cACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:4325146/1‑66 (MQ=255)
                                                                cACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:4695927/1‑66 (MQ=255)
                                                                cACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGa           >  1:5216082/1‑66 (MQ=255)
                                                                cACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATa         >  1:2144495/1‑68 (MQ=255)
                                                                cACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTAAAGCCGTTGCCGAACGTAGGCTGGa           >  1:3206706/1‑66 (MQ=255)
                                                                 aCGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGc     <  1:272331/71‑1 (MQ=255)
                                                                  cGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCg    >  1:2472219/1‑71 (MQ=255)
                                                                   gA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGt   >  1:2886272/1‑71 (MQ=255)
                                                                   gA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGt   <  1:5626522/71‑1 (MQ=255)
                                                                    a‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCg                        <  1:2427619/49‑1 (MQ=255)
                                                                    a‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGtt  <  1:322042/71‑1 (MQ=255)
                                                                    a‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGtt  <  1:1874956/71‑1 (MQ=255)
                                                                     |                                                                       
GCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTT  >  Exported/2773612‑2773750

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