Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,314,660 C→T 24.6% L305L (CTC→CTT hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,314,6600CT24.6% 123.1 / 32.6 97L305L (CTC→CTThisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (47/26);  new base T (11/13);  total (58/39)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGATCATC  >  minE/1314595‑1314715
                                                                 |                                                       
gAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCa                                                      <  1:310515/69‑1 (MQ=255)
   tttACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATac                                                   >  1:701701/1‑69 (MQ=255)
   tttACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATac                                                   >  1:2292637/1‑69 (MQ=255)
   ttaacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATac                                                   >  1:1461699/4‑69 (MQ=255)
     tacctaccAGCGCGTGGGGGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:1295359/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGctct                                                        <  1:346804/62‑1 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2325396/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:1906631/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2100395/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:1411684/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2115493/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2287082/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2390391/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2467964/1‑69 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 <  1:2519250/69‑1 (MQ=255)
     tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATacac                                                 >  1:2602999/1‑69 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:1020409/68‑1 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:2630644/68‑1 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:1861729/68‑1 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:3046647/68‑1 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:16866/68‑1 (MQ=255)
       cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:1613683/68‑1 (MQ=255)
         taccAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCAt                                              >  1:1325888/1‑68 (MQ=255)
         taccAGCGCGTGGGCGATCTCAGCACCGAAAAGAACGAGCACTTCTTCCGATCGCTCTCATACACCATg                                             >  1:1980713/1‑69 (MQ=255)
          accAGCGCTTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:2906739/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:1631995/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:107856/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:23920/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:81620/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:3279815/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:2993358/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:678700/1‑69 (MQ=255)
          accAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATgg                                            >  1:1728738/1‑69 (MQ=255)
             aGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACAcc                                                <  1:2092041/62‑1 (MQ=255)
             aGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                                          >  1:460720/1‑68 (MQ=255)
             aGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                                          >  1:252309/1‑68 (MQ=255)
             aGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                                          >  1:1725017/1‑68 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACg                                    >  1:2926639/1‑68 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:2428061/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:2829735/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:3307337/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:2327383/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:953520/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:718420/1‑69 (MQ=255)
                   tGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                                   >  1:1380635/1‑69 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:1010773/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:3382371/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:2337218/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:70166/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:855511/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:1491344/59‑1 (MQ=255)
                     ggCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:186985/59‑1 (MQ=255)
                      gCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATggg                                           <  1:952533/58‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:2891237/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:3222942/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:1026886/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:917852/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:1326324/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:2289844/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:187208/69‑1 (MQ=255)
                                          aTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:154612/69‑1 (MQ=255)
                                          aTCGAGCACTGCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTaaa            <  1:492390/69‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:600721/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:640459/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:756108/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:3232289/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:3226409/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:958478/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:2239827/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:1092663/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:128007/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:1541444/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:2227285/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:2505487/68‑1 (MQ=255)
                                               gCACTTCTTCCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACg        <  1:2993713/68‑1 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:2836492/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:3111404/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:1300758/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:1332162/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:1437212/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:733557/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:178172/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:2929457/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:2535184/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGat      >  1:3374266/1‑65 (MQ=255)
                                                    tcttcCGTTCGCTTTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGa       >  1:31965/1‑64 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACCCCTGAAAACCAAAGGTAAAAACGatcatc  >  1:609244/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:566816/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:3140858/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:546955/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:2342740/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:1666756/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:3048015/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:150423/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:3292324/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:1105213/1‑69 (MQ=255)
                                                    tcttcCGTTCGCTCTCATACACCAAGGGCGAGACGCTACACCTTAAAACCAAAGGTAAAAACGa       >  1:1254214/1‑64 (MQ=255)
                                                    tcttcAGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGatcatc  >  1:2684736/1‑69 (MQ=255)
                                                                 |                                                       
GAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAGGTAAAAACGATCATC  >  minE/1314595‑1314715

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: