Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,897,915:1 +A 100% coding (747/1374 nt) mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,897,9151.A94.7% 122.5 / 0.7 38coding (747/1374 nt)mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (0/36);  total (0/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGC  >  minE/2897872‑2897955
                                            |                                        
ggAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTc                   <  1:1632396/67‑1 (MQ=255)
ggAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGAGAAAAAGCTGACCACCGATGCTTc                   <  1:1940762/67‑1 (MQ=255)
                   gggtgAGCAGGGTGACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2243689/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:556608/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2270709/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2302903/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:326972/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:416403/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:426120/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:472414/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:509253/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:552569/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2213807/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:611832/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:649031/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:763305/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:836160/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:859767/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:905314/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:954483/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:973208/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1835551/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1335170/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1373002/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1461813/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1462239/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:152082/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1710327/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1778931/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1811443/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1022229/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1921931/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:1972711/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2047717/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2059765/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2156856/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:219510/66‑1 (MQ=255)
                   gggtgAGCAGGGTCACTGGCAGGCGAAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGc  <  1:2199717/66‑1 (MQ=255)
                                            |                                        
GGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGC  >  minE/2897872‑2897955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: