Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,207 G→A 31.1% D141D (GAC→GAT yqeF ← predicted acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,2070GA31.1% 121.5 / 62.6 128D141D (GAC→GATyqeFpredicted acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (63/25);  new base A (19/21);  total (82/46)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 128 total aligned reads displayed.

ATTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGCGCG  >  minE/1932141‑1932272
                                                                  |                                                                 
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:2375423/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:1889772/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:3057567/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:3061197/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:1802992/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:3108552/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:1664612/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:314517/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:3187658/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:356680/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:142161/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:2626060/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:130647/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:2391514/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:478307/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:475833/69‑1 (MQ=255)
aTTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCACACAACCCATAATGCACAAGACTATc                                                                 <  1:3174929/69‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:2553251/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:2845843/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:668813/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:1150170/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:312881/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:2543224/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:1586788/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:3241209/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:1493596/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:2649524/68‑1 (MQ=255)
      ggcggTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTcaacca                                                            <  1:3367182/68‑1 (MQ=255)
        cggTTACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTATc                                                                 <  1:2689061/61‑1 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTgg                                                       >  1:1677136/1‑55 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGc                                                      >  1:2176401/1‑56 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGc                                                      >  1:2366528/1‑56 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGc                                                      >  1:395387/1‑56 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGc                                                      >  1:1298596/1‑56 (MQ=255)
                        aTAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGc                                                      >  1:1154466/1‑56 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCtgtg                                       >  1:1598044/1‑65 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      <  1:2881355/66‑1 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:2718852/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      <  1:2710370/66‑1 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      <  1:329537/66‑1 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:362515/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:66232/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:97078/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      <  1:2506034/66‑1 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:1835253/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:1879325/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:2232310/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      >  1:1928338/1‑66 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGc                                      <  1:1042489/66‑1 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCAAAg                                            >  1:2059450/1‑60 (MQ=255)
                              aTTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACAAACTGGCTATTGCCAAGCtg                                         >  1:1601934/1‑63 (MQ=255)
                                 gAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGt                                >  1:1465548/1‑69 (MQ=255)
                                 gAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGt                                >  1:2550454/1‑69 (MQ=255)
                                 gAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGt                                >  1:1132177/1‑69 (MQ=255)
                                 gAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGt                                >  1:2067484/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:415673/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:1377303/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:840871/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:420671/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:1163811/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:2466589/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:2405195/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:2260869/1‑69 (MQ=255)
                                          ccACAACCCATCATGCACAAGACTGTCAACAAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtc                       >  1:1885313/1‑69 (MQ=255)
                                          ccACAACCCATCATCCACAAGACTGTCAACCAACTGCCTAt                                                   >  1:1796323/1‑41 (MQ=38)
                                             caACCCCTCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:3999/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:2830759/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:1493016/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:1430663/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:3253715/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:2996779/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:48090/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:2228118/1‑67 (MQ=255)
                                             caACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtca                      >  1:1020548/1‑67 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCTACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:588843/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:1318673/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:620523/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:34464/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:3368303/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:1497417/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:1555699/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:3062949/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:889360/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2736835/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2070222/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2127979/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2607865/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2351238/1‑56 (MQ=255)
                                              aaCCCATCATGCACAAGACTATCAACCAACTGGCTATTGCCAAGCTGTGCGCCGGt                                >  1:2347218/1‑56 (MQ=255)
                                                      atGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtcagtcag                 >  1:1101517/1‑63 (MQ=255)
                                                      atGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTAtcagtcag                 >  1:1661313/1‑63 (MQ=255)
                                                        gCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACAtgtg         >  1:2132774/1‑69 (MQ=255)
                                                         cacaAGACTGTCAACCAACTGTCTATTGCCGAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATgtgg        <  1:3049532/69‑1 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTgcgc   >  1:3388482/1‑68 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGcgcg  >  1:1039429/1‑69 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGcgcg  >  1:3253581/1‑69 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGcgcg  >  1:2975928/1‑69 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGcgcg  >  1:2888652/1‑69 (MQ=255)
                                                               aCTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGcgcg  >  1:2170787/1‑69 (MQ=255)
                                                                cTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGGGCGGCTATCAGTCAGAACATgtgg        <  1:1747730/62‑1 (MQ=255)
                                                                  |                                                                 
ATTTTCGGCGGTGACACCAATATGATAATCATTGAAGGCATCCCACAACCCATCATGCACAAGACTGTCAACCAACTGGCTATTGCCAAGCTGTGCACCGGTGCGGCTATCAGTCAGAACATGTGGTGCGCG  >  minE/1932141‑1932272

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: