Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,824,288 G→A 8.1% D452D (GAC→GAT hypF ← carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,824,2880GA 8.1% 234.1 / 8.7 74D452D (GAC→GAThypFcarbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base G (29/39);  new base A (3/3);  total (32/42)
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.58e-01

ATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGC  >  minE/1824219‑1824349
                                                                     |                                                             
atgCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGt                                                              <  1:845121/71‑1 (MQ=255)
  gCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCa                                                            >  1:1109051/1‑71 (MQ=255)
  gCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCa                                                            >  1:720887/1‑71 (MQ=255)
   cAGTACCGTTTTCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCa                                                            <  1:4606304/70‑1 (MQ=255)
    aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATg                                                          >  1:2077876/1‑71 (MQ=255)
    aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATg                                                          >  1:162727/1‑71 (MQ=255)
       aCCGTTTGCGTAGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGca                                                        >  1:2617037/1‑70 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcaga                                                      <  1:1331171/71‑3 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcaca                                                      <  1:3210985/71‑1 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcaca                                                      <  1:2415235/71‑1 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcac                                                       <  1:1376538/70‑1 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcac                                                       <  1:230331/70‑1 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGcac                                                       <  1:1158461/70‑1 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGca                                                        >  1:1153447/1‑69 (MQ=255)
        ccGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGca                                                        >  1:4301469/1‑69 (MQ=255)
             tGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGt                                                  <  1:2299363/70‑1 (MQ=255)
             tGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGt                                                  <  1:4631037/70‑1 (MQ=255)
               cGTCGGCAGATTCATTTCACGCGCCCACTGGCTTGAGACATAGCCCGGATGCGCGTCATGCACAACGTAt                                                <  1:1767945/70‑1 (MQ=255)
               cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTAtt                                               >  1:632069/1‑71 (MQ=255)
               cGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCAGGATGCGCGTCATGCACAACGTAtt                                               <  1:198718/71‑1 (MQ=255)
                 tCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGc                                             <  1:4398569/71‑1 (MQ=255)
                 tCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGc                                             <  1:4244354/71‑1 (MQ=255)
                 tCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGc                                             <  1:3046037/71‑1 (MQ=255)
                   ggCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCgg                                           >  1:3083874/1‑71 (MQ=255)
                   ggCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCg                                            <  1:1039020/70‑1 (MQ=255)
                        aTTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTa                                       >  1:1450710/1‑70 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCg                                            <  1:2928106/63‑1 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCg                                            <  1:3287459/63‑1 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCg                                            <  1:3219127/63‑1 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCg                                            <  1:4409102/63‑1 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTaaaa                                    >  1:3745033/1‑71 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTaa                                      <  1:2288248/69‑1 (MQ=255)
                          tCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTaa                                      <  1:3126627/69‑1 (MQ=255)
                           cATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGt                                                  <  1:2465780/56‑1 (MQ=255)
                           cATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGt                                                  <  1:4502781/56‑1 (MQ=255)
                               tCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCg                                 >  1:4366508/1‑69 (MQ=255)
                                cACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                 >  1:113487/1‑68 (MQ=255)
                                  cgcgCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                            >  1:2231729/1‑71 (MQ=255)
                                   gcgcCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATg                           <  1:4426735/71‑1 (MQ=255)
                                   gcgcCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATg                           <  1:3321476/71‑1 (MQ=255)
                                    cgcCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGt                          <  1:4258615/71‑1 (MQ=255)
                                      cccACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTaaaa                                    >  1:3580727/1‑59 (MQ=255)
                                       ccACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                       >  1:4468023/1‑71 (MQ=255)
                                         aCTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                      >  1:1496647/1‑70 (MQ=255)
                                         aCTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                      >  1:4745326/1‑70 (MQ=255)
                                          gtggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATg                           <  1:893082/63‑1 (MQ=255)
                                            ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                   >  1:3746463/1‑70 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                   <  1:4662884/69‑1 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                   >  1:3049647/1‑69 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                   <  1:2591730/69‑1 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                  <  1:1175588/70‑1 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                  <  1:2062473/70‑1 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                  <  1:397388/70‑1 (MQ=255)
                                             gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATc                  <  1:2118032/70‑1 (MQ=255)
                                                 gagaCATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCg             >  1:3802699/1‑71 (MQ=255)
                                                 gagaCATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCg             >  1:3983930/1‑71 (MQ=255)
                                                    aCATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTaa          >  1:3111742/1‑71 (MQ=255)
                                                    aCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTTGATGttt                        <  1:3424063/57‑1 (MQ=255)
                                                    aCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                        <  1:2575433/57‑1 (MQ=255)
                                                    aCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                        <  1:3171997/57‑1 (MQ=255)
                                                     cATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGCTGtt                         >  1:1491502/1‑55 (MQ=255)
                                                     cATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGt                          <  1:2495351/54‑1 (MQ=255)
                                                       tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                        >  1:2758085/1‑54 (MQ=255)
                                                       tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                        >  1:2952266/1‑54 (MQ=255)
                                                       tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                        >  1:571452/1‑54 (MQ=255)
                                                       tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGc       <  1:2082073/71‑1 (MQ=255)
                                                         gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAg                              >  1:1174192/1‑46 (MQ=255)
                                                         gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGa                             >  1:4271988/1‑47 (MQ=255)
                                                         gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGa                             >  1:2965329/1‑47 (MQ=255)
                                                            cGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTcgc  <  1:4310852/71‑1 (MQ=255)
                                                            cGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTcgc  <  1:2734229/71‑1 (MQ=255)
                                                             ggATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                            >  1:108199/1‑44 (MQ=255)
                                                             ggATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                            >  1:1213881/1‑44 (MQ=255)
                                                                tGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGt                          <  1:2852777/43‑1 (MQ=255)
                                                                  cgcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt            >  1:2046531/1‑55 (MQ=255)
                                                                  cgcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt            >  1:867957/1‑55 (MQ=255)
                                                                     gTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGc       <  1:839674/57‑1 (MQ=255)
                                                                     |                                                             
ATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGC  >  minE/1824219‑1824349

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: