Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,261 C→A 9.4% A21A (GCG→GCT hydN ← formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,826,2610CA 9.4% 233.6 / 4.5 71A21A (GCG→GCThydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base A (5/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTC  >  minE/1826192‑1826330
                                                                     |                                                                     
gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACg                                                                      >  1:2184299/1‑71 (MQ=255)
  tACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCt                                                                    <  1:3196766/71‑1 (MQ=255)
      cgGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACtt                                                                <  1:843693/71‑1 (MQ=255)
            aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAgg                                                          <  1:340263/71‑1 (MQ=255)
            aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                          >  1:2277659/1‑71 (MQ=255)
            aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                          >  1:158531/1‑71 (MQ=255)
             aaaGTTTCCGGGGTCAGCGATGAACAGTCCTGATTTTCCTGATGAGAAACAGCGCACGCTACTTCACAgg                                                          >  1:2851255/1‑70 (MQ=255)
               aGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTAc                                                       >  1:3113173/1‑71 (MQ=255)
               aGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                                           >  1:4043734/1‑67 (MQ=255)
                 tttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACg                                                      <  1:2726939/70‑1 (MQ=255)
                 tttCCGGGGTCAGCGATGCACAGTCCTGATTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACg                                                      <  1:110761/69‑1 (MQ=255)
                   tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGc                                                    <  1:1558496/70‑1 (MQ=255)
                   tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGc                                                    <  1:1274855/70‑1 (MQ=255)
                      ggggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAAcc                                                <  1:2179737/71‑1 (MQ=255)
                          tCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCAc                                                                       <  1:257502/44‑1 (MQ=255)
                          tCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAAc                                                 <  1:2993132/66‑1 (MQ=255)
                                aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCaa                                                  >  1:1007957/1‑59 (MQ=255)
                                aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAAcc                                                <  1:3346515/61‑1 (MQ=255)
                                aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTAATTCACAGGTACGGCAACCAATACATTTa                                      >  1:3282187/1‑71 (MQ=255)
                                  gCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTcaca                                                            >  1:2153449/1‑47 (MQ=255)
                                  gCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                                    >  1:591218/1‑71 (MQ=255)
                                   cacaGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCaa                                              >  1:2279165/1‑60 (MQ=255)
                                   cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTc                                   >  1:1911597/1‑71 (MQ=255)
                                    acaGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTa                                                                   >  1:92698/1‑38 (MQ=255)
                                    acaGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg                                  >  1:1977199/1‑71 (MQ=255)
                                     caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACa                                          >  1:2551913/4‑62 (MQ=255)
                                      aGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATAc                                           <  1:3663530/60‑1 (MQ=255)
                                        tCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACgg                                                     >  1:1371467/1‑48 (MQ=255)
                                        tCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGt                               >  1:2465707/1‑70 (MQ=255)
                                         ccTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATc                              >  1:2025834/1‑70 (MQ=255)
                                          cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCa                                               >  1:1776352/1‑52 (MQ=255)
                                          cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                            >  1:3722088/1‑71 (MQ=255)
                                          cTGATTTTCCTGATGAGAAACCACTCACGCTACTTCACAGGTACgg                                                     >  1:196872/1‑46 (MQ=39)
                                          cTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                                    >  1:3217621/1‑47 (MQ=255)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg                            >  1:418110/1‑69 (MQ=255)
                                           tGATTTTCCTGATGAGTAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCa                             <  1:3630492/69‑1 (MQ=255)
                                           tGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCaa                                              >  1:2706893/1‑52 (MQ=255)
                                           tGATTTTCCTGATGAGAAACCGCGAACGCTACTTCACAGGTACGGCAACCAATACAttt                                       <  1:178181/59‑1 (MQ=255)
                                            gATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACAtt                                        >  1:2293343/1‑57 (MQ=255)
                                            gATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTa                                      <  1:2537935/59‑1 (MQ=255)
                                            gATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                           >  1:2662702/1‑70 (MQ=255)
                                            gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTa                                      <  1:2375934/59‑1 (MQ=255)
                                             aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                           <  1:4064495/69‑1 (MQ=255)
                                             aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                           <  1:3131823/69‑1 (MQ=255)
                                              ttttCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAat                        <  1:1333237/71‑1 (MQ=255)
                                               tttCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAat                        >  1:2002547/1‑70 (MQ=255)
                                               tttCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAat                        >  1:246115/1‑70 (MQ=255)
                                                ttCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatg                       <  1:1890953/70‑1 (MQ=255)
                                                   cTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAtgatga                   <  1:3990270/71‑1 (MQ=255)
                                                    tgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATa                                            <  1:585440/45‑1 (MQ=255)
                                                     gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATa                                            >  1:3579189/1‑44 (MQ=16)
                                                         agaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                 >  1:1496109/1‑67 (MQ=255)
                                                         agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt             >  1:2057456/1‑71 (MQ=255)
                                                         agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt             >  1:251758/1‑71 (MQ=255)
                                                          gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt            >  1:811174/1‑71 (MQ=255)
                                                           aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                     <  1:1169741/61‑1 (MQ=255)
                                                           aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt            >  1:2894175/1‑70 (MQ=255)
                                                           aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt            >  1:1822717/1‑70 (MQ=255)
                                                           aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtc           <  1:3534467/71‑1 (MQ=255)
                                                             aCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcat         >  1:3215940/1‑71 (MQ=255)
                                                              ctgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCg               <  1:557322/62‑1 (MQ=255)
                                                              ccGCGCACGCTACTTCACAGGTACGGCAACCAATACAtt                                        <  1:4057081/39‑1 (MQ=255)
                                                               cgcgCACGCTACTTCACAGGTACGGCAACCAATACAtt                                        <  1:3038027/38‑1 (MQ=255)
                                                               cacgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcat         >  1:1757785/3‑69 (MQ=255)
                                                                gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat      >  1:2577913/1‑68 (MQ=255)
                                                                 cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                            <  1:1236999/48‑1 (MQ=255)
                                                                 cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa     <  1:3502761/71‑1 (MQ=255)
                                                                 cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa     <  1:1422001/71‑1 (MQ=255)
                                                                 cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa     <  1:801564/71‑1 (MQ=255)
                                                                   cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAt    >  1:1935472/1‑70 (MQ=255)
                                                                    aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc  <  1:3257930/71‑1 (MQ=255)
                                                                    aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc  <  1:2618879/71‑1 (MQ=255)
                                                                     cGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt            >  1:1075433/1‑60 (MQ=255)
                                                                     |                                                                     
GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTC  >  minE/1826192‑1826330

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: