Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,472,315 +C intergenic (‑594/‑102) yfaL ← / → nrdA adhesin/ribonucleoside diphosphate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,472,3121.C100.0% 145.0 / NA 40intergenic (‑591/‑105)yfaL/nrdAadhesin/ribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (0/40);  total (0/40)

AAGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCAAACATACC  >  minE/1472264‑1472324
                                                 |            
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:60854/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:30054/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1049102/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:31526/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:359046/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:409935/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:411559/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:432801/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:575276/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:607212/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:313807/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:64470/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:658790/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:753141/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:807047/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:864246/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:867852/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:899890/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:971941/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:227888/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1085155/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1124144/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1125884/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1157068/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1164737/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:117004/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1175300/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1185101/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1190462/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1198928/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:146426/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:157979/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:172424/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:180736/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:219094/62‑1 (MQ=255)
aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:275371/62‑1 (MQ=255)
 aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:632443/61‑1 (MQ=255)
 aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:789427/61‑1 (MQ=255)
 aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:20906/61‑1 (MQ=255)
 aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:932360/61‑1 (MQ=255)
                                                 |            
AAGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCAAACATACC  >  minE/1472264‑1472324

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