Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,256 G→A 32.3% A23V (GCG→GTG)  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,2560GA32.3% 76.9 / 53.0 62A23V (GCG→GTG) hydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base G (26/16);  new base A (10/10);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAA  >  minE/1826187‑1826326
                                                                     |                                                                      
aCATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcgc                                                                       >  1:3304531/1‑71 (MQ=255)
aCATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcgc                                                                       >  1:2569633/1‑71 (MQ=255)
  aTGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCgcgc                                                                     <  1:437912/71‑1 (MQ=255)
     gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACg                                                                  >  1:2184299/1‑71 (MQ=255)
       tACGCGGTAAAAAAGTTTCCGGGGTCAGCTATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCa                                                                    <  1:3552862/67‑1 (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:158531/1‑71 (MQ=255)
                 aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                      >  1:2277659/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:1274855/70‑1 (MQ=255)
                        tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGc                                                <  1:1558496/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:3131823/69‑1 (MQ=255)
                                                  aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                       <  1:4064495/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:251758/1‑71 (MQ=255)
                                                              agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt         >  1:2057456/1‑71 (MQ=255)
                                                               gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt        >  1:811174/1‑71 (MQ=255)
                                                                aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                 <  1:1169741/61‑1 (MQ=255)
                                                                aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt        >  1:1822717/1‑70 (MQ=255)
                                                                aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt        >  1:2894175/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)
                                                                     gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat  >  1:2577913/1‑68 (MQ=255)
                                                                     |                                                                      
ACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAA  >  minE/1826187‑1826326

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: