Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,225 A→G 19.6% I204V (ATC→GTC)  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,2250AG19.6% 94.7 / 15.0 66I204V (ATC→GTC) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (16/37);  new base G (11/2);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 9.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATT  >  minE/2449157‑2449293
                                                                    |                                                                    
aCGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTa                                                                      <  1:1982651/69‑1 (MQ=255)
aCGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTa                                                                      <  1:1566981/69‑1 (MQ=255)
                  aaCCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGaa                                                    <  1:1489544/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:1832340/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:1693596/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:2357125/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:1634217/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:2472533/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:2779690/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:1378014/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:2913066/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              <  1:183824/69‑1 (MQ=255)
                        cTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAg                                              >  1:1281545/1‑69 (MQ=255)
                         taaTGGAGTATCCGGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:1219441/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:1939207/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:1982254/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:610277/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:2314702/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:240872/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:2899712/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:2941252/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:365019/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:386263/3‑69 (MQ=255)
                         taaTGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAgg                                             >  1:1454070/3‑69 (MQ=255)
                                gTATCCCTGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:2986567/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:2965946/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:1628894/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:382884/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:1669526/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:1190009/68‑1 (MQ=255)
                                gTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGa                                       <  1:3125507/68‑1 (MQ=255)
                                  aTCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGaa                                                    >  1:2392245/1‑53 (MQ=255)
                                   tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                   >  1:402225/1‑69 (MQ=255)
                                   tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                   >  1:453375/1‑69 (MQ=255)
                                   tCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAAtt                                   >  1:1900569/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:3403622/69‑1 (MQ=255)
                                                     gATCTCAGTGATGTTGTCGGTCAGGAGCAAGGAAAGCGAGGACTGGAAATTACCGCCGCTGGCGGGcac                 <  1:1823179/69‑1 (MQ=255)
                                                      aTCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                <  1:184523/69‑1 (MQ=255)
                                                      aTCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                >  1:1295169/1‑69 (MQ=255)
                                                      aTCTCAGTGATGTTATCGGGCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                <  1:2927343/69‑1 (MQ=255)
                                                           aGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGcaca                <  1:621709/64‑1 (MQ=255)
                                                             tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:2761995/68‑1 (MQ=255)
                                                             tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:2694696/68‑1 (MQ=255)
                                                             tGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:858766/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:3252177/1‑68 (MQ=255)
                                                               aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa        >  1:2379630/1‑68 (MQ=255)
                                                               aTGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTa        >  1:3147232/1‑68 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGaaa                                     >  1:294686/1‑37 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:1636780/64‑1 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:3283280/64‑1 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:397725/64‑1 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:3030291/64‑1 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:1012154/64‑1 (MQ=255)
                                                                 gTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:27551/64‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:2124729/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:2185835/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:1614186/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:359806/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:3248986/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:3161025/63‑1 (MQ=255)
                                                                  ttATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCttt          <  1:2751905/63‑1 (MQ=255)
                                                                    aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt  >  1:827088/1‑69 (MQ=255)
                                                                    aTCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGAtt  >  1:2520094/1‑69 (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:1035586/1‑69 (MQ=255)
                                                                    |                                                                    
ACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGCTACAACATGATCTCAGTGATGTTATCGGTCAGGAACAAGGAAAGCGAGGACTGGAAATTACCGCTGCTGGCGGGCACAACCTTTTACTGATT  >  minE/2449157‑2449293

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: