Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,858,845 G→A 100% S373S (TCG→TCA amiB → N‑acetylmuramoyl‑l‑alanine amidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,858,8450GA100.0% 164.6 / NA 52S373S (TCG→TCAamiBN‑acetylmuramoyl‑l‑alanine amidase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (45/7);  total (45/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTT  >  minE/2858794‑2858881
                                                   |                                    
cGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATc                        <  1:1115413/66‑1 (MQ=255)
cGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATc                        <  1:927571/66‑1 (MQ=255)
cGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATc                        <  1:1835370/66‑1 (MQ=255)
                   ccAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGt     <  1:1891372/66‑1 (MQ=255)
                   ccAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGt     <  1:1506016/66‑1 (MQ=255)
                   ccAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGt     <  1:1470841/66‑1 (MQ=255)
                   ccAGAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGt     <  1:616423/66‑1 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATTCCATCAGTACTGGTCGAAACCGGtttt  >  1:1264448/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACt                  >  1:133624/1‑50 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAAc         >  1:2112085/1‑59 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:56446/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:2027091/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1986198/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1949954/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1918623/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:2060955/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:212933/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:345470/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:366508/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:476035/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:479846/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:498826/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1855568/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:622575/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:687887/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:72062/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:739979/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:957114/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:990780/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1583238/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1014689/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1104932/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1138199/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1159601/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1212855/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1327944/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1337423/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1377116/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1518522/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:156059/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1573150/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1865952/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:160266/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1615956/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:167438/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1722265/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1797612/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1834828/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1836002/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGtttt  >  1:1014114/1‑66 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGGTCGAAACCGGttt   >  1:192778/1‑65 (MQ=255)
                      gAACACGCCAGCCTTGGCGTTCTGCGCTCACCGGATATCCCATCAGTACTGCTCGAAACCGGtttt  >  1:471141/1‑66 (MQ=255)
                                                   |                                    
CGAAATACATAAACGTCGACCAGAACACGCCAGCCTTGGCGTTCTGCGCTCGCCGGATATCCCATCAGTACTGGTCGAAACCGGTTTT  >  minE/2858794‑2858881

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: