Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,278,440 T→G 41.9% V103V (GTA→GTC 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,4400TG41.9% 18.2 / 76.2 74V103V (GTA→GTCbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (36/7);  new base G (17/14);  total (53/21)
Fisher's exact test for biased strand distribution p-value = 9.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAG  >  minE/2278374‑2278495
                                                                  |                                                       
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:2340097/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:2450457/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:1905095/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:575522/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:267662/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3351368/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3263562/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:1697370/1‑69 (MQ=255)
cTGGACAAGCGCAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGAc                                                       >  1:3236452/1‑69 (MQ=255)
         cgcAATAGGGAACTGTCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:1126375/1‑68 (MQ=255)
         cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:962531/1‑68 (MQ=255)
         cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:1087193/1‑68 (MQ=255)
         cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:2484857/1‑68 (MQ=255)
         cgcAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCaa                                               >  1:2746433/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:1878934/1‑68 (MQ=255)
         cgcAATAGGGAACTGACTAAGATCAGCAGAACCACTCACCCCCGCCAAATTCAGGCGGACGCCATCaa                                               >  1:2387653/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:2439129/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:1196357/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:1941097/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:765186/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:1874111/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:2690333/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAAtt                                             >  1:3034357/1‑63 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAacac                                        >  1:3087180/1‑68 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAaca                                         >  1:1147402/1‑67 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:312474/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2281550/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:660882/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2625955/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:3072583/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:2308404/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:2176952/1‑69 (MQ=255)
                gggAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACAcc                                       >  1:227042/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:1298129/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1344724/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1418288/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:739045/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1115028/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2981969/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:1738159/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:2428266/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2427410/69‑1 (MQ=255)
                                                  ccGCCAAATTCAGGCGGACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATaaaaa     <  1:2424262/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)
                                                                  |                                                       
CTGGACAAGCGCAATAGGGAACTGGCTAAGGTCGGCAGAACCACTCACCCCCGCCAGATTCAGGCGTACGCCATCAATTAACACCAACACATGACTGGCATTTGTACCGCGAATAAAAATAG  >  minE/2278374‑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: