Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,236 A→G 17.4% E207E (GAA→GAG yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,449,2360AG17.4% 115.0 / 13.8 74E207E (GAA→GAGyifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (22/39);  new base G (11/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

AACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCGGG  >  minE/2449175‑2449304
                                                             |                                                                    
aaCCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGaa                                                               <  1:1489544/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:2357125/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:1832340/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:1693596/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:183824/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:1634217/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:1378014/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:2472533/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:2779690/69‑1 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         >  1:1281545/1‑69 (MQ=255)
      cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                                         <  1:2913066/69‑1 (MQ=255)
       taaTGGAGTATCCGGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:1219441/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:240872/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:1982254/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:2314702/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:365019/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:386263/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:1939207/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:2899712/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:2941252/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:610277/3‑69 (MQ=255)
       taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                                        >  1:1454070/3‑69 (MQ=255)
              gTATCCCTGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:2986567/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:1669526/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:1190009/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:2965946/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:3125507/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:382884/68‑1 (MQ=255)
              gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                                  <  1:1628894/68‑1 (MQ=255)
                aTCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGaa                                                               >  1:2392245/1‑53 (MQ=255)
                 tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                              >  1:1900569/1‑69 (MQ=255)
                 tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                              >  1:453375/1‑69 (MQ=255)
                 tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                              >  1:402225/1‑69 (MQ=255)
                   ccGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTAc                                            >  1:2744156/1‑69 (MQ=255)
                   ccGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTAc                                            >  1:835885/1‑69 (MQ=255)
                      ggCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCgc                                         >  1:312495/1‑69 (MQ=255)
                      ggCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCgc                                         >  1:247327/1‑69 (MQ=255)
                                   gATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAGGACTGGAAATTACCGCCGCTGGCGGGcac                            <  1:1823179/69‑1 (MQ=255)
                                   gATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAGGACTGGAAATTACCGCCGCTGGCGGGcac                            <  1:3403622/69‑1 (MQ=255)
                                    aTCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                           >  1:1295169/1‑69 (MQ=255)
                                    aTCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                           <  1:184523/69‑1 (MQ=255)
                                    aTCTCAGTGATGTTATCGGGCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                           <  1:2927343/69‑1 (MQ=255)
                                         aGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                           <  1:621709/64‑1 (MQ=255)
                                           tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:858766/68‑1 (MQ=255)
                                           tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:2761995/68‑1 (MQ=255)
                                           tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:2694696/68‑1 (MQ=255)
                                            gATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:1622984/67‑1 (MQ=255)
                                             aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa                   >  1:2959060/1‑68 (MQ=255)
                                             aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa                   >  1:2655532/1‑68 (MQ=255)
                                             aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa                   >  1:3147232/1‑68 (MQ=255)
                                             aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa                   >  1:3252177/1‑68 (MQ=255)
                                             aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa                   >  1:2379630/1‑68 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGaaa                                                >  1:294686/1‑37 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:27551/64‑1 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:3030291/64‑1 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:1012154/64‑1 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:1636780/64‑1 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:3283280/64‑1 (MQ=255)
                                               gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:397725/64‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:3161025/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:2185835/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:359806/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:2124729/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:3248986/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:1614186/63‑1 (MQ=255)
                                                ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt                     <  1:2751905/63‑1 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:2371566/1‑69 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:214785/1‑69 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:1318171/1‑69 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:827088/1‑69 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:2520094/1‑69 (MQ=255)
                                                  aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt             >  1:1035586/1‑69 (MQ=255)
                                                           gAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGgccgcc     <  1:2996963/68‑1 (MQ=255)
                                                           gAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGgccgcc     <  1:1331074/68‑1 (MQ=255)
                                                             aCAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCggg  <  1:2198720/69‑1 (MQ=255)
                                                             aCAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCggg  <  1:1106718/69‑1 (MQ=255)
                                                             aCAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCggg  <  1:3429419/69‑1 (MQ=255)
                                                             aCAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCggg  <  1:3236963/69‑1 (MQ=255)
                                                             aCAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCggg  <  1:2893902/69‑1 (MQ=255)
                                                             |                                                                    
AACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATTGGGCCGCCGGG  >  minE/2449175‑2449304

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: