Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,036,505 G→A 100% A632A (GCG→GCA rsxC → fused predicted 4Fe‑4S ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,036,5050GA100.0% 100.7 / NA 36A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/17);  total (19/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036441‑1036572
                                                                |                                                                   
aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                               >  1:2274923/1‑71 (MQ=21)
 aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacgc                                                                  <  1:1037891/67‑1 (MQ=25)
    agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                                             <  1:3219208/69‑1 (MQ=25)
                        gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       >  1:2252643/1‑71 (MQ=14)
                         cAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCg                                      <  1:83373/71‑1 (MQ=17)
                            aGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGc                                   >  1:3426415/1‑71 (MQ=255)
                               cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                                <  1:3782626/71‑1 (MQ=255)
                                 ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               <  1:1968734/70‑1 (MQ=255)
                                  cgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:1531595/71‑1 (MQ=255)
                                  cgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                              >  1:2790062/1‑70 (MQ=255)
                                         gCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                      >  1:3079773/1‑71 (MQ=255)
                                          ccGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                     >  1:1468238/1‑71 (MQ=255)
                                           cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:1258739/1‑55 (MQ=255)
                                           cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:1473371/71‑1 (MQ=255)
                                           cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:1723424/71‑1 (MQ=255)
                                            gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   >  1:3765531/1‑71 (MQ=255)
                                                ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:960419/70‑1 (MQ=255)
                                                  gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             >  1:595214/1‑71 (MQ=255)
                                                  gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc              >  1:2973990/1‑70 (MQ=255)
                                                  gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc              >  1:2978573/1‑70 (MQ=255)
                                                   cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa                                 >  1:2643875/1‑50 (MQ=255)
                                                    ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa           <  1:185730/71‑1 (MQ=255)
                                                     cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:1167438/71‑1 (MQ=255)
                                                      gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                           <  1:1509975/53‑1 (MQ=255)
                                                        gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:483527/49‑1 (MQ=255)
                                                         ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       <  1:2917696/70‑1 (MQ=255)
                                                         ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCAAGAACAACAGGTCGATCCGCGCAAAgcc      >  1:2580786/1‑71 (MQ=255)
                                                          cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:1763999/1‑45 (MQ=255)
                                                          cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:1484028/1‑45 (MQ=255)
                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                                <  1:1482823/43‑1 (MQ=255)
                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc    >  1:1052873/1‑71 (MQ=255)
                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc    >  1:3049128/1‑71 (MQ=255)
                                                            aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:1428698/46‑1 (MQ=255)
                                                            aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:1944937/46‑1 (MQ=255)
                                                             aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc      >  1:2590476/1‑67 (MQ=255)
                                                             aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  <  1:1849377/71‑1 (MQ=255)
                                                             aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  >  1:2329866/1‑71 (MQ=255)
                                                                |                                                                   
AACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036441‑1036572

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: