Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 776,453 A→T 7.1% intergenic (+825/+40) icd → / ← minE isocitrate dehydrogenase, specific for NADP+/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE776,4530AT 7.1% 162.1 / 4.1 56intergenic (+825/+40)icd/minEisocitrate dehydrogenase, specific for NADP+/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base T (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGG  >  minE/776385‑776521
                                                                    |                                                                    
agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                     >  1:3474811/1‑70 (MQ=255)
    tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                >  1:3260366/1‑71 (MQ=255)
    tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                >  1:1916563/1‑71 (MQ=255)
      aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                              >  1:3172501/1‑71 (MQ=255)
      aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                              >  1:2012210/1‑71 (MQ=255)
        cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                   >  1:3552542/1‑64 (MQ=255)
        cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                   >  1:2507747/1‑64 (MQ=255)
         aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                           <  1:548973/71‑1 (MQ=255)
         aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                           <  1:2832013/71‑1 (MQ=255)
            aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATctct                                                               >  1:1460710/1‑64 (MQ=255)
             aTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                         >  1:2649312/1‑69 (MQ=255)
                 tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                                     <  1:897572/69‑1 (MQ=255)
                   gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                 >  1:3638431/1‑71 (MQ=255)
                    ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                <  1:453763/71‑1 (MQ=255)
                      gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                 >  1:1756953/1‑68 (MQ=255)
                      gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                              <  1:1898275/71‑1 (MQ=255)
                       cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTAACTCTGCCTTGAAGATAAATgc                                              >  1:1001135/1‑70 (MQ=255)
                         ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                       >  1:1751126/1‑59 (MQ=255)
                         ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                       >  1:1396019/1‑59 (MQ=255)
                         ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                       >  1:3593590/1‑59 (MQ=255)
                          tACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                          >  1:624095/1‑71 (MQ=255)
                             cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTg                                                              <  1:2647049/48‑1 (MQ=255)
                              acaACGTGGATGAGAGGGATGCAAAACTCAAGGCAGAGTTATct                                                                 <  1:3398427/44‑1 (MQ=25)
                              acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATct                                                                 <  1:2011866/44‑1 (MQ=38)
                              acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                       >  1:628945/1‑70 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                     <  1:3080871/39‑1 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                       >  1:25068/1‑69 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                       >  1:3098370/1‑69 (MQ=255)
                               caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACa                                     <  1:2438200/71‑1 (MQ=255)
                                aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                    <  1:2248623/55‑1 (MQ=255)
                                aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                    >  1:2950024/1‑71 (MQ=255)
                                aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                    <  1:69960/71‑1 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                               <  1:1651166/70‑1 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                               <  1:1646827/70‑1 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                              <  1:189615/71‑1 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                              >  1:1526558/1‑71 (MQ=255)
                                      gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                              >  1:2871981/1‑71 (MQ=255)
                                       aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                 >  1:1876172/1‑67 (MQ=255)
                                       aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                             <  1:1864710/71‑1 (MQ=255)
                                        tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTa                            >  1:1787073/1‑71 (MQ=255)
                                         gagaGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                           >  1:1618567/1‑71 (MQ=255)
                                          agagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                           <  1:878240/70‑1 (MQ=255)
                                              ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                       >  1:2795716/1‑70 (MQ=255)
                                               gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                  <  1:854733/58‑1 (MQ=255)
                                                aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                  <  1:373709/57‑1 (MQ=255)
                                                 tGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctc                   >  1:2747042/1‑71 (MQ=255)
                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                       >  1:3161510/1‑49 (MQ=255)
                                                   aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtctt                 <  1:2518680/71‑1 (MQ=255)
                                                    aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtctt                 <  1:1927952/70‑1 (MQ=255)
                                                        cTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCt            <  1:2257086/71‑1 (MQ=255)
                                                           aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt           <  1:389740/69‑1 (MQ=255)
                                                               cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTa     <  1:393042/71‑1 (MQ=255)
                                                               cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTa     <  1:1937541/71‑1 (MQ=255)
                                                               cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTa     <  1:3434511/71‑1 (MQ=255)
                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc         >  1:1790938/1‑66 (MQ=255)
                                                                agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt      >  1:2497095/1‑69 (MQ=255)
                                                                  agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg  <  1:2385945/71‑1 (MQ=255)
                                                                  agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg  <  1:20749/71‑1 (MQ=255)
                                                                  agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg  <  1:1394815/71‑1 (MQ=255)
                                                                    aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                     <  1:1148983/50‑1 (MQ=255)
                                                                    aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                     <  1:1262463/50‑1 (MQ=255)
                                                                    |                                                                    
AGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGG  >  minE/776385‑776521

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: