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

CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCC  >  minE/1036438‑1036568
                                                                   |                                                               
cGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                               >  1:3342904/1‑70 (MQ=25)
     ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                                         <  1:244050/71‑1 (MQ=21)
           gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                     >  1:2219579/1‑69 (MQ=14)
           gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                     >  1:3090483/1‑69 (MQ=14)
           gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                     >  1:2215140/1‑69 (MQ=14)
                           gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAAAAGGCTAATGcc                                   <  1:544677/71‑1 (MQ=9)
                                         gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGGGCCAGaacaa                       >  1:4034404/1‑69 (MQ=255)
                                         gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                       >  1:2380731/1‑69 (MQ=255)
                                          aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                    <  1:1855895/71‑1 (MQ=255)
                                            gCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                        <  1:3702383/65‑1 (MQ=255)
                                              cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                <  1:1618434/71‑1 (MQ=255)
                                                cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt              >  1:2075116/1‑71 (MQ=255)
                                                cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt              >  1:2307219/1‑71 (MQ=255)
                                                cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt              >  1:3025922/1‑71 (MQ=255)
                                                 tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                          >  1:1925555/1‑58 (MQ=255)
                                                 tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                          >  1:525494/1‑58 (MQ=255)
                                                   ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaac                      <  1:193619/60‑1 (MQ=255)
                                                   ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaac                      <  1:1799812/60‑1 (MQ=255)
                                                   ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc            <  1:3878506/70‑1 (MQ=255)
                                                    tGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc          >  1:147071/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg         >  1:3448730/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg         >  1:772001/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc          <  1:370080/70‑1 (MQ=255)
                                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc            >  1:1110896/1‑62 (MQ=255)
                                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc             >  1:2839008/1‑61 (MQ=255)
                                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc             >  1:2834421/1‑61 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                     >  1:2169879/1‑52 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc  >  1:287327/1‑71 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc  >  1:3227180/1‑71 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc  >  1:2696879/1‑71 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc  >  1:3767166/1‑71 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgcc  <  1:214834/71‑1 (MQ=255)
                                                              aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                       >  1:1529234/1‑48 (MQ=255)
                                                                aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                       >  1:1686034/1‑46 (MQ=255)
                                                                   |                                                               
CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCC  >  minE/1036438‑1036568

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: