Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,295 C→G 27.2% S10T (AGT→ACT) ‡ 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,2950CG27.2% 219.4 / 95.1 136S10T (AGT→ACT) ‡hydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base C (62/37);  new base G (17/20);  total (79/57)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 138 total aligned reads displayed.

CAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATT  >  minE/1826229‑1826362
                                                                  |                                                                   
caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                                                 <  1:1102946/68‑1 (MQ=255)
caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                                                 <  1:4541114/68‑1 (MQ=255)
 aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgca                                                                >  1:649410/3‑70 (MQ=255)
 aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                                                 <  1:2550428/68‑1 (MQ=255)
  atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                                               >  1:4393609/2‑69 (MQ=255)
  atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                                               >  1:1623890/2‑69 (MQ=255)
  atCCTGAATCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                                               <  1:1165970/70‑3 (MQ=255)
     cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                                                            >  1:4335482/1‑71 (MQ=255)
     cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                                                            >  1:3502090/1‑71 (MQ=255)
     cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg                                                            >  1:1523518/1‑69 (MQ=255)
      tGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                                                                    <  1:3074446/62‑1 (MQ=255)
        aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                                                          <  1:2285167/70‑1 (MQ=255)
        aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg                                                                  >  1:515577/1‑62 (MQ=255)
        aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                                                           <  1:1264433/69‑1 (MQ=255)
        aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc                                                           <  1:4796910/69‑1 (MQ=255)
            tCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:4631497/1‑71 (MQ=255)
              atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCACCCGATACATTTCGTCGCATCCGCaataata                                                    >  1:4067795/2‑69 (MQ=255)
              atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                                                    <  1:5204141/69‑2 (MQ=255)
                gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac                                                 >  1:1356230/1‑70 (MQ=255)
                gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac                                                 >  1:3092677/1‑70 (MQ=255)
                 atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                >  1:860411/1‑71 (MQ=255)
                 atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                >  1:5350741/1‑71 (MQ=255)
                 atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                 <  1:2500988/70‑1 (MQ=255)
                 atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                >  1:4810690/1‑71 (MQ=255)
                  tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg                                               >  1:1061345/1‑71 (MQ=255)
                  tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg                                               >  1:3193442/1‑71 (MQ=255)
                   gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:4546756/1‑64 (MQ=255)
                   gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg                                              >  1:4285982/1‑71 (MQ=255)
                     gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt                                             >  1:3715697/1‑70 (MQ=255)
                     gaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt                                            <  1:4195368/71‑1 (MQ=255)
                      aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtc                                           <  1:1646897/71‑1 (MQ=255)
                      aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                                                            <  1:5528802/54‑1 (MQ=255)
                        actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat                                         >  1:1878049/4‑71 (MQ=255)
                        aCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                                                          >  1:1485272/1‑54 (MQ=255)
                        aCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                                                          >  1:3727298/1‑54 (MQ=255)
                         ctgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatg                                        <  1:363866/69‑2 (MQ=255)
                          tgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAat                                                        <  1:863635/53‑1 (MQ=255)
                           gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                                                    >  1:685906/1‑56 (MQ=255)
                           gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataat                                                     >  1:4456837/1‑56 (MQ=255)
                           gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat                                      >  1:3055161/1‑68 (MQ=255)
                           gcgcAAGCCACTTCACAGGTACGGCAACCGATACAATTCGTCGCATCCGCAATAATAAACCgg                                              >  1:2367061/1‑63 (MQ=255)
                           acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca                                       >  1:2827556/2‑70 (MQ=255)
                           acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca                                       >  1:788160/2‑70 (MQ=255)
                           acgcACGCTACTTCACAGGTACGGCAACCAATACAATTACTCCCGTCAGCa                                                          >  1:1079576/2‑51 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTcgc                                                                 <  1:3769048/43‑1 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTcgc                                                                 <  1:975174/43‑1 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTcgc                                                                 <  1:1759731/43‑1 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGt                                                               >  1:2429850/1‑45 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGt                                                               >  1:4551204/1‑45 (MQ=255)
                            cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                 >  1:2286834/1‑59 (MQ=255)
                              cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                   <  1:3055557/71‑1 (MQ=255)
                               aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc                                  <  1:4349942/71‑1 (MQ=255)
                               aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc                                  <  1:1406533/71‑1 (MQ=255)
                                cGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAtgatga                                                   <  1:220667/53‑1 (MQ=255)
                                 gCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtct                                 >  1:4964328/1‑70 (MQ=255)
                                   tACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                   >  1:2931810/1‑66 (MQ=255)
                                     cTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                                                    <  1:3698722/47‑2 (MQ=255)
                                     cTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                                                    <  1:683848/47‑2 (MQ=255)
                                     cTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                                                    <  1:1849232/47‑2 (MQ=255)
                                     cTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                 >  1:2833369/1‑50 (MQ=255)
                                     cTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGc                            <  1:4611331/71‑1 (MQ=255)
                                      ttCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                <  1:780530/50‑1 (MQ=255)
                                       tCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaa                                                         >  1:1186528/1‑40 (MQ=255)
                                       tCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCaa                          >  1:875849/1‑71 (MQ=255)
                                        cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg                                              >  1:1649699/1‑50 (MQ=255)
                                        cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg                                              >  1:2277003/1‑50 (MQ=255)
                                        cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACTGTTCATCAAATTCTCCAGCAAt                         >  1:2634725/1‑71 (MQ=255)
                                        cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                         >  1:4893344/1‑71 (MQ=255)
                                        cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                         >  1:3200815/1‑71 (MQ=255)
                                         acaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                         <  1:1570595/70‑1 (MQ=255)
                                            ggTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:2382824/1‑39 (MQ=255)
                                            ggTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:4438090/1‑39 (MQ=255)
                                            ggTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:2881310/1‑39 (MQ=255)
                                            ggTACGGCAACCAATACATTTACTCGCGTCAGCAatgat                                                     >  1:4887213/1‑39 (MQ=255)
                                              tACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                   <  1:3427/71‑1 (MQ=255)
                                               aCGGCAACCAATACATTTACTCGCGTCAGCAAttatga                                                   <  1:4235802/38‑1 (MQ=255)
                                               aCGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                   <  1:3380428/70‑1 (MQ=255)
                                                cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                  <  1:3899912/70‑1 (MQ=255)
                                                cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                  <  1:2694846/70‑1 (MQ=255)
                                                cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                  <  1:3680/70‑1 (MQ=255)
                                                 ggCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGAc                      <  1:2180224/65‑1 (MQ=255)
                                                 ggCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGAc                      <  1:2658611/65‑1 (MQ=255)
                                                 ggCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGAc                      <  1:266824/65‑1 (MQ=255)
                                                  gCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                >  1:3059631/1‑70 (MQ=255)
                                                  gCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                >  1:3394854/1‑70 (MQ=255)
                                                     aCCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATgcg            >  1:5187152/1‑71 (MQ=255)
                                                       cAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCa                           >  1:5682124/1‑54 (MQ=255)
                                                       cAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGcc          >  1:3609962/1‑71 (MQ=255)
                                                         aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                   <  1:4153070/60‑1 (MQ=255)
                                                         aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGa        >  1:533537/1‑71 (MQ=255)
                                                          tACATTTACTCGCGTCAGCAATGATGAAACGGTTAATCAAATtc                                  <  1:379588/44‑1 (MQ=255)
                                                           aCATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa      >  1:3074050/1‑71 (MQ=255)
                                                            cATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtct                                 >  1:5587444/1‑43 (MQ=255)
                                                             aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                     <  1:4338875/38‑1 (MQ=255)
                                                             aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                     <  1:3890578/38‑1 (MQ=255)
                                                             aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                     <  1:1935836/38‑1 (MQ=255)
                                                               ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                   >  1:3844548/1‑38 (MQ=255)
                                                               ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACAtt  >  1:1775688/1‑71 (MQ=255)
                                                                  cTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGc                            <  1:3270990/42‑1 (MQ=255)
                                                                  cTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGc                            <  1:980644/42‑1 (MQ=255)
                                                                  |                                                                   
CAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATT  >  minE/1826229‑1826362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: