Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,330,586 A→G Y334H (TAC→CAC)  cpsG ← phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,330,5860AG100.0% 166.1 / NA 44Y334H (TAC→CAC) cpsGphosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/31);  total (13/31)

ACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTT  >  minE/1330522‑1330642
                                                                |                                                        
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:110420/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:92175/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:659648/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:4118295/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:3393807/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:3228423/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:2995200/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:2798150/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:1125220/71‑1 (MQ=255)
acCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGAtggtg                                                    <  1:1453081/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:835522/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:589606/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:538808/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:3583347/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:2420789/71‑1 (MQ=255)
    gTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGgcg                                                <  1:1505248/71‑1 (MQ=255)
                 cagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGccac                                   <  1:4050042/71‑1 (MQ=255)
                 cagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGccac                                   <  1:3972550/71‑1 (MQ=255)
                 cagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGccac                                   <  1:4327311/71‑1 (MQ=255)
                 cagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGccac                                   <  1:2129553/71‑1 (MQ=255)
                     cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:160346/71‑1 (MQ=255)
                     cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:2367447/71‑1 (MQ=255)
                     cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:3483988/71‑1 (MQ=255)
                     cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:1499432/71‑1 (MQ=255)
                     cACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:3472198/71‑1 (MQ=255)
                              aTGCCGCTGTCGCAGTAAGCGAAATCACTGAAGTGAtggtg                                                    <  1:1047926/41‑1 (MQ=38)
                                            gTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATa                               <  1:4262708/48‑1 (MQ=255)
                                            gTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATAc        >  1:113996/1‑71 (MQ=255)
                                             tAAGCGAAATCACGTAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACg       <  1:2653178/71‑1 (MQ=255)
                                             tAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACg       <  1:1967048/71‑1 (MQ=255)
                                             tAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACg       <  1:2708622/71‑1 (MQ=255)
                                             tAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACg       <  1:397612/71‑1 (MQ=255)
                                                 cGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt   >  1:373178/1‑71 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCCCCATAGATGGCGTCTTCCTTGCGCATACGtt     >  1:2067014/1‑67 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCtt                    >  1:3162627/1‑52 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTcc                  >  1:3690887/1‑54 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTgcg             >  1:3035634/1‑59 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTgc              >  1:3419601/1‑58 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt  >  1:434340/1‑70 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt  >  1:2326734/1‑70 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt  >  1:623966/1‑70 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt  >  1:1933229/1‑70 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt  >  1:1862053/1‑70 (MQ=255)
                                                   aaaTCACGGAAGTGATGGTGGGCGCTCATATCGCCACCAt                                >  1:2635342/1‑40 (MQ=38)
                                                                |                                                        
ACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTT  >  minE/1330522‑1330642

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: