Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,318,210 G→A 33.0% L135F (CTC→TTC) ‡ yiiT ← stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,318,2100GA33.0% 108.0 / 67.4 127L135F (CTC→TTC) ‡yiiTstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (35/50);  new base A (24/18);  total (59/68)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 127 total aligned reads displayed.

TGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGT  >  minE/2318157‑2318275
                                                     |                                                                 
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:1938670/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:38755/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:1722458/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:512373/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:2474422/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:2582260/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:1889580/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:3376532/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:1968154/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:2812690/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:1993602/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                                    >  1:2188824/1‑68 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACt                                                       >  1:331857/1‑66 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:3248113/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:776519/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:558360/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:1625997/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:2392176/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:212330/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                                    >  1:2325765/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATACGCACTCAt                                                    >  1:3366452/1‑69 (MQ=255)
tGACAGCGGAGCATTACTGGTAAGCGGTGACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCa                                                     >  1:1447297/1‑68 (MQ=255)
  aCAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCt                                                  >  1:1466921/1‑69 (MQ=255)
    aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                                                <  1:3157704/69‑1 (MQ=255)
    aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                                                <  1:702766/69‑1 (MQ=255)
    aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                                                <  1:1432227/69‑1 (MQ=255)
        gAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGa                                            <  1:203813/69‑1 (MQ=255)
         aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                                           <  1:1885328/69‑1 (MQ=255)
         aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                                           <  1:2991962/69‑1 (MQ=255)
         aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                                           <  1:979244/69‑1 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCaa                                                                >  1:958562/1‑42 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGTTCATccc                                     >  1:1185307/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:1818323/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:2341770/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:1045051/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:1135818/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:2012536/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:563942/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:259645/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:2436504/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                                     >  1:1708931/1‑69 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATcc                                      >  1:991363/1‑68 (MQ=255)
               aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATcc                                      >  1:1460590/1‑68 (MQ=255)
                  ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                                   <  1:3025806/68‑1 (MQ=255)
                  ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                                   <  1:1298341/68‑1 (MQ=255)
                  ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                                   <  1:893472/68‑1 (MQ=255)
                  ggTAAGCGGATACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                                   <  1:99010/68‑1 (MQ=255)
                   gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCcgcg                                  >  1:3099240/1‑68 (MQ=255)
                   gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGa                                 >  1:1813008/1‑69 (MQ=255)
                   gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGa                                 >  1:3360639/1‑69 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:1373795/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:2820287/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:1281162/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:1311216/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:512472/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:71199/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:643074/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:591391/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:568853/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:2052908/67‑1 (MQ=255)
                     aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:2379561/67‑1 (MQ=255)
                      aGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:2425916/66‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:384003/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:459810/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:3233883/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:1029141/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:3128580/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:1873406/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:1423617/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:204170/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:1418964/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:1310622/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:2622665/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:2696547/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:2732019/69‑1 (MQ=255)
                          gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:279447/69‑1 (MQ=255)
                           ttACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg                          <  1:2251773/68‑1 (MQ=255)
                             aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc                        >  1:278587/1‑68 (MQ=255)
                             aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc                        >  1:2529629/1‑68 (MQ=255)
                             aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc                        >  1:726889/1‑68 (MQ=255)
                             aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc                        >  1:1705759/1‑68 (MQ=255)
                                tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:1235985/56‑1 (MQ=255)
                                tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:1084975/56‑1 (MQ=255)
                                tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:3059815/56‑1 (MQ=255)
                                tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                                 <  1:580748/56‑1 (MQ=255)
                                 aTCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCaa                   >  1:622222/1‑69 (MQ=255)
                                 aTCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCaa                   >  1:1799699/1‑69 (MQ=255)
                                 aTCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCaa                   >  1:2561588/1‑69 (MQ=255)
                                          cGGCACGATGAGCAAATCCGCACTCATCTTATTGat                                           <  1:2477096/36‑1 (MQ=255)
                                             cACGATGAGCAAATC‑‑CACTCATCTTATTGATCATCCCGCGa                                 <  1:87237/41‑1 (MQ=37)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:583092/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:585790/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:1863545/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:1992895/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:2547539/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:802510/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:1290304/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:2585740/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:3194878/69‑1 (MQ=255)
                                                  tGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGt  <  1:1004454/69‑1 (MQ=255)
                                                     |                                                                 
TGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAAACGGTTGATAAATGAGT  >  minE/2318157‑2318275

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: