Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,286,261 A→T 22.0% pseudogene (379/1053 nt)
pseudogene (379/705 nt)
yeeL ←
yeeL ←
ECK1975:JW1961+JW5325:b4497; hypothetical protein
ECK1975:JW1961:b1980; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,286,2610AT22.0% 98.5 / 23.9 59pseudogene (379/1053 nt)
pseudogene (379/705 nt)
yeeL
yeeL
ECK1975:JW1961+JW5325:b4497; hypothetical protein
ECK1975:JW1961:b1980; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base A (16/30);  new base T (2/11);  total (18/41)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

AATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTA  >  minE/1286201‑1286318
                                                            |                                                         
aaTTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGttt                                                   <  1:123934/69‑1 (MQ=255)
aaTTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGttt                                                   <  1:3136665/69‑1 (MQ=255)
aaTTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGttt                                                   <  1:2701481/69‑1 (MQ=255)
aaTTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGttt                                                   <  1:1875762/69‑1 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:692183/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:3210882/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:1341991/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:1470468/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:2983526/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:2194589/1‑68 (MQ=255)
           ataCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCtt                                         >  1:2079086/1‑68 (MQ=255)
              ccAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTa                                     >  1:925975/1‑69 (MQ=255)
              ccAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTa                                     >  1:1110041/1‑69 (MQ=255)
              ccAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTa                                     >  1:687139/1‑69 (MQ=255)
              ccAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTa                                     >  1:3071495/1‑69 (MQ=255)
              ccAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTa                                     >  1:2765897/1‑69 (MQ=255)
                aacaacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACg                                   <  1:1469781/69‑1 (MQ=255)
                aacaacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACg                                   <  1:1125994/69‑1 (MQ=255)
                aacaacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACg                                   <  1:2312764/69‑1 (MQ=255)
                aacaacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACg                                   <  1:2273072/69‑1 (MQ=255)
                aacaacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACg                                   <  1:2125603/69‑1 (MQ=255)
                   aacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCAt                                >  1:76410/1‑69 (MQ=255)
                   aacCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCAt                                >  1:2964739/1‑69 (MQ=255)
                    acCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:226118/69‑1 (MQ=255)
                    acCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:2296919/69‑1 (MQ=255)
                    acCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:1626712/69‑1 (MQ=255)
                    acCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:197511/69‑1 (MQ=255)
                    acCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:1202699/69‑1 (MQ=255)
                     ccTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:32376/68‑1 (MQ=255)
                        ggaaaTTCGTTTTTGCTTTCAACTGACTGATAAATATCATCGTTTTCAGGTTCTTTTTACGCATcgcc                            >  1:554190/3‑68 (MQ=255)
                        ggaaaTTCGTTTTTGCTTTCAACTGACTGATAAATATCATCGTTTTCAGGTTCTTTTTACGCATcgcc                            >  1:3170287/3‑68 (MQ=255)
                           aaTTCGTTTTTGCTTTCAATTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:69612/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:2930279/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:326362/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:3198658/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:3084311/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:161872/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:2720064/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:2679243/67‑1 (MQ=255)
                           aaTTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:2105243/67‑1 (MQ=255)
                             ttCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAAt                       >  1:691152/1‑68 (MQ=255)
                             ttCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAAt                       >  1:2073533/1‑68 (MQ=255)
                             ttCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATg                      <  1:2349564/69‑1 (MQ=255)
                                     tGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGAt               <  1:1389554/68‑1 (MQ=255)
                                     tGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGAt               <  1:1192714/68‑1 (MQ=255)
                                      gCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATc                               <  1:77135/51‑1 (MQ=255)
                                        tttCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                          <  1:648643/54‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:32214/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:339392/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:3426393/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:1335133/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:562554/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:2969632/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:2964526/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:1805576/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:2641996/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:1879271/67‑1 (MQ=255)
                                                   tgaTAAATATCATCGTTTTCAGGTTCTTTTTACGCATCGCCTCTATACAAAGATCCGGCGTACCGTa  <  1:1837672/67‑1 (MQ=255)
                                                    gaTAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGa                <  1:3348770/52‑1 (MQ=255)
                                                            |                                                         
AATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTA  >  minE/1286201‑1286318

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