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,2720AG99.3% 561.2 / ‑5.0 146L50L (CTT→CTCuspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (77/68);  minor base T (0/1);  total (77/69)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 146 total aligned reads displayed.

AGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATAAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAGGG  >  minE/2605213‑2605337
                                                           |                                                                 
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:2739716/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:941146/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:578813/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:2912115/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:253134/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:2353398/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:1903498/71‑1 (MQ=255)
aGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:1058313/71‑1 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:755186/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:2456797/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:1315726/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:543984/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:2676540/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAgg                                                    >  1:2789532/1‑71 (MQ=255)
    cATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     >  1:2255871/1‑70 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATTAGCCCGGTGTATAg                                                     <  1:3727199/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:3712142/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:3637697/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:3090853/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:4159927/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:1021927/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:4294535/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:691443/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:2413174/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:2364403/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:2809906/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:1711013/68‑1 (MQ=255)
      tgatgTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGCGTATAg                                                     <  1:564740/68‑1 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:3171388/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:3381676/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:3676676/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:1498482/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:2956365/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:2187920/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:74336/1‑71 (MQ=255)
             tctcTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGagtag                                           >  1:1280671/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGc                                                                >  1:2622872/1‑44 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGgtgt                                                         >  1:2164489/1‑51 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGTTGTATAgg                                                    >  1:1896683/1‑56 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTAc                                      <  1:4379386/70‑1 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTAc                                      <  1:288438/70‑1 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTAc                                      <  1:812464/70‑1 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTAc                                      >  1:3792250/1‑70 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                                     >  1:2856438/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                                     >  1:1241564/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                                     >  1:3797006/1‑71 (MQ=255)
                   cAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACa                                     >  1:2217748/1‑71 (MQ=255)
                               tatGCATATCACCCAGATTCACATCAATGAGCCCGGTGta                                                        <  1:4106166/38‑1 (MQ=38)
                               tCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGAt                          >  1:328505/1‑70 (MQ=255)
                               tCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGAt                          >  1:1707049/1‑70 (MQ=255)
                               tCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGAt                          >  1:1550258/1‑70 (MQ=255)
                               tCTGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGAt                          >  1:1263790/1‑70 (MQ=255)
                                 tGCATATCACCCAGATTCACATCAATGAGCCCGGTGTATAg                                                     <  1:1824814/41‑1 (MQ=255)
                                      atCACCCAGATTCGCATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:4034563/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTCCGTGGATCAGAGaaa                  >  1:3390960/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCaga                      >  1:1193372/1‑67 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:3141374/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:2334963/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:1630126/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:3250279/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:4128579/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:1323661/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:3054941/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:1094133/1‑71 (MQ=255)
                                      atCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGaaa                  >  1:885517/1‑71 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAAc                 >  1:3939258/1‑71 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAAc                 <  1:2051726/71‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAAc                 <  1:1986815/71‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAAc                 <  1:1540750/71‑1 (MQ=255)
                                       tCACCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAAc                 <  1:1233095/71‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:1283183/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:1337436/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:3527256/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:2774295/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:1911878/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:1593164/70‑1 (MQ=255)
                                         aCCCAGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACt                <  1:1426197/70‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATATGTCAGAGTAGTTTACATCTACGTGGAGCAGAGAAACTTTCGc           <  1:4029030/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:67694/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:3558464/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:936707/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:11721/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:3222969/71‑1 (MQ=255)
                                             aGATTCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGc           <  1:1325319/71‑1 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:1950393/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:208521/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:828065/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:133031/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:222070/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:1413329/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAtt        >  1:1226775/1‑71 (MQ=255)
                                                ttCACATCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCAt         >  1:1597258/1‑70 (MQ=255)
                                                   acaTCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGt      <  1:4364258/70‑1 (MQ=255)
                                                   acaTCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGt      <  1:2229364/70‑1 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTAc                                >  1:4079142/1‑41 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAggg  >  1:249057/1‑71 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAggg  >  1:4260980/1‑71 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAggg  >  1:852872/1‑71 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAggg  >  1:2813100/1‑71 (MQ=255)
                                                      tCAATGAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAgg   >  1:1587958/1‑70 (MQ=255)
                                                           |                                                                 
AGTGCATGATGTGTCTCTTCAGAGATGCGTTTCTGCATATCACCCAGATTCACATCAATAAGCCCGGTGTATAGGTCAGAGTAGTTTACATCTACGTGGATCAGAGAAACTTTCGCATTGTAGGG  >  minE/2605213‑2605337

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: