Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,315 A→C 7.8% R3R (CGT→CGG 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,3150AC 7.8% 261.2 / 14.8 103R3R (CGT→CGGhydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base A (46/49);  new base C (4/4);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 111 total aligned reads displayed.

ATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTT  >  minE/1826246‑1826382
                                                                     |                                                                   
atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                                    >  1:5350741/1‑71 (MQ=255)
atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                                    >  1:860411/1‑71 (MQ=255)
atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                                     <  1:2500988/70‑1 (MQ=255)
atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                                    >  1:2564220/1‑71 (MQ=255)
atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                                    >  1:4810690/1‑71 (MQ=255)
 tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg                                                                   >  1:1061345/1‑71 (MQ=255)
 tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg                                                                   >  1:3193442/1‑71 (MQ=255)
  gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg                                                                  >  1:4285982/1‑71 (MQ=255)
   agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt                                                                 >  1:4197328/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)
       actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat                                                             <  1:846580/68‑1 (MQ=255)
       actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat                                                             >  1:1878049/4‑71 (MQ=255)
        ctgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatg                                                            <  1:363866/69‑2 (MQ=255)
          gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat                                                          >  1:3596661/1‑68 (MQ=255)
          gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat                                                          >  1:3055161/1‑68 (MQ=255)
          gcgcAAGCCACTTCACAGGTACGGCAACCGATACAATTCGTCGCATCCGCAATAATAAACCgg                                                                  >  1:2367061/1‑63 (MQ=255)
          acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca                                                           >  1:4714803/2‑70 (MQ=255)
          acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca                                                           >  1:788160/2‑70 (MQ=255)
          acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca                                                           >  1:2827556/2‑70 (MQ=255)
           cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                                     >  1:2286834/1‑59 (MQ=255)
             cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                                       <  1:4459305/71‑1 (MQ=255)
             cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                                       <  1:3055557/71‑1 (MQ=255)
              aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc                                                      <  1:4349942/71‑1 (MQ=255)
              aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc                                                      <  1:1406533/71‑1 (MQ=255)
                  tACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                                       >  1:2931810/1‑66 (MQ=255)
                    cTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                                                                     >  1:2833369/1‑50 (MQ=255)
                    cTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGc                                                <  1:4611331/71‑1 (MQ=255)
                     ttCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                                                                    <  1:780530/50‑1 (MQ=255)
                      tCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCaa                                              >  1:875849/1‑71 (MQ=255)
                       cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg                                                                  >  1:2277003/1‑50 (MQ=255)
                       cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACTGTTCATCAAATTCTCCAGCAAt                                             >  1:2634725/1‑71 (MQ=255)
                       cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCaa                                              >  1:761459/1‑70 (MQ=255)
                       cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                                             >  1:4893344/1‑71 (MQ=255)
                       cacaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                                             >  1:3200815/1‑71 (MQ=255)
                        acaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCGAGCAAt                                             <  1:4551119/70‑1 (MQ=255)
                        acaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAt                                             <  1:1570595/70‑1 (MQ=255)
                             tACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                                       <  1:1846860/71‑1 (MQ=255)
                             tACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                                       <  1:3427/71‑1 (MQ=255)
                              aCGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                                       <  1:3380428/70‑1 (MQ=255)
                               cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                                      <  1:3680/70‑1 (MQ=255)
                               cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                                      <  1:3899912/70‑1 (MQ=255)
                               cGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGtt                                      <  1:2694846/70‑1 (MQ=255)
                                ggCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGAc                                          <  1:266824/65‑1 (MQ=255)
                                ggCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGAc                                          <  1:2658611/65‑1 (MQ=255)
                                 gCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    >  1:3059631/1‑70 (MQ=255)
                                 gCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    >  1:1607714/1‑70 (MQ=255)
                                 gCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    >  1:3394854/1‑70 (MQ=255)
                                    aCCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATgcg                                >  1:5187152/1‑71 (MQ=255)
                                      cAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGcc                              >  1:3609962/1‑71 (MQ=255)
                                        aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGt                                       <  1:4153070/60‑1 (MQ=255)
                                        aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGa                            >  1:533537/1‑71 (MQ=255)
                                        aTACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGa                            >  1:5391877/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)
                                           cATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTcc                                                   >  1:2136589/1‑45 (MQ=255)
                                            aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                                         <  1:3890578/38‑1 (MQ=255)
                                            aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                                         <  1:4338875/38‑1 (MQ=255)
                                            aTTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa                                                         <  1:1935836/38‑1 (MQ=255)
                                              ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt                                                       >  1:3844548/1‑38 (MQ=255)
                                              ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAAt                                                        >  1:4324407/1‑37 (MQ=255)
                                              ttACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACAtt                      >  1:1775688/1‑71 (MQ=255)
                                                 cTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGc                                                <  1:980644/42‑1 (MQ=255)
                                                  tCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    <  1:1984563/53‑1 (MQ=255)
                                                  tCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    <  1:228687/53‑1 (MQ=255)
                                                   cgcgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    <  1:3071083/52‑1 (MQ=255)
                                                    gcgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTcaca                  <  1:3063555/69‑1 (MQ=255)
                                                    gcgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTcaca                  <  1:4505748/69‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCa                                               <  1:1916927/37‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    <  1:5405472/48‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGcc                              <  1:4544233/54‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGcc                              <  1:1102610/54‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa                          <  1:1397349/58‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa                          <  1:5242618/58‑1 (MQ=255)
                                                       tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa                          <  1:3686111/58‑1 (MQ=255)
                                                             atgatgAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAg       <  1:2140160/71‑1 (MQ=255)
                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTa                                     <  1:2309971/40‑1 (MQ=255)
                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTa                                     <  1:5614095/40‑1 (MQ=255)
                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa                          <  1:1178277/51‑1 (MQ=255)
                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATa                          <  1:3464499/51‑1 (MQ=255)
                                                              tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGt      <  1:2305325/71‑1 (MQ=255)
                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTaa                                    >  1:4153284/1‑39 (MQ=255)
                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATTCGCCGATACATTCACAAATCa             >  1:130738/1‑62 (MQ=255)
                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATc              >  1:3472620/1‑61 (MQ=255)
                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGttt    >  1:3257437/1‑71 (MQ=255)
                                                                atgaAACGGTTCATCAAATTCTCCAGCAATGAAAGt                                       >  1:4018804/1‑36 (MQ=255)
                                                                 tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGtttt   <  1:875600/71‑1 (MQ=255)
                                                                 tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGtttt   >  1:439443/1‑71 (MQ=255)
                                                                 tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGtttt   >  1:1793718/1‑71 (MQ=255)
                                                                  gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGcc         >  1:96935/1‑64 (MQ=255)
                                                                  gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGcc         >  1:4317757/1‑64 (MQ=255)
                                                                  gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGcc         >  1:2679566/1‑64 (MQ=255)
                                                                  gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGcc         >  1:2916635/1‑64 (MQ=255)
                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACa                        >  1:3546946/1‑48 (MQ=255)
                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACa                        >  1:3843278/1‑48 (MQ=255)
                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACa                        >  1:1137824/1‑48 (MQ=255)
                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGttttt  >  1:2400256/1‑70 (MQ=255)
                                                                   aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGttttt  >  1:5298949/1‑70 (MQ=255)
                                                                     |                                                                   
ATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTT  >  minE/1826246‑1826382

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