Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,264 T→C 37.3% V20V (GTA→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,2640TC37.3% 141.4 / 149.5 150V20V (GTA→GTGhydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (56/38);  new base C (31/25);  total (87/63)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 151 total aligned reads displayed.

TAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCT  >  minE/1826201‑1826331
                                                               |                                                                   
tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                                              <  1:612773/71‑1 (MQ=255)
     aaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTa                                                         >  1:2947631/1‑71 (MQ=255)
         ttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACgg                                                      >  1:737801/1‑70 (MQ=255)
          tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                                              >  1:1112419/1‑61 (MQ=255)
             ggggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAAc                                                  >  1:5412158/1‑70 (MQ=255)
             ggggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGTCAAcc                                                 <  1:3590309/71‑1 (MQ=255)
               ggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCaa                                               <  1:2852499/71‑1 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACTCACGCTACTACACAGGTACGGCAACCAATACa                                           >  1:508002/1‑71 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                                         >  1:3961321/1‑57 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                                         >  1:3008819/1‑57 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                                               >  1:4262611/1‑67 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                                               >  1:4462017/1‑67 (MQ=255)
                   aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa                                           >  1:1007127/1‑71 (MQ=255)
                     cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                                         <  1:4151053/71‑1 (MQ=255)
                     cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                                         >  1:3295559/1‑71 (MQ=255)
                     cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                                         <  1:1007139/71‑1 (MQ=255)
                       aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTa                                       >  1:1894282/1‑71 (MQ=255)
                       aTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                           <  1:4132488/51‑1 (MQ=255)
                       aTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTa                                       <  1:1973093/71‑1 (MQ=255)
                         gCACAGTCCTGATTTTGCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                                         >  1:3382647/1‑51 (MQ=38)
                         gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                                            >  1:1436344/1‑48 (MQ=255)
                         gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                                            >  1:4692973/1‑48 (MQ=255)
                          cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTAc                                                        >  1:2238228/1‑51 (MQ=255)
                            caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                  <  1:4541114/68‑1 (MQ=255)
                            caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                  <  1:3192130/68‑1 (MQ=255)
                             aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAgg                                                           <  1:1078367/43‑1 (MQ=12)
                             aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgca                                 >  1:649410/3‑70 (MQ=255)
                             aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                                  <  1:2550428/68‑1 (MQ=255)
                              gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                           <  1:1449875/44‑1 (MQ=255)
                              gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                                     >  1:4059290/1‑66 (MQ=255)
                              atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                >  1:1623890/2‑69 (MQ=255)
                              atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                >  1:4393609/2‑69 (MQ=255)
                              atCCTGAATCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                                <  1:1165970/70‑3 (MQ=255)
                               tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                                     <  1:1982038/49‑1 (MQ=255)
                               tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                                     <  1:5149652/49‑1 (MQ=255)
                               tCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACg                                                       >  1:2740211/1‑47 (MQ=16)
                                 cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                             >  1:880577/1‑71 (MQ=255)
                                 cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                             >  1:4335482/1‑71 (MQ=255)
                                 cTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAt                                          >  1:1394768/1‑58 (MQ=255)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:4574773/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:386318/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:168476/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:5273999/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:3370833/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                                         >  1:5233976/1‑43 (MQ=12)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAAccga                                               <  1:1103046/53‑3 (MQ=255)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAAccga                                               <  1:2051226/53‑3 (MQ=255)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                                          >  1:4947987/1‑58 (MQ=255)
                                 cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg                             >  1:1523518/1‑69 (MQ=255)
                                   gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                                     >  1:3681554/1‑61 (MQ=255)
                                   gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                                     >  1:3799546/1‑61 (MQ=255)
                                   gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                            >  1:5298515/1‑70 (MQ=255)
                                   gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                                                    <  1:4620319/46‑1 (MQ=16)
                                   gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                                                    <  1:2131427/46‑1 (MQ=16)
                                   gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc                            >  1:4850430/1‑70 (MQ=255)
                                    aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                           >  1:3737640/1‑70 (MQ=255)
                                    aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                           <  1:2285167/70‑1 (MQ=255)
                                    aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                            <  1:1264433/69‑1 (MQ=255)
                                    aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc                            <  1:3267245/69‑1 (MQ=255)
                                    aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCa                           <  1:3064599/70‑1 (MQ=255)
                                       ttCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                                     <  1:2303942/41‑1 (MQ=255)
                                          atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                     <  1:5204141/69‑2 (MQ=255)
                                            gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac                  >  1:3092677/1‑70 (MQ=255)
                                            gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac                  >  1:1356230/1‑70 (MQ=255)
                                             atgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAttt                                        >  1:1321704/1‑48 (MQ=19)
                                             atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaaa                  <  1:2500988/70‑1 (MQ=255)
                                             atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                 >  1:4810690/1‑71 (MQ=255)
                                             atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc                 >  1:2564220/1‑71 (MQ=255)
                                               gagaAACCGCGCACGCTACTTCACAGGTACGGCaa                                                   >  1:1746774/1‑35 (MQ=255)
                                                agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt              >  1:4197328/1‑71 (MQ=255)
                                                 gaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaa                   >  1:113282/1‑65 (MQ=255)
                                                  aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtc            <  1:1646897/71‑1 (MQ=255)
                                                  aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                             <  1:5528802/54‑1 (MQ=255)
                                                    actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat          >  1:1878049/4‑71 (MQ=255)
                                                    actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat          <  1:846580/68‑1 (MQ=255)
                                                    aCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                           >  1:3727298/1‑54 (MQ=255)
                                                     ctgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatg         <  1:437939/69‑2 (MQ=255)
                                                     ctgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatg         <  1:363866/69‑2 (MQ=255)
                                                      tgcgcAAGCCACTTCACAGGTACGGCAACCGATACAtt                                         <  1:3410947/37‑1 (MQ=25)
                                                      tgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAat                         <  1:863635/53‑1 (MQ=255)
                                                       gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata                     >  1:685906/1‑56 (MQ=255)
                                                       gcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatgat       >  1:3596661/1‑68 (MQ=255)
                                                       acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca        >  1:4714803/2‑70 (MQ=255)
                                                       acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca        >  1:788160/2‑70 (MQ=255)
                                                       acgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcatca        >  1:2827556/2‑70 (MQ=255)
                                                       acgcACGCTACTTCACAGGTACGGCAACCAATACAATTACTCCCGTCAGCa                           >  1:1079576/2‑51 (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:2350963/1‑45 (MQ=255)
                                                        cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGt                                >  1:4551204/1‑45 (MQ=255)
                                                        cgcAAGCCACTTCACAGGTACGGCAACCGATACAtttc                                       >  1:4546304/1‑37 (MQ=25)
                                                        cgcAAGCCACTTCACAGGTACGGCAACCGATACAtttc                                       >  1:46261/1‑37 (MQ=25)
                                                         gcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgatg                     <  1:227899/55‑1 (MQ=255)
                                                         gcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataataa                    <  1:2872217/56‑3 (MQ=255)
                                                          cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt    <  1:3055557/71‑1 (MQ=255)
                                                          cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt    <  1:4459305/71‑1 (MQ=255)
                                                           aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc   <  1:4349942/71‑1 (MQ=255)
                                                            cGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAtgatga                    <  1:220667/53‑1 (MQ=255)
                                                             gCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtct  >  1:4964328/1‑70 (MQ=255)
                                                               |                                                                   
TAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCT  >  minE/1826201‑1826331

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: