Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,088,530 T→C 34.7% N296S (AAC→AGC)  yhaJ ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,088,5300TC34.7% 29.3 / 23.8 46N296S (AAC→AGC) yhaJpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base C (5/11);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

TAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAG  >  minE/2088467‑2088595
                                                               |                                                                 
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:217753/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:540606/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:2029723/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:503209/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:1603090/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:1755699/1‑69 (MQ=255)
tAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTaaaa                                                              >  1:1732910/1‑69 (MQ=255)
           ttGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGaa                                                    >  1:3188100/1‑68 (MQ=255)
           ttGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAAt                                                   >  1:1416312/1‑69 (MQ=255)
           ttGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAAt                                                   >  1:1835889/1‑69 (MQ=255)
           ttGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAAt                                                   >  1:1726482/1‑69 (MQ=255)
                        tCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGacac                                      <  1:667540/69‑1 (MQ=255)
                        tCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGacac                                      <  1:1717713/69‑1 (MQ=255)
                            cccTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAg                                   <  1:462871/68‑1 (MQ=255)
                            cccTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAg                                   <  1:31946/68‑1 (MQ=255)
                            cccTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAg                                   <  1:2444937/68‑1 (MQ=255)
                            cccTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAg                                   <  1:2613862/68‑1 (MQ=255)
                                       aaCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTc                        <  1:1406345/68‑1 (MQ=255)
                                       aaCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTc                        <  1:2240388/68‑1 (MQ=255)
                                           ataataAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCa                   <  1:2677064/69‑1 (MQ=255)
                                           ataataAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTCCCCCCa                   <  1:3285322/69‑1 (MQ=255)
                                            taataaATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:2204017/69‑1 (MQ=255)
                                                 aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:469123/64‑1 (MQ=255)
                                                 aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:3260228/64‑1 (MQ=255)
                                                 aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:549974/64‑1 (MQ=255)
                                                 aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:235389/64‑1 (MQ=255)
                                                 aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:1947124/64‑1 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:2757861/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:1382427/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:503778/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:1484292/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:1485221/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:1699453/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:454171/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:2134144/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:3088635/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGcc                          <  1:30429/52‑2 (MQ=255)
                                                     ttattTTCCGCTAAAAAGTTTGGGAATTTCCCCAGACACCAGGATTTCGcc                          <  1:1184250/51‑2 (MQ=255)
                                                      tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTcc                       >  1:658538/1‑54 (MQ=255)
                                                      tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga        >  1:2307818/1‑69 (MQ=255)
                                                      tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga        >  1:2882405/1‑69 (MQ=255)
                                                      tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga        >  1:2811477/1‑69 (MQ=255)
                                                      tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga        >  1:1046123/1‑69 (MQ=255)
                                                          ttCCGTTAAAAAGTTTGGGAATTTCCCGCAGCCACCAGGATTTCGCTTCCCCCATACTGTCACGACGcc    >  1:2954089/1‑69 (MQ=255)
                                                          ttCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGcc    >  1:2221668/1‑69 (MQ=255)
                                                             cGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAg  <  1:3155856/68‑1 (MQ=255)
                                                               ttAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                  <  1:128825/50‑1 (MQ=255)
                                                               |                                                                 
TAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAG  >  minE/2088467‑2088595

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: