Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,605,272 A→G L50L (CTT→CTC)  uspA ← universal stress global response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,605,2720AG100.0% 483.1 / NA 127L50L (CTT→CTC) uspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (58/69);  total (58/69)

Only 100 of 127 total aligned reads displayed.

AGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATAAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAA  >  minE/2605213‑2605321
                                                           |                                                 
aGTGCATGATGTGTCTCTTCAGAGTTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:1388361/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:1476304/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:136717/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:181584/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:1835467/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:3096182/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:2672930/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:2795230/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:2669075/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:381854/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:1143351/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:2461782/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:627068/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:2407611/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGccc                                              >  1:794199/1‑65 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGcac                                              >  1:362523/1‑63 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2920764/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:1419113/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2848855/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:1739231/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:3082238/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:1206327/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2597066/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2422572/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2072764/70‑1 (MQ=255)
           tgtCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:1677955/70‑1 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTAc                      >  1:2383036/1‑70 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1801660/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1900068/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:2463235/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1713728/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1952762/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1978410/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:2231484/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:166/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1654568/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1622526/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:682777/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:580855/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:984694/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:3275572/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1281817/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1360516/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:3013571/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:908605/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:28115/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1528999/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:155356/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1582940/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                     >  1:1650123/1‑71 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACAt                    >  1:3289325/1‑69 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:45576/1‑70 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:656838/1‑70 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:76463/1‑70 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:1011363/1‑70 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:10381/1‑70 (MQ=255)
                      agaTGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATc                   >  1:1338646/1‑70 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1178558/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1456522/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:910986/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:3016812/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1364353/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:3140000/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:3254915/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1254774/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1231797/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:231579/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1164770/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:113290/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:563895/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:572569/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1106111/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:109439/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1029603/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2890251/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2082357/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:842884/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2057411/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:903140/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2025319/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1932972/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2480523/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1785116/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2565166/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2578606/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1579109/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1578197/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1521186/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1511902/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1470952/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2898109/71‑1 (MQ=255)
                         tGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGGCAGAGTAGTTTACATCTACg               <  1:1466274/71‑1 (MQ=255)
                          gCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1025465/70‑1 (MQ=255)
                           cGTTTCTGCATATCACCCAGATTCACATCTATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2077147/69‑1 (MQ=255)
                                     tatCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:2561155/59‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa  <  1:871651/70‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa  <  1:1418962/70‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa  <  1:1415277/70‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCagag                              <  1:2051066/36‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACg               <  1:1179278/51‑1 (MQ=255)
                                                           |                                                 
AGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATAAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAA  >  minE/2605213‑2605321

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: ‑