Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,425 C→G 24.7% N545K (AAC→AAG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,4250CG24.7% 166.9 / 51.4 89N545K (AAC→AAGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (37/30);  new base G (18/4);  total (55/34)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTG  >  minE/386356‑386494
                                                                     |                                                                     
aGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAaca                                                                      >  1:2280740/1‑71 (MQ=255)
aGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAaca                                                                      >  1:2785264/1‑71 (MQ=255)
aGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAaca                                                                      >  1:3332413/1‑71 (MQ=255)
aGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAaca                                                                      >  1:15669/1‑71 (MQ=255)
  cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    <  1:5171616/71‑1 (MQ=255)
  cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    <  1:3681459/71‑1 (MQ=255)
  cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    <  1:4234938/71‑1 (MQ=255)
  cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    >  1:4569652/1‑71 (MQ=255)
   gcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAaca                                                                      <  1:2365723/68‑1 (MQ=255)
    cTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActct                                                                  >  1:1258156/1‑71 (MQ=255)
    cTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActct                                                                  >  1:1434803/1‑71 (MQ=255)
     tgtgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    >  1:4495123/1‑68 (MQ=255)
     tgtgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                    >  1:2042437/1‑68 (MQ=255)
       tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                                                 <  1:4539334/69‑1 (MQ=255)
       tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGACCCGCTGGGCGGCATTAACACTCtg                                                                 <  1:5365051/69‑1 (MQ=255)
           tACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCg                                                           <  1:4755316/71‑1 (MQ=255)
             cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtgtg                                                               >  1:2113715/1‑65 (MQ=255)
               ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                                       >  1:5405944/1‑71 (MQ=255)
                ctgTGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGtt                                                      >  1:4246239/3‑71 (MQ=255)
                 cGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                     <  1:1226240/71‑1 (MQ=255)
                   tGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                   <  1:461061/71‑1 (MQ=255)
                   tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTgg                                                              >  1:3887993/1‑60 (MQ=255)
                   tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTgg                                                              >  1:1617627/1‑60 (MQ=255)
                   tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGttcg                                                    <  1:4717093/70‑3 (MQ=255)
                   tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCgg                                                   >  1:631643/1‑71 (MQ=255)
                      ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                <  1:143491/71‑1 (MQ=255)
                      ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                >  1:3791246/1‑71 (MQ=255)
                          aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActc                                                                   <  1:2358133/48‑1 (MQ=255)
                           cTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                                                 >  1:2041384/1‑49 (MQ=255)
                           cTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCt                                                         <  1:3705264/57‑1 (MQ=255)
                            ttCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                                       <  1:4154420/58‑1 (MQ=255)
                            ttCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCaa                                          >  1:3665008/1‑71 (MQ=255)
                             tcctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                                         >  1:568502/1‑71 (MQ=255)
                              cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                     >  1:1398909/1‑58 (MQ=255)
                              cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                 >  1:405259/1‑62 (MQ=255)
                                tgtaTCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTAt                                                >  1:469441/4‑61 (MQ=255)
                                 ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                 >  1:5199583/1‑59 (MQ=255)
                                 ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                >  1:2478397/1‑60 (MQ=255)
                                 ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                >  1:4503111/1‑60 (MQ=255)
                                  taTCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                     <  1:4894737/69‑1 (MQ=255)
                                  cATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAtt                                               <  1:675295/60‑1 (MQ=255)
                                   aTCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGAt                                    >  1:3166269/1‑70 (MQ=255)
                                   aTCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATg                                   >  1:5089241/1‑71 (MQ=255)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                         >  1:2601680/1‑48 (MQ=21)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                         >  1:5307376/1‑48 (MQ=21)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                         >  1:4005202/1‑48 (MQ=21)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCgg                                                   >  1:21154/1‑54 (MQ=25)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGc                                  >  1:512031/1‑71 (MQ=255)
                                    tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGc                                  <  1:3424525/71‑1 (MQ=255)
                                    tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                   <  1:4409241/54‑1 (MQ=255)
                                    tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                   <  1:5247843/54‑1 (MQ=255)
                                    tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                                         <  1:2615699/64‑1 (MQ=255)
                                    tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAg                                      <  1:3803587/67‑1 (MQ=255)
                                    tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGc                                  >  1:3092143/1‑71 (MQ=255)
                                      aGGTCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTg                                >  1:4025979/1‑71 (MQ=19)
                                      aGGGCGTGGTCGATCCGCTGGGCGGGATCAACACCCTGTGGCCGCTGTTCGGCATCGCCAACCAGATGCTg                                >  1:598713/1‑71 (MQ=255)
                                      aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                 >  1:4872730/1‑54 (MQ=255)
                                      aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAtt                                               >  1:2362930/1‑56 (MQ=255)
                                      aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTg                                >  1:764011/1‑71 (MQ=255)
                                      aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTg                                >  1:892867/1‑71 (MQ=255)
                                          cGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTGGCAg                            >  1:5376531/1‑71 (MQ=25)
                                          cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                >  1:3416126/1‑35 (MQ=255)
                                          cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                >  1:315017/1‑35 (MQ=255)
                                          cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                >  1:2767976/1‑35 (MQ=255)
                                          cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                >  1:2761479/1‑35 (MQ=255)
                                           gTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAgg                           >  1:5403478/1‑71 (MQ=255)
                                           gTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAgg                           <  1:705262/71‑1 (MQ=255)
                                              gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                     >  1:5442634/1‑42 (MQ=255)
                                              gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                     >  1:793803/1‑42 (MQ=255)
                                              gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                     >  1:1523083/1‑42 (MQ=255)
                                              gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGt                                                  >  1:4778075/1‑45 (MQ=255)
                                               ttgATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAAc                                         >  1:943000/3‑53 (MQ=16)
                                               ttgATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                               <  1:3797837/61‑1 (MQ=19)
                                                cGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTg                                                    >  1:4795395/1‑41 (MQ=255)
                                                 gATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGt                                                  >  1:3596631/1‑42 (MQ=19)
                                                 gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                <  1:1805577/44‑1 (MQ=255)
                                                 gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGgc                     <  1:2984461/71‑1 (MQ=255)
                                                     cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                  <  1:5492450/70‑1 (MQ=255)
                                                     cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                  <  1:5444251/70‑1 (MQ=255)
                                                      gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                     >  1:146458/1‑50 (MQ=255)
                                                      gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                     >  1:2204434/1‑50 (MQ=255)
                                                      gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                               >  1:4957401/1‑56 (MQ=16)
                                                      gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                               >  1:1378796/1‑56 (MQ=16)
                                                      gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                        <  1:2398335/63‑1 (MQ=255)
                                                      gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                  <  1:146810/69‑1 (MQ=255)
                                                      gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                  <  1:2758244/69‑1 (MQ=255)
                                                       cTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                        <  1:490666/62‑1 (MQ=255)
                                                          ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATg              <  1:4941133/69‑1 (MQ=255)
                                                          ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGc             <  1:1693235/70‑1 (MQ=255)
                                                                 ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTgg                               <  1:3298866/45‑1 (MQ=255)
                                                                 ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGcc     <  1:4829369/71‑1 (MQ=255)
                                                                  aaaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATg                                   >  1:4734756/2‑40 (MQ=255)
                                                                   aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt   >  1:1956186/1‑71 (MQ=255)
                                                                   aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt   >  1:3141143/1‑71 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGc                              >  1:1600008/1‑43 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                             >  1:976486/1‑44 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                             >  1:891339/1‑44 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                             >  1:2459854/1‑44 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCtgt        >  1:5453234/1‑65 (MQ=255)
                                                                    acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtg  >  1:22249/1‑71 (MQ=255)
                                                                     |                                                                     
AGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTG  >  minE/386356‑386494

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: