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,5050GA96.0% 63.2 / ‑3.3 25A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/6);  total (18/7)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036437‑1036572
                                                                    |                                                                   
gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:2655948/1‑71 (MQ=21)
    aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                               <  1:2318068/71‑1 (MQ=34)
     aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGc                                                              >  1:2663386/1‑71 (MQ=21)
      ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   <  1:978317/65‑1 (MQ=25)
       cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTg                                                            >  1:2441878/1‑71 (MQ=21)
        agaagtACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                                                           <  1:2637761/71‑1 (MQ=255)
          aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                         >  1:1773035/1‑71 (MQ=21)
            gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                         >  1:714160/1‑69 (MQ=14)
                                      cgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:316779/71‑1 (MQ=255)
                                          gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                           >  1:2617595/1‑69 (MQ=255)
                                          gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                           >  1:369504/1‑69 (MQ=255)
                                          gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                           >  1:2212434/1‑69 (MQ=255)
                                               cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:1362665/71‑1 (MQ=255)
                                                      gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             >  1:524159/1‑71 (MQ=255)
                                                      gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             >  1:873716/1‑71 (MQ=255)
                                                          gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                              >  1:1417007/1‑50 (MQ=255)
                                                           tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg        >  1:2595571/1‑71 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       >  1:2460743/1‑71 (MQ=255)
                                                             ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             >  1:48219/1‑48 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               <  1:2428524/45‑1 (MQ=255)
                                                               aaaGCACGCAAGCTGGAACAGGAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   <  1:1274984/56‑1 (MQ=14)
                                                               aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc              >  1:770345/1‑61 (MQ=255)
                                                               aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccg     >  1:2032220/1‑70 (MQ=255)
                                                                 aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:2198109/44‑1 (MQ=255)
                                                                 aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGCACAACAGGTCGATCCGCGCAAAGccgccg  >  1:2274208/1‑71 (MQ=255)
                                                                 aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  >  1:2933068/1‑71 (MQ=255)
                                                                 aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  >  1:2959899/1‑71 (MQ=255)
                                                                    |                                                                   
GCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036437‑1036572

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: