Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 505,424 T→C G128G (GGT→GGC ybhE → 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE505,4240TC100.0% 139.2 / NA 38G128G (GGT→GGCybhE6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/33);  total (5/33)

CGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCG  >  minE/505366‑505489
                                                          |                                                                 
cgcgTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:4917966/71‑1 (MQ=255)
cgcgTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:4662851/71‑1 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:3400886/65‑1 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:4315342/65‑1 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:4031363/65‑1 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:2294248/65‑1 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtc                                                 >  1:838281/1‑71 (MQ=255)
      tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtc                                                 >  1:3577657/1‑71 (MQ=255)
       ggAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGcc                                                       <  1:3029947/64‑1 (MQ=255)
        taaGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1818357/70‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1856577/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:2907878/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:2240042/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:4097800/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1990072/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:4154927/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:728170/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1358738/71‑1 (MQ=255)
        gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1005542/71‑1 (MQ=255)
              ggCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATAtctc                                               <  1:1247236/65‑1 (MQ=255)
                ccTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAat                                                    <  1:1601017/58‑1 (MQ=255)
                  tGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTcc                                                          >  1:3205062/1‑50 (MQ=255)
                  tGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAatat                                                  >  1:716091/1‑58 (MQ=255)
                  tGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAatat                                                  >  1:4063910/1‑58 (MQ=255)
                    ccAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCg                                   <  1:1786657/71‑1 (MQ=255)
                            cgtcgtCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                           <  1:3673094/71‑1 (MQ=255)
                            cgtcgtCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                           <  1:60895/71‑1 (MQ=255)
                            cgtcgtCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                           <  1:2010326/71‑1 (MQ=255)
                            cgtcgtCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                           <  1:4155230/71‑1 (MQ=255)
                                gtcgATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGt                       <  1:1308977/71‑1 (MQ=255)
                                gtcgATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGt                       <  1:3239708/71‑1 (MQ=255)
                                          tCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa             <  1:431456/71‑1 (MQ=255)
                                          tCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa             <  1:1802508/71‑1 (MQ=255)
                                          tCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa             <  1:1633096/71‑1 (MQ=255)
                                          tCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa             <  1:2521539/71‑1 (MQ=255)
                                          tCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa             <  1:3067738/71‑1 (MQ=255)
                                                     gACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCg  <  1:1809295/71‑1 (MQ=255)
                                                     gACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCg  <  1:2533347/71‑1 (MQ=255)
                                                          |                                                                 
CGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCG  >  minE/505366‑505489

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: