Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,299,396 A→G 6.3% intergenic (‑321/+169) cobU ← / ← yeeX bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,299,3960AG 6.3% 176.3 / 2.9 63intergenic (‑321/+169)cobU/yeeXbifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAAT  >  minE/1299327‑1299465
                                                                     |                                                                     
gCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGa                                                                       >  1:900087/1‑70 (MQ=255)
gCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGa                                                                       >  1:1031785/1‑70 (MQ=255)
   aCTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGa                                                                       >  1:3796610/1‑67 (MQ=255)
    cTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAg                                                                    <  1:1349356/69‑1 (MQ=255)
    cTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAg                                                                    <  1:1373339/69‑1 (MQ=255)
    cTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAg                                                                    <  1:501264/69‑1 (MQ=255)
           aacaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCt                                                           <  1:736260/71‑1 (MQ=255)
            acaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGa                                                                   >  1:884664/1‑62 (MQ=255)
            acaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTg                                                          >  1:1549515/1‑71 (MQ=255)
            acaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTg                                                          >  1:1776204/1‑71 (MQ=255)
              aacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGa                                                                       <  1:795617/56‑1 (MQ=255)
                ccAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATc                                                                 >  1:2742889/1‑60 (MQ=255)
                ccAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATc                                                                 >  1:3360739/1‑60 (MQ=255)
                ccAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTaa                                                       >  1:2418353/1‑70 (MQ=255)
                  aGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTaa                                                       >  1:295951/1‑68 (MQ=255)
                   gCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAg                                                                    <  1:698979/54‑1 (MQ=255)
                    cAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTgcg                                                             >  1:1043765/1‑60 (MQ=255)
                     aaaCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAAc                                                 <  1:2359040/71‑1 (MQ=255)
                      aaCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACt                                                <  1:358104/71‑1 (MQ=255)
                       aCGGGCATTGTGTTCTGAAATCGAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACtt                                               <  1:3273185/71‑1 (MQ=255)
                           gCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTaa                                                       >  1:750387/1‑59 (MQ=255)
                           gCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaa                                           <  1:1654604/71‑1 (MQ=255)
                              ttGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAg                                                                    <  1:844953/43‑1 (MQ=255)
                               tgtgTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGAt                                                                  >  1:764052/1‑44 (MQ=255)
                                gtgtTCTGAAATCCAGGCTATTGATTCAACGTCAGCGa                                                                       >  1:2811974/1‑38 (MQ=255)
                                 tgtTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTg                                                               <  1:581214/45‑1 (MQ=255)
                                   ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACg                                    >  1:753072/1‑70 (MQ=255)
                                   ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGt                                   <  1:3902214/71‑1 (MQ=255)
                                    tCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAAc                                                 >  1:2582390/1‑56 (MQ=255)
                                     cTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa                                 >  1:3636786/1‑71 (MQ=255)
                                     cTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa                                 >  1:1557396/1‑71 (MQ=255)
                                       gAAATCCAGGCTATTGATTAAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTaa                                 >  1:2106251/1‑69 (MQ=255)
                                         aaTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAAt                                >  1:3110752/1‑68 (MQ=255)
                                          aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGaa                            >  1:2986970/1‑71 (MQ=255)
                                           tCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAAt                           <  1:1399187/71‑1 (MQ=255)
                                            ccAGGCTATTGATTCAACGTCAGCGACAGATCTgcgc                                                            >  1:2085122/1‑37 (MQ=255)
                                            ccAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAaac                                        >  1:3733922/1‑57 (MQ=255)
                                             cAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCt                                                           >  1:2128583/1‑37 (MQ=255)
                                             cAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa                                  <  1:1853309/62‑1 (MQ=255)
                                              aGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaa                                           <  1:14849/52‑1 (MQ=255)
                                                gCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATcc                      <  1:2646719/71‑1 (MQ=255)
                                                 cTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGa                             <  1:3770518/63‑1 (MQ=255)
                                                   aTTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCgg                    >  1:2517659/1‑70 (MQ=255)
                                                    ttGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGt                                             <  1:450261/44‑1 (MQ=255)
                                                    ttGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCg                     <  1:3277954/68‑1 (MQ=255)
                                                     tGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt                 >  1:1415213/1‑71 (MQ=255)
                                                      ggttCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAtg                >  1:513489/3‑71 (MQ=255)
                                                      gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtg                >  1:387079/1‑71 (MQ=255)
                                                      gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt                 <  1:1958453/70‑1 (MQ=255)
                                                       aTTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAg                              <  1:3181737/56‑1 (MQ=255)
                                                        ttCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTaaa                                          <  1:156618/43‑1 (MQ=255)
                                                        ttCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt                        <  1:2362783/61‑1 (MQ=255)
                                                         tCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg             >  1:1656376/1‑71 (MQ=255)
                                                          cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa                                  >  1:2355526/1‑49 (MQ=255)
                                                          cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTa                                  >  1:2454941/1‑49 (MQ=255)
                                                          cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt                        >  1:3089553/1‑59 (MQ=255)
                                                          cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg             <  1:299404/70‑1 (MQ=255)
                                                           agcGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc           <  1:1515239/69‑1 (MQ=255)
                                                           aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggt              >  1:2573165/1‑68 (MQ=255)
                                                           aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGc           >  1:1999654/1‑71 (MQ=255)
                                                             cGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc           >  1:2213787/1‑69 (MQ=255)
                                                               tCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAAt                                <  1:546458/46‑1 (MQ=255)
                                                               tCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACt       <  1:176595/71‑1 (MQ=255)
                                                               tCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGCGagc       <  1:2394593/71‑3 (MQ=255)
                                                                  gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtat   <  1:3822283/71‑2 (MQ=255)
                                                                  gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtat   <  1:3785583/71‑2 (MQ=255)
                                                                   cGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaa   >  1:1357433/1‑71 (MQ=255)
                                                                    gACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatgca  >  1:1481928/1‑68 (MQ=255)
                                                                     |                                                                     
GCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAAT  >  minE/1299327‑1299465

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: