Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,066,835 C→G 42.2% V219L (GTT→CTT)  ydhB ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,066,8350CG42.2% 36.7 / 84.4 83V219L (GTT→CTT) ydhBpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (42/6);  new base G (21/14);  total (63/20)
Fisher's exact test for biased strand distribution p-value = 8.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ATCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCGCGC  >  minE/1066787‑1066899
                                                |                                                                
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTAt                                                    >  1:242289/1‑63 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:1161682/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:978293/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:731875/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:343399/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:3303098/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:2970478/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:2384896/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:2326170/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:23041/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:2133012/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:1933542/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:1895111/1‑69 (MQ=255)
aTCCTGCCGAGATACAGGTTACCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCa                                              >  1:433271/1‑69 (MQ=255)
      ccGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:3021249/69‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:1082630/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:2998814/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:1265069/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:124510/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:152107/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:607244/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:642831/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:691434/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:2069861/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:803423/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:1906275/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:1927246/67‑1 (MQ=255)
        gagaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAgg                                        <  1:87336/67‑1 (MQ=255)
         agaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTaa                                     <  1:1979081/69‑1 (MQ=255)
         agaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTaa                                     <  1:2418623/69‑1 (MQ=255)
         agaTACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTaa                                     <  1:3069031/69‑1 (MQ=255)
              cAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACg                                 >  1:3346171/1‑68 (MQ=255)
                       gatgatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAAc                                            <  1:9658/48‑1 (MQ=255)
                        atgatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCg                      <  1:2198207/69‑1 (MQ=255)
                        atgatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCg                      <  1:1397697/69‑1 (MQ=255)
                          gatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTa                                                     >  1:177743/1‑36 (MQ=255)
                          gatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTa                                                     >  1:849466/1‑36 (MQ=255)
                          gatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTa                                                     >  1:3408426/1‑36 (MQ=255)
                          gatTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGtt                    >  1:2350869/1‑69 (MQ=255)
                          gatTACCAGTCAGGCACGACAACTCTTTTTTGGtt                                                      >  1:2468504/1‑35 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:1601676/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:954551/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:137553/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:943445/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:419793/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:449595/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:586422/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:808310/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:804311/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2963494/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2905055/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2822710/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:1592164/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:1560874/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:1896302/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2324995/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:14817/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2338369/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcgc                 >  1:2348122/1‑67 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTcg                  >  1:2936496/1‑66 (MQ=255)
                               ccAGTCAGGCACGACAAGTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGAGTTcgc                 >  1:2446048/1‑67 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:161350/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:3279803/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:943532/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:564726/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:1752870/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:3119311/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:2162053/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:738034/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTCAGGTAGCGTTCGCGAGgtg           >  1:2587455/1‑69 (MQ=255)
                                   tCAGGCACGACAACTCTTTTTTGGTTAGCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGgtg           >  1:1759356/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1184706/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1551653/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1553907/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:3131646/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:2222861/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:810312/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:2248189/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1675509/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:464417/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1641772/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:1452828/1‑69 (MQ=255)
                                            aCAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCgcgc  >  1:2368649/1‑69 (MQ=255)
                                                |                                                                
ATCCTGCCGAGATACAGGTTGCCGATGATTCCCAGTCAGGCACGACAACTCTTTTTTGGTTATCCAGCAACCAGGTAATACGTTTAGGTAGCGTTCGCGAGGTGTCTTCGCGC  >  minE/1066787‑1066899

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