Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,505,809 C→A 12.7% G510G (GGG→GGT yfbS ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,505,8090CA12.7% 136.0 / 6.5 62G510G (GGG→GGTyfbSpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base A (5/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGC  >  minE/1505747‑1505871
                                                              |                                                              
aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:985875/1‑69 (MQ=255)
aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:2361729/1‑69 (MQ=255)
aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:3166209/1‑69 (MQ=255)
aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:403168/1‑69 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:896376/1‑69 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:1496723/1‑69 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:592760/69‑1 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:3315647/1‑69 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:2889013/69‑1 (MQ=255)
    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:2610844/69‑1 (MQ=255)
      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:858879/69‑1 (MQ=255)
      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1407476/69‑1 (MQ=255)
      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1722437/69‑1 (MQ=255)
      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGAAAAGACAAACAAACACCCCAGCATCata                                                    <  1:1507605/69‑1 (MQ=255)
       ccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:941142/68‑1 (MQ=255)
                tATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGt                                          <  1:1416547/69‑1 (MQ=255)
                tATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGt                                          <  1:2435173/69‑1 (MQ=255)
                 aTTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:153997/69‑1 (MQ=255)
                 aTTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1891557/69‑1 (MQ=255)
                  ttAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1102302/68‑1 (MQ=255)
                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGc                                              >  1:1085710/1‑61 (MQ=255)
                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAacc                                      >  1:2598292/1‑69 (MQ=255)
                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAacc                                      >  1:1505080/1‑69 (MQ=255)
                       gaTAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1015040/52‑1 (MQ=255)
                          aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:1738058/69‑1 (MQ=255)
                          aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:2917957/69‑1 (MQ=255)
                           aaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:1571196/68‑1 (MQ=255)
                             tAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAAt                             >  1:2951204/1‑69 (MQ=255)
                              aGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATg                            <  1:1335371/69‑1 (MQ=255)
                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1393920/50‑1 (MQ=255)
                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:3197563/50‑1 (MQ=255)
                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:2720259/50‑1 (MQ=255)
                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1131323/50‑1 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:3273754/37‑1 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1908795/37‑1 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCa                        >  1:2069350/1‑65 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:648902/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:504616/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:421036/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:713820/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1030110/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:323409/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:781215/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2050538/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2642875/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:196487/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2240270/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1955936/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2521224/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2672126/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1815605/1‑69 (MQ=255)
                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATAAg                    >  1:880651/1‑69 (MQ=255)
                                              acaagacaaACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTc             <  1:2837577/68‑1 (MQ=255)
                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:1551177/69‑1 (MQ=255)
                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:1829722/69‑1 (MQ=255)
                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:2158003/69‑1 (MQ=255)
                                                   acaaacaaacaaCCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:468041/68‑1 (MQ=255)
                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2798308/1‑69 (MQ=255)
                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:695049/1‑69 (MQ=255)
                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2944992/1‑69 (MQ=255)
                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2139193/1‑69 (MQ=255)
                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2982651/1‑69 (MQ=255)
                                                        caaacaCCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAgcgc  <  1:3138759/69‑1 (MQ=255)
                                                              |                                                              
AACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGC  >  minE/1505747‑1505871

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