Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,255 C→T 32.5% A23A (GCG→GCA) ‡ 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,2550CT32.5% 161.8 / 94.0 117A23A (GCG→GCA) ‡hydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base C (49/30);  new base T (21/17);  total (70/47)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 122 total aligned reads displayed.

GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTC  >  minE/1826192‑1826321
                                                               |                                                                  
gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcacgcacg                                                             >  1:1488008/1‑71 (MQ=255)
gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcacgcacg                                                             >  1:3602715/1‑71 (MQ=255)
gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcacgcacg                                                             >  1:4261569/1‑71 (MQ=255)
         tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                                    <  1:612773/71‑1 (MQ=255)
            aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                 <  1:4945897/71‑1 (MQ=255)
             aaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGt                                                <  1:3883524/71‑1 (MQ=255)
              aaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTa                                               >  1:2947631/1‑71 (MQ=255)
                 tttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGt                                                <  1:1239893/67‑1 (MQ=255)
                  ttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACgg                                            >  1:737801/1‑70 (MQ=255)
                   tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAAcca                                                                   >  1:96347/1‑45 (MQ=255)
                   tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                                    >  1:1112419/1‑61 (MQ=255)
                      ggggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGTCAAcc                                       <  1:3590309/71‑1 (MQ=255)
                        ggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCaa                                     <  1:2852499/71‑1 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAAcca                                                                   >  1:4041080/1‑36 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACTCACGCTACTACACAGGTACGGCAACCAATACa                                 >  1:508002/1‑71 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                               >  1:3961321/1‑57 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                                     >  1:4262611/1‑67 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                                     >  1:4462017/1‑67 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATAc                                  <  1:3326430/70‑1 (MQ=255)
                            aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa                                 >  1:1007127/1‑71 (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:966424/51‑1 (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:4513815/1‑48 (MQ=255)
                                  gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                                  >  1:1436344/1‑48 (MQ=255)
                                   cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTAc                                              >  1:2244745/1‑51 (MQ=255)
                                   cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTAc                                              >  1:2238228/1‑51 (MQ=255)
                                     caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                        <  1:3192130/68‑1 (MQ=255)
                                     caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                        <  1:1102946/68‑1 (MQ=255)
                                      aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAgg                                                 <  1:2021130/43‑1 (MQ=12)
                                      aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAgg                                                 <  1:1078367/43‑1 (MQ=12)
                                      aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgca                       >  1:649410/3‑70 (MQ=255)
                                      aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc                        <  1:2550428/68‑1 (MQ=255)
                                      aGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                 <  1:1912060/45‑1 (MQ=255)
                                       gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                                 <  1:1449875/44‑1 (MQ=255)
                                       gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                           >  1:4059290/1‑66 (MQ=255)
                                       gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                           >  1:133643/1‑66 (MQ=255)
                                       atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                      >  1:4393609/2‑69 (MQ=255)
                                       atCCTGAATCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat                      <  1:1165970/70‑3 (MQ=255)
                                        tCCTGATTTTCCTTATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                           <  1:3634267/49‑1 (MQ=255)
                                        tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                           <  1:5149652/49‑1 (MQ=255)
                                        tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                                           <  1:1982038/49‑1 (MQ=255)
                                        tCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACg                                             >  1:2740211/1‑47 (MQ=16)
                                          cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                   >  1:880577/1‑71 (MQ=255)
                                          cTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAt                                >  1:1394768/1‑58 (MQ=255)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:4594554/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:3370833/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:4574773/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:5273999/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:5239275/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:168476/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                               >  1:386318/1‑43 (MQ=12)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                                >  1:1124460/1‑58 (MQ=255)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                                >  1:2844205/1‑58 (MQ=255)
                                          cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg                   >  1:1523518/1‑69 (MQ=255)
                                           tGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                           <  1:3074446/62‑1 (MQ=255)
                                            gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                           >  1:3799546/1‑61 (MQ=255)
                                            gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt                           >  1:3681554/1‑61 (MQ=255)
                                            gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                  >  1:5298515/1‑70 (MQ=255)
                                            gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                                          <  1:4620319/46‑1 (MQ=16)
                                            gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                                          <  1:4354397/46‑1 (MQ=16)
                                            gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                                          <  1:2131427/46‑1 (MQ=16)
                                            gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc                  >  1:4850430/1‑70 (MQ=255)
                                             aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCa                                          >  1:1634268/1‑45 (MQ=255)
                                             aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCaa                >  1:5036326/1‑71 (MQ=255)
                                             aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                 >  1:3737640/1‑70 (MQ=255)
                                             aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                 <  1:2285167/70‑1 (MQ=255)
                                             aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAt                                >  1:2524375/1‑55 (MQ=255)
                                             aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg                         >  1:515577/1‑62 (MQ=255)
                                             aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc                  <  1:602642/69‑1 (MQ=255)
                                             aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc                  <  1:4796910/69‑1 (MQ=255)
                                             aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCa                 <  1:3064599/70‑1 (MQ=255)
                                                 tCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat            >  1:4631497/1‑71 (MQ=255)
                                                   atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCACCCGATACATTTCGTCGCATCCGCaataata           >  1:4067795/2‑69 (MQ=255)
                                                   atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaataata           <  1:5204141/69‑2 (MQ=255)
                                                   atgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAataa             <  1:2683741/67‑3 (MQ=255)
                                                     gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac        >  1:1356230/1‑70 (MQ=255)
                                                     gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac        >  1:3092677/1‑70 (MQ=255)
                                                      atgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAttt                              >  1:1321704/1‑48 (MQ=19)
                                                      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)
                                                      atgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAAc       >  1:2564220/1‑71 (MQ=255)
                                                       tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg      >  1:3193442/1‑71 (MQ=255)
                                                       tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg      >  1:1061345/1‑71 (MQ=255)
                                                        gagaAACCGCGCACGCTACTTCACAGGTACGGCaa                                         >  1:1746774/1‑35 (MQ=255)
                                                        gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat            >  1:4546756/1‑64 (MQ=255)
                                                        gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACgg     >  1:4285982/1‑71 (MQ=255)
                                                         agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt    >  1:4197328/1‑71 (MQ=255)
                                                          gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt    >  1:3715697/1‑70 (MQ=255)
                                                          gaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGaa         >  1:113282/1‑65 (MQ=255)
                                                          gaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt   <  1:4195368/71‑1 (MQ=255)
                                                           aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtc  <  1:1646897/71‑1 (MQ=255)
                                                           aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg                   <  1:5528802/54‑1 (MQ=255)
                                                             aCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa                 >  1:1485272/1‑54 (MQ=255)
                                                               |                                                                  
GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTC  >  minE/1826192‑1826321

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: