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% 67.2 / NA 29A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/7);  total (22/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGC  >  minE/1036445‑1036570
                                                            |                                                                 
agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGTGCCTAAGCACGCAAGc                                                            >  1:1777683/1‑68 (MQ=11)
 gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                         >  1:2933544/1‑70 (MQ=25)
 gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                         >  1:3476829/1‑70 (MQ=25)
 gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                         >  1:3256425/1‑70 (MQ=25)
  aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                       >  1:748096/1‑71 (MQ=21)
                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                     <  1:1287536/71‑1 (MQ=14)
                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                     >  1:2348974/1‑71 (MQ=14)
                       aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                  <  1:3145025/71‑1 (MQ=25)
                              cttCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                            >  1:165538/3‑70 (MQ=255)
                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                            >  1:1486325/1‑70 (MQ=255)
                                tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                         >  1:845794/1‑71 (MQ=255)
                                tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                         <  1:1290464/71‑1 (MQ=255)
                                  gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                         >  1:1990168/1‑69 (MQ=255)
                                  gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                         >  1:2178448/1‑69 (MQ=255)
                                  gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                         >  1:175059/1‑69 (MQ=255)
                                     gCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                    <  1:206580/71‑1 (MQ=255)
                                          tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                     <  1:2134332/49‑1 (MQ=255)
                                           aTTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc              >  1:975858/1‑71 (MQ=255)
                                           aTTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc              >  1:1738303/1‑71 (MQ=255)
                                            ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                  <  1:913465/66‑1 (MQ=255)
                                                 cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                   >  1:2598176/1‑60 (MQ=255)
                                                 cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCgcgc          >  1:2627698/1‑69 (MQ=255)
                                                 cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa         >  1:279810/1‑70 (MQ=255)
                                                 cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa         >  1:2418075/1‑70 (MQ=255)
                                                    gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc     >  1:1519599/1‑71 (MQ=255)
                                                     ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc    >  1:1269659/1‑71 (MQ=255)
                                                     ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc     >  1:1086865/1‑70 (MQ=255)
                                                       aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                       >  1:1236808/1‑50 (MQ=255)
                                                       aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc  <  1:3290844/71‑1 (MQ=255)
                                                            |                                                                 
AGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGC  >  minE/1036445‑1036570

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: