Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,274,991 A→G 31.3% Q101R (CAA→CGA)  yodA → conserved metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,274,9910AG31.3% 103.1 / 60.0 118Q101R (CAA→CGA) yodAconserved metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (34/47);  new base G (19/18);  total (53/65)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 118 total aligned reads displayed.

TGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAAACAA  >  minE/1274934‑1275055
                                                         |                                                                
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:828405/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:3188831/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:2651616/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:3216348/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:255025/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:1945292/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:1985674/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:2487247/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:1200910/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:24220/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:647612/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:2063462/69‑1 (MQ=255)
tGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCga                                                       <  1:823132/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1634266/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:2411671/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1483772/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1451698/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:3094676/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1405973/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1239836/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:439460/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:7594/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:78347/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                                                   <  1:1358915/69‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    <  1:2762984/68‑1 (MQ=255)
    gATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    <  1:1479297/68‑1 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:1524784/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:2965257/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:446633/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:1180766/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:2518486/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:1978889/1‑66 (MQ=255)
      taaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:524003/1‑66 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:2235855/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:237822/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:23422/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:527716/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:2545389/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:3215114/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:2941932/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:1026094/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:3142383/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:311668/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:1028198/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:1244232/63‑1 (MQ=255)
      taaaaccaaaacATTTGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:2010501/63‑1 (MQ=255)
      taaaaccaaaacATATGCTGAAATTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCga                                                       <  1:1414598/63‑1 (MQ=255)
       aaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatgat                                                 >  1:79171/1‑68 (MQ=255)
       aaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatgat                                                 >  1:1253383/1‑68 (MQ=255)
       aaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGATATGAtt                                                >  1:3410071/1‑69 (MQ=255)
       aaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGAtt                                                >  1:110127/1‑69 (MQ=255)
       aaaactaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGAtt                                                >  1:359513/1‑69 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:1049709/68‑1 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:3185841/68‑1 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:1398409/68‑1 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:1548962/68‑1 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:300029/68‑1 (MQ=255)
          actaaaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGc                                             >  1:2916760/1‑69 (MQ=255)
              aaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:2556237/1‑58 (MQ=255)
              aaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:2623321/1‑58 (MQ=255)
              aaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:217681/1‑58 (MQ=255)
              aaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:672958/1‑58 (MQ=255)
              aaacATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgat                                                    >  1:3430050/1‑58 (MQ=255)
                      gCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCgt                                                       <  1:2136443/47‑2 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:1172975/68‑1 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:2276892/68‑1 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:2197176/68‑1 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:2085407/68‑1 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:2818990/68‑1 (MQ=255)
                         gAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAt                               <  1:1621422/68‑1 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:1010053/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:1232529/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:483285/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:463005/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:450169/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:258529/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:672489/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:1908812/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:2726853/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:273841/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:719921/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:750818/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:1291302/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:1912695/1‑69 (MQ=255)
                          aaaTTAAAGATTATTATCGCAAAGGTTATGCGACAGATATCGAAATGATTGGCATTGAGGACGGCAttg                             >  1:2004338/1‑69 (MQ=255)
                             ttAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAttgttg                          <  1:3025248/69‑1 (MQ=255)
                             ttAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCAttgttg                          <  1:1682885/69‑1 (MQ=255)
                               aaaGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGaa                        >  1:1370561/1‑69 (MQ=255)
                                     tattatCACAAAGGTTATGCAACAGATATCGTGATGATTGGCATTg                                         <  1:1211390/46‑1 (MQ=255)
                                        tatCACAAAGGTTATGCAACAGATATCGAGATGATTgg                                              <  1:2918431/38‑1 (MQ=255)
                                            acaAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAaata           >  1:425358/1‑69 (MQ=255)
                                            acaAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAaata           >  1:1560970/1‑69 (MQ=255)
                                            acaAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAaata           >  1:942532/1‑69 (MQ=255)
                                            acaAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAaat            <  1:431889/68‑1 (MQ=255)
                                                    tATGCAACAGATATCGAGATGATTGGCATTGAGGACTGCATTGTTGAATTCCATAGAAATAATGAaac    >  1:274975/1‑68 (MQ=255)
                                                    tATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGaaa     >  1:128384/1‑67 (MQ=255)
                                                    tATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAaac    >  1:575972/1‑68 (MQ=255)
                                                    tATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAaac    >  1:1422792/1‑68 (MQ=255)
                                                    tATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAaac    >  1:1467048/1‑68 (MQ=255)
                                                     aTGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAaacaa  >  1:530460/1‑69 (MQ=255)
                                                         |                                                                
TGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGCATTGTTGAATTCCATAGAAATAATGAAACAA  >  minE/1274934‑1275055

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: