Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 579,722 C→T 22.3% Y37Y (TAC→TAT clpS → regulatory protein for ClpA substrate specificity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE579,7220CT22.3% 109.2 / 22.9 80Y37Y (TAC→TATclpSregulatory protein for ClpA substrate specificity
Reads supporting (aligned to +/- strand):  ref base C (37/25);  new base T (6/12);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGC  >  minE/579672‑579789
                                                  |                                                                   
gcgcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGtt                                                   <  1:3259888/69‑1 (MQ=255)
gcgcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGtt                                                   <  1:1908401/69‑1 (MQ=255)
     aaaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGa                                              >  1:3068745/1‑69 (MQ=255)
     aaaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGa                                              >  1:181632/1‑69 (MQ=255)
     aaaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGa                                              >  1:2190711/1‑69 (MQ=255)
     aaaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGa                                              >  1:1967968/1‑69 (MQ=255)
      aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:3126810/1‑69 (MQ=255)
      aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:3432397/1‑69 (MQ=255)
      aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:2603292/1‑69 (MQ=255)
      aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:1889819/1‑69 (MQ=255)
      aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:2021916/1‑69 (MQ=255)
      aaaCCGCAATCTATGTATAAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGAc                                             >  1:85192/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:1787133/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:33366/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:3285528/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:735140/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:1255136/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:2293609/1‑36 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTaca                                                                    >  1:2089646/1‑35 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATTGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:3280148/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTAtt                                                >  1:2049903/1‑55 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTAtt                                                >  1:3146429/1‑55 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTAtt                                                >  1:2051455/1‑55 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTAtt                                                >  1:1207400/1‑55 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTAtt                                                >  1:500267/1‑55 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:2472953/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:290173/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:2601490/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:3240550/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaaa                                  >  1:1177562/1‑69 (MQ=255)
                 tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACaaa                                   >  1:1180579/1‑68 (MQ=255)
                 tATGTAAAAAGTGATATTAGTCAATGATGATTacac                                                                   >  1:8000/1‑36 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2662567/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2506309/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:1686237/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:269190/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2693543/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2745278/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2772207/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:3035136/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2442705/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:1649272/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:1852120/69‑1 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACTAAAAttc                              >  1:918627/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:2813972/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:841318/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:2817833/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:3130668/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:3183289/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:548911/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:6248/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:652471/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:3283437/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:101615/1‑69 (MQ=255)
                     tataAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              >  1:2320053/1‑69 (MQ=255)
                      ataAAGTGATATTAGTCAATGATGATTATACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttc                              <  1:2896146/68‑1 (MQ=255)
                                 ttAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttct                        <  1:2459976/63‑1 (MQ=255)
                                 ttAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttct                        <  1:3147141/63‑1 (MQ=255)
                                 ttAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttct                        <  1:1388721/63‑1 (MQ=255)
                                 ttAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatga                   >  1:407401/1‑68 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:1049977/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:3370515/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:706028/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:836278/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:1474449/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:2866011/68‑1 (MQ=255)
                                   aGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgatg                 <  1:1214313/68‑1 (MQ=255)
                                       aaTGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGaa            >  1:2117367/1‑69 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1996999/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1480471/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:49105/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1387860/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1599118/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1768790/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:1287206/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:819344/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:2306790/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:2577558/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:3156980/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAg              <  1:2813746/63‑1 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGt         >  1:1302939/1‑68 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGt         >  1:1417128/1‑68 (MQ=255)
                                           atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGt         >  1:480277/1‑68 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:924720/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:1225106/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:1155130/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:1153455/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:1092465/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:607127/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:440475/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:407892/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:39546/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:3429057/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:1783904/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:206273/69‑1 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:3272500/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:2612486/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:3209266/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  >  1:2732798/1‑69 (MQ=255)
                                                 acacTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGc  <  1:2789882/69‑1 (MQ=255)
                                                  |                                                                   
GCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGC  >  minE/579672‑579789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: