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,5050GA98.5% 98.8 / ‑6.1 36A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/15);  minor base T (0/1);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCC  >  minE/1036438‑1036571
                                                                   |                                                                  
cGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                  <  1:2701650/70‑1 (MQ=25)
 ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCa                                                                <  1:4036829/71‑1 (MQ=21)
    aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGc                                                             <  1:3962325/71‑1 (MQ=21)
       agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTg                                                           >  1:2971616/1‑70 (MQ=25)
                                  cgccgGCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                               <  1:230443/71‑1 (MQ=255)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                    >  1:1218509/1‑65 (MQ=25)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                    >  1:317251/1‑65 (MQ=25)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                    >  1:2995833/1‑65 (MQ=25)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                    >  1:2044296/1‑65 (MQ=25)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                    >  1:2339135/1‑65 (MQ=25)
                                       tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                          >  1:2702210/1‑71 (MQ=255)
                                         gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                        <  1:2133574/71‑1 (MQ=255)
                                          aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                       <  1:922999/71‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                  <  1:1277736/71‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                  <  1:551810/71‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                  <  1:1687438/71‑1 (MQ=255)
                                                cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                     >  1:3047006/1‑67 (MQ=255)
                                                cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                     >  1:2957731/1‑67 (MQ=255)
                                                   ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc               <  1:501150/70‑1 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg            >  1:1761236/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc             <  1:1694660/70‑1 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc             <  1:2159575/70‑1 (MQ=255)
                                                        cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa         <  1:3063562/71‑1 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                 >  1:2086710/1‑62 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                 >  1:2648241/1‑62 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                 >  1:608168/1‑62 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaaa        <  1:842776/71‑1 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaaa        <  1:2485706/71‑1 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa          >  1:2494944/1‑69 (MQ=255)
                                                          tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                          <  1:3780606/52‑1 (MQ=255)
                                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcg              >  1:1825534/1‑63 (MQ=255)
                                                           gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc      >  1:1902017/1‑71 (MQ=255)
                                                             cAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                              <  1:150018/45‑1 (MQ=255)
                                                               aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc  >  1:2808960/1‑71 (MQ=255)
                                                               aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc  >  1:3485178/1‑71 (MQ=255)
                                                               aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc  >  1:3680434/1‑71 (MQ=255)
                                                               aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc  >  1:401656/1‑71 (MQ=255)
                                                                   |                                                                  
CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCC  >  minE/1036438‑1036571

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: