Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,278,404 G→A 33.7% D115D (GAC→GAT btuB ← vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,278,4040GA33.7% 56.0 / 35.9 71D115D (GAC→GATbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (33/14);  new base A (18/6);  total (51/20)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGC  >  minE/2278340‑2278472
                                                                |                                                                    
cAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCgg                                                                  <  1:1899446/69‑1 (MQ=255)
cAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCgg                                                                  <  1:1413371/69‑1 (MQ=255)
cAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCgg                                                                  <  1:3427642/69‑1 (MQ=255)
   cGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAg                                                               <  1:749399/69‑1 (MQ=255)
   cGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAg                                                               <  1:1596690/69‑1 (MQ=255)
        cgcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:2923044/69‑1 (MQ=255)
        cgcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:1904169/69‑1 (MQ=255)
        cgcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:1161841/69‑1 (MQ=255)
        cgcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:1268772/69‑1 (MQ=255)
        cgcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:561699/69‑1 (MQ=255)
         gcgGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCa                                                          <  1:1611373/68‑1 (MQ=255)
                     tatatTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTATGGTCgg                                                                  <  1:557064/48‑1 (MQ=255)
                             acacGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCa                                                          <  1:2612304/48‑1 (MQ=255)
                             acacGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCa                                                          <  1:1129110/48‑1 (MQ=255)
                               acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                   <  1:1681071/69‑1 (MQ=255)
                               acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                   <  1:709851/69‑1 (MQ=255)
                               acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                   <  1:89359/69‑1 (MQ=255)
                                cGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                   <  1:1756036/68‑1 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:575522/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:3351368/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:2340097/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:3263562/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:1697370/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:2450457/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:267662/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:3236452/1‑69 (MQ=255)
                                  cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                >  1:1905095/1‑69 (MQ=255)
                                           cgcAATAGGGAACTGTCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                        >  1:1126375/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                        >  1:962531/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                        >  1:2746433/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                        >  1:1087193/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                        >  1:2484857/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGCCTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCAt                           >  1:3236659/1‑65 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACt                                                        >  1:196886/1‑36 (MQ=25)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:2509478/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:2387653/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:1908508/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:1878934/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:1778857/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:1732736/1‑68 (MQ=255)
                                           cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                        >  1:41625/1‑68 (MQ=255)
                                              aaTAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTaa                    <  1:2267532/69‑1 (MQ=255)
                                                  gtgAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1396892/3‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGc                              >  1:2449884/1‑55 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaaat                      >  1:1933855/1‑61 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:1941097/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:1874111/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:3034357/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:2690333/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:1196357/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:765186/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                      >  1:2439129/1‑63 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAacac                 >  1:3087180/1‑68 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAaca                  >  1:1147402/1‑67 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:660882/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1117434/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1427925/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:145630/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1475264/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1655858/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1692203/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:1899359/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:312474/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:3072583/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:2176952/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:2625955/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:227042/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:2281550/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:2308404/1‑69 (MQ=255)
                                                  gggAACTGGCTAAGGTCGGCAGAACCACTCACCACCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                >  1:727746/1‑69 (MQ=255)
                                                      aCTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAaca            <  1:1625952/69‑1 (MQ=255)
                                                                gTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGc  >  1:1473256/1‑69 (MQ=255)
                                                                gTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGc  >  1:2198532/1‑69 (MQ=255)
                                                                |                                                                    
CAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGC  >  minE/2278340‑2278472

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: