Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,007,627 T→A 5.5% S92C (AGC→TGC)  hybA ← hydrogenase 2 4Fe‑4S ferredoxin‑type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,007,6270TA 5.5% 387.5 / 10.6 128S92C (AGC→TGC) hybAhydrogenase 2 4Fe‑4S ferredoxin‑type component
Reads supporting (aligned to +/- strand):  ref base T (73/48);  new base A (4/3);  total (77/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 129 total aligned reads displayed.

GCATACACTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTGTTC  >  minE/2007561‑2007693
                                                                  |                                                                  
gCATACACTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                >  1:2740212/1‑71 (MQ=255)
 cATACACTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                >  1:796900/1‑70 (MQ=255)
 cATACACTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                >  1:700498/1‑70 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTgg                                                                  >  1:208589/1‑64 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTgg                                                                  >  1:453004/1‑64 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCcaca                                                            <  1:1659095/70‑1 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCcaca                                                            <  1:4072505/70‑1 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCacac                                                           >  1:5180645/1‑71 (MQ=255)
     cacTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCacac                                                           >  1:4990902/1‑71 (MQ=255)
       cTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                <  1:4169036/64‑1 (MQ=255)
       cTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                <  1:2455171/64‑1 (MQ=255)
       cTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTc                                                                <  1:3834793/64‑1 (MQ=255)
         gTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCca                                                              >  1:123257/1‑64 (MQ=255)
         gTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCca                                                              >  1:583697/1‑64 (MQ=255)
            tCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAat                                                     >  1:1922696/1‑70 (MQ=255)
            tCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAat                                                     <  1:3679890/70‑1 (MQ=255)
            tCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAat                                                     <  1:2755543/70‑1 (MQ=255)
                    tACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTAtt                                            >  1:3777984/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:86870/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:4414956/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:3768433/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:715364/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:728938/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:2288846/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTg                                           >  1:1253572/1‑71 (MQ=255)
                     aCGCGTAGCCGTTCTCCTCCTGGTCTTTCTTGACCCCTGTGCCGCt                                                                    >  1:2390552/1‑46 (MQ=255)
                      cgcgTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTgg                                          <  1:1532906/71‑1 (MQ=255)
                         gTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        >  1:4398538/1‑70 (MQ=255)
                             ccGTTCTCCTCCTGGTCTTTGATGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACg                                    <  1:511323/70‑1 (MQ=255)
                                ttCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCca                                                              >  1:3652396/1‑41 (MQ=255)
                                ttCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGa                                >  1:4503022/1‑71 (MQ=255)
                                ttCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGa                                >  1:3639473/1‑71 (MQ=255)
                                 tctcCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGa                                >  1:4575679/1‑70 (MQ=255)
                                  ctcctcCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCg                                 >  1:176224/1‑68 (MQ=255)
                                  ctcctcCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCg                                 >  1:5135054/1‑68 (MQ=255)
                                    cctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGt                            >  1:5052280/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCca                                                              >  1:4625239/1‑36 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGtt                           >  1:2664443/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGtt                           >  1:3974240/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGtt                           >  1:4129492/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGtt                           >  1:2120948/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGtt                           >  1:2929327/1‑71 (MQ=255)
                                     ctcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATAGGCGACAGttt                          >  1:4600148/1‑71 (MQ=255)
                                      tcctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGttt                          >  1:2932556/1‑71 (MQ=255)
                                       cctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTg                         >  1:3259800/1‑71 (MQ=255)
                                        ctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGttt                          <  1:3236423/69‑1 (MQ=255)
                                        ctGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGttt                          <  1:5078841/69‑1 (MQ=255)
                                         tGGTCTTTGTTGACCCCTGTGCCGCTGGTCCacac                                                           >  1:1715503/1‑35 (MQ=255)
                                         tGGTCTTTGTTGACCCCTGTGCCGCTGGTCCacac                                                           >  1:2150741/1‑35 (MQ=255)
                                         tGGTCTTTGTTGACCCCTGTGCCGCTGGTCCacac                                                           >  1:2518740/1‑35 (MQ=255)
                                          ggTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGttt                          <  1:67874/67‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTAtt                                            >  1:3772947/1‑45 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:3984165/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:4428008/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:4280687/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:4296999/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:2069837/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:3629422/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:1146408/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGta                                        <  1:2251694/49‑1 (MQ=255)
                                              tttGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCgttgt                  >  1:2941235/1‑71 (MQ=255)
                                                tgttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCgttgtt                 >  1:3548612/1‑70 (MQ=255)
                                                tgttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                >  1:5169941/1‑71 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                >  1:479860/1‑70 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                <  1:4936254/70‑1 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                >  1:1726420/1‑70 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                >  1:3031139/1‑70 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                <  1:3705081/70‑1 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                <  1:2261781/70‑1 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                <  1:3018262/70‑1 (MQ=255)
                                                 gttgACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTc                <  1:4065921/70‑1 (MQ=255)
                                                        ccTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTc                       <  1:1033773/56‑1 (MQ=255)
                                                        ccTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTc                       <  1:1053524/56‑1 (MQ=255)
                                                        ccTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTc                       <  1:1042796/56‑1 (MQ=255)
                                                        ccTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTc                       <  1:4423047/56‑1 (MQ=255)
                                                         cTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTc       >  1:628342/1‑71 (MQ=255)
                                                          tgtgCCGCTGGTCCACACCTGAATGATGTTATTgg                                          <  1:3668394/35‑1 (MQ=255)
                                                          tgtgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      >  1:5150588/1‑71 (MQ=255)
                                                          tgtgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      >  1:3275087/1‑71 (MQ=255)
                                                          tgtgCCGCAGGTCCACACCTGAATGATGTTATTGGt                                         <  1:1890643/36‑1 (MQ=255)
                                                          tgtgCCGCAGGTCCACACCTGAATGATGTTATTGGt                                         <  1:860174/36‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      <  1:1183958/69‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      <  1:3528751/69‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      <  1:4287267/69‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      <  1:184521/69‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCt      <  1:1310209/69‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:3339384/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     >  1:187873/1‑70 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:605860/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:3069636/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:3605759/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:3290397/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTg     <  1:4342753/70‑1 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTGt    >  1:3691625/1‑71 (MQ=255)
                                                            tgCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTGt    >  1:4343722/1‑71 (MQ=255)
                                                            tgCCGCAGGTCCACACCTGAATGATGTTATTGGtat                                       >  1:471527/1‑36 (MQ=255)
                                                            tgCCGCAGGTCCACACCTGAATGATGTTATTGGtat                                       >  1:2770477/1‑36 (MQ=255)
                                                            tgCCGCAGGTCCACACCTGAATGATGTTATTGGtat                                       >  1:1085712/1‑36 (MQ=255)
                                                              ccGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTGTTc  >  1:4540990/1‑71 (MQ=255)
                                                                 cTGGTCCACACCTGAATGATGTAATTTGTATACGGCGACAGTTTGTc                       <  1:4920526/47‑1 (MQ=39)
                                                                  |                                                                  
GCATACACTGTTTCTTAATGTACGCGTAGCCGTTCTCCTCCTGGTCTTTGTTGACCCCTGTGCCGCTGGTCCACACCTGAATGATGTTATTGGTATACGGCGACAGTTTGTCGTTGTTCGACCAGGTCTGTTC  >  minE/2007561‑2007693

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: