Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,278,430 G→A 39.7% L107L (CTG→TTG)  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,4300GA39.7% 36.9 / 51.2 78L107L (CTG→TTG) btuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (36/11);  new base A (17/14);  total (53/25)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAG  >  minE/2278371‑2278495
                                                           |                                                                 
acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                                          <  1:89359/69‑1 (MQ=255)
acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                                          <  1:1681071/69‑1 (MQ=255)
acGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                                          <  1:709851/69‑1 (MQ=255)
 cGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCg                                                          <  1:1756036/68‑1 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:2450457/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3236452/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3263562/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3351368/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:575522/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:1697370/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:2340097/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:267662/1‑69 (MQ=255)
   cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:1905095/1‑69 (MQ=255)
            cgcAATAGGGAACTGTCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:1126375/1‑68 (MQ=255)
            cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:2746433/1‑68 (MQ=255)
            cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:2484857/1‑68 (MQ=255)
            cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:962531/1‑68 (MQ=255)
            cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:1087193/1‑68 (MQ=255)
            cgcAATAGGGAACTGCCTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCAt                                                  >  1:3236659/1‑65 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:41625/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:1732736/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:2509478/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:1778857/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:2387653/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:1878934/1‑68 (MQ=255)
            cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:1908508/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:2439129/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:2690333/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:1874111/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:3034357/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:765186/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:1196357/1‑63 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAacac                                        >  1:3087180/1‑68 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAaca                                         >  1:1147402/1‑67 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2625955/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:3072583/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:312474/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:660882/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2281550/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2176952/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:145630/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:1427925/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:1899359/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:227042/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2308404/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:1117434/1‑69 (MQ=255)
                   gggAACTGGCTAAGGTCGGCAGAACCACTCACCACCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:727746/1‑69 (MQ=255)
                       aCTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAaca                                   <  1:1625952/69‑1 (MQ=255)
                                 gTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGc                         >  1:2198532/1‑69 (MQ=255)
                                 gTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGc                         >  1:1473256/1‑69 (MQ=255)
                                      cAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTg                    <  1:2657703/69‑1 (MQ=255)
                                                 aCCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAAt          >  1:62894/1‑68 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1344724/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1115028/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1418288/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:739045/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1298129/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2427410/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1738159/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2981969/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1775335/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2577875/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:246808/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:228064/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2424262/69‑1 (MQ=255)
                                                     ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2428266/69‑1 (MQ=255)
                                                        ccAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAg  <  1:1520356/69‑1 (MQ=255)
                                                        ccAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAg  <  1:1710064/69‑1 (MQ=255)
                                                        ccAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAg  <  1:218153/69‑1 (MQ=255)
                                                        ccAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAg  <  1:934990/69‑1 (MQ=255)
                                                           |                                                                 
ACGCTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAG  >  minE/2278371‑2278495

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: