Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,364,103 A→G W456R (TGG→CGG)  trkH ← potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,364,1030AG100.0% 229.6 / NA 64W456R (TGG→CGG) trkHpotassium transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/35);  total (29/35)

GTAAAGAGCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCATTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGCCCCAGGTTAT  >  minE/2364039‑2364168
                                                                |                                                                 
gTAAAGAGCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGc                                                             >  1:363429/1‑71 (MQ=255)
   aaGAGCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCAc                                                          >  1:1907670/1‑71 (MQ=255)
   aaGAGCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCAc                                                          >  1:3026853/1‑71 (MQ=255)
       gCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCgg                                                       >  1:65271/1‑70 (MQ=255)
        cACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGt                                                     <  1:1802840/71‑1 (MQ=255)
          ccAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCggg                                                      >  1:2839930/1‑68 (MQ=255)
           cAGCATTGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCTGGATCCGTTTAGCCACCGGGTTCa                                                  <  1:444492/71‑1 (MQ=255)
           cAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCa                                                  <  1:1365987/71‑1 (MQ=255)
             gCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTa                                                               >  1:693163/1‑56 (MQ=255)
              cAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATAc                                               >  1:1081831/1‑71 (MQ=255)
              cAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATAc                                               >  1:1929473/1‑71 (MQ=255)
                  gTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTgg                                            >  1:1904041/1‑70 (MQ=255)
                   tAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGTATCCGTTTAGCCACCGGGTTCATACTGGTa                                          <  1:2660718/71‑1 (MQ=255)
                   tAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGc                                                             >  1:2478966/1‑52 (MQ=255)
                   tAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGt                                           >  1:89762/1‑70 (MQ=255)
                    aaaGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTa                                          <  1:1909740/70‑1 (MQ=255)
                     aaGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGtt                                                    <  1:922689/59‑1 (MQ=255)
                     aaGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGtt                                                    <  1:1231318/59‑1 (MQ=255)
                      aGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATAc                                               <  1:1361812/63‑1 (MQ=255)
                      aGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATAc                                               <  1:1819585/63‑1 (MQ=255)
                      aGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAg                                       <  1:1100819/71‑1 (MQ=255)
                        aCCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGtt                                     <  1:625955/71‑1 (MQ=255)
                         ccTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTa                                    <  1:2435750/71‑1 (MQ=255)
                         ccTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTa                                    >  1:2632621/1‑71 (MQ=255)
                         ccTCGAGACGACAAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTa                                    >  1:74314/1‑71 (MQ=255)
                          cTCGAGACTACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTAt                                   <  1:1082955/71‑1 (MQ=255)
                          cTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCgg                                                       <  1:1683739/51‑1 (MQ=255)
                          cTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTa                                    <  1:961510/70‑1 (MQ=255)
                          cTCGAGACGACCAAACAGCATGATGGCAATCAGGATCCGTTTAGCCACCgg                                                       <  1:365753/51‑1 (MQ=39)
                          cTCGAGAAGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTAt                                   >  1:208008/1‑71 (MQ=255)
                           tCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTgg                                            <  1:770110/61‑1 (MQ=255)
                           tCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTAt                                   <  1:687994/70‑1 (MQ=255)
                            cGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGtt                                                                 >  1:489950/1‑39 (MQ=255)
                            cGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTgg                                            >  1:306522/1‑60 (MQ=255)
                            cGAGACGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTAt                                   >  1:2133223/1‑69 (MQ=255)
                             gagaCGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGt                                           >  1:958647/1‑60 (MQ=255)
                             gagaCGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTAt                                   >  1:2588321/1‑68 (MQ=255)
                              agaCGAGCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAg                                >  1:393877/1‑70 (MQ=255)
                              agaCGACCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGc                               <  1:1329107/71‑1 (MQ=255)
                                acgacCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCgg                                                       <  1:1209594/45‑1 (MQ=255)
                                 cgacCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAATTATCAGCAAc                            <  1:2299089/71‑1 (MQ=255)
                                  gacCAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCa                                                  <  1:2422837/48‑1 (MQ=255)
                                     cAAACAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGtt                                     <  1:2914412/58‑1 (MQ=255)
                                       aaCAGCATGTTGGCAATCAGGATCCGTTTAGCCAc                                                          >  1:1552860/1‑35 (MQ=255)
                                       aaCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCa                                                  <  1:2131037/43‑1 (MQ=255)
                                       aaCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCa                                                  <  1:127214/43‑1 (MQ=255)
                                       aaCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTaaa                                        >  1:1273487/1‑53 (MQ=255)
                                       aaCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGc                       >  1:1412182/1‑70 (MQ=255)
                                        aCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCa                                 >  1:2595376/1‑59 (MQ=255)
                                        aCAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCa                     <  1:857116/71‑1 (MQ=255)
                                         cAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCaa                    <  1:2087146/71‑1 (MQ=255)
                                         cAGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCaa                    <  1:1837527/71‑1 (MQ=255)
                                          aGCATGTTGGCAATCAGGATCCGTTTAGCCACCgg                                                       >  1:1694438/1‑35 (MQ=255)
                                          aGCATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCAAACTGGTAAAGTTATCAGCAACCACGCCAAg                   >  1:473996/1‑71 (MQ=255)
                                           gCATGTTGGCTTTCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGtt                                     <  1:904986/52‑1 (MQ=255)
                                            cATGTTGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGcc                 <  1:483827/71‑1 (MQ=255)
                                                ttGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTa                                          <  1:1261787/42‑1 (MQ=255)
                                                ttGGCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTgg             >  1:1861097/1‑71 (MQ=255)
                                                   gCAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGcc                      >  1:1294361/1‑59 (MQ=255)
                                                    cAATCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGcccc         <  1:255197/71‑1 (MQ=255)
                                                      aTCAGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCa                              <  1:2083583/48‑1 (MQ=255)
                                                         aGGATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGCCCCAGGt     <  1:1664755/70‑1 (MQ=255)
                                                           gATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGCCCCAGGTTAt  <  1:1310642/71‑1 (MQ=255)
                                                           gATCCGTTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGCCCCAGGTTAt  >  1:447222/1‑71 (MQ=255)
                                                                |                                                                 
GTAAAGAGCACCAGCAATGTAAAGACCTCGAGACGACCAAACAGCATGTTGGCAATCAGGATCCATTTAGCCACCGGGTTCATACTGGTAAAGTTATCAGCAACCACGCCAAGCCCTGGCCCCAGGTTAT  >  minE/2364039‑2364168

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: