Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,773,670 Δ1 bp 25.5% intergenic (+59/‑124) 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,6700T.25.5% 102.0 / 45.9 51intergenic (+59/‑124)aceA/aceKisocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (17/21);  new base . (5/8);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAG  >  Exported/2773604‑2773735
                                                                  |                                                                   
aGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCgg                                                                  <  1:2431421/70‑1 (MQ=255)
    cTAAGCAACAACAACAGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCAc                                                             >  1:2642731/1‑71 (MQ=255)
        gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATg                                                        >  1:2750278/1‑71 (MQ=255)
        gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATg                                                        >  1:4513233/1‑71 (MQ=255)
        gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATg                                                        >  1:1488604/1‑71 (MQ=255)
        gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATg                                                        >  1:4473628/1‑71 (MQ=255)
        gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATg                                                        >  1:2082521/1‑71 (MQ=255)
           acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCAc                                                             <  1:651087/64‑1 (MQ=255)
           acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                     <  1:551396/71‑1 (MQ=255)
            caacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                     >  1:5141301/1‑70 (MQ=255)
            caacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTg                                                    <  1:1795518/71‑1 (MQ=255)
             aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGa                                                   >  1:2894040/1‑71 (MQ=255)
             aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGa                                                   >  1:3714733/1‑71 (MQ=255)
              acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGAt                                                  >  1:499976/1‑71 (MQ=255)
                        cTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                                     <  1:2311645/43‑1 (MQ=255)
                          gACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgcg                                     >  1:5337156/1‑71 (MQ=255)
                            cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                   >  1:4337777/1‑71 (MQ=255)
                            cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                   <  1:310363/71‑1 (MQ=255)
                                aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCt                                >  1:4861730/1‑70 (MQ=255)
                                 ggCCGGATAAGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTa                              <  1:2980922/70‑1 (MQ=11)
                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGc                                                <  1:1754404/54‑1 (MQ=255)
                                 ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCtt                               <  1:702953/70‑1 (MQ=255)
                                 gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCtt                               <  1:3613819/68‑1 (MQ=1)
                                 gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTa                              <  1:3007464/69‑1 (MQ=11)
                                 gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTa                              <  1:895741/69‑1 (MQ=11)
                                 gaccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTa                              <  1:3307294/69‑1 (MQ=11)
                                 agCCGGATAAGGCGCTCGCGCCGCATCCGGCAACGGTGCCAGATGCCTGATGCGACGCTTGCGCGTCTTAt                             <  1:1428683/70‑1 (MQ=14)
                                  gCCGGATAAGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTAt                             <  1:4359086/70‑1 (MQ=11)
                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTAt                             <  1:1913176/71‑1 (MQ=255)
                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTAt                             <  1:2817935/71‑1 (MQ=255)
                                  accGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTAt                             <  1:1166876/70‑1 (MQ=11)
                                   ccGGATAAGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCa                           <  1:3793475/71‑1 (MQ=11)
                                   ccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATc                            <  1:4271156/71‑1 (MQ=11)
                                    cGGATAAGGCGCTCGCGCCGCATCCGGCAACGGTGCCAGATGCCTGATGCGACGCTTGCGCGTCTTATCag                          >  1:1015364/1‑70 (MQ=14)
                                     ggATAAGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTAt                             <  1:5187437/67‑1 (MQ=11)
                                     ggATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCag                          >  1:455277/1‑70 (MQ=1)
                                       aTAAGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGcc                       >  1:3814546/1‑71 (MQ=11)
                                         aaGGCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTa                     >  1:234268/1‑71 (MQ=11)
                                          aGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCt                      <  1:3737784/70‑1 (MQ=255)
                                           ggCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTa                     >  1:3284257/1‑69 (MQ=11)
                                           ggCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACa                   <  1:3713134/71‑1 (MQ=11)
                                           ggCGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACa                   <  1:169435/71‑1 (MQ=11)
                                           ggCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTAc                    <  1:423264/71‑1 (MQ=255)
                                             cGTTTACGCCGCATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTAt                             >  1:4935700/1‑59 (MQ=11)
                                               ttCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGc                                                <  1:214201/40‑1 (MQ=40)
                                               ttCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGcc                >  1:4729389/1‑71 (MQ=255)
                                               ttCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGcc                <  1:5093041/71‑1 (MQ=255)
                                                  aCGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGtt             >  1:4200124/1‑71 (MQ=255)
                                                  aCGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGtt             >  1:27977/1‑71 (MQ=255)
                                                   cgccgcATCCGGCAACGGTGCCGACTGCCTGATGCGACGCTTGCGCGTc                                 <  1:4656924/49‑1 (MQ=11)
                                                            cGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGt    <  1:2852359/70‑1 (MQ=255)
                                                              gCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAg  <  1:1052830/70‑1 (MQ=255)
                                                                  |                                                                   
AGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAG  >  Exported/2773604‑2773735

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: