Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,581,666 G→A 100% A523A (GCC→GCT ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,581,6660GA84.8% 40.0 / 0.2 26A523A (GCC→GCTligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (0/22);  total (0/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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: Variant not supported by required number of reads on each strand.

AGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGAGCAAAGGT  >  minE/1581612‑1581716
                                                      |                                                  
aGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACtt                                       <  1:1380634/68‑1 (MQ=255)
aGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACtt                                       <  1:2896921/68‑1 (MQ=255)
            cgcttccagcgTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAg                          <  1:709063/69‑1 (MQ=255)
               ttccagcgTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTg                        <  1:658826/68‑1 (MQ=255)
                                   ccgccAGACCTCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2722038/69‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:926913/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1059227/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:905658/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:627407/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:534656/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:430054/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:369726/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:3277937/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:3156848/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:3059079/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2701767/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2647867/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2221721/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2129867/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:2022302/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1884441/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1866148/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1446427/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1324574/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:1142545/68‑1 (MQ=255)
                                     gccAGACCTGCTGCGGTAGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGGGCAAAggt  <  1:113413/68‑1 (MQ=255)
                                                      |                                                  
AGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGAGCAAAGGT  >  minE/1581612‑1581716

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: