Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,582,830 G→C 40.9% T135T (ACC→ACG ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,582,8300GC40.9% 31.8 / 64.5 66T135T (ACC→ACGligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base G (27/12);  new base C (15/12);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

ATGGCGCGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGC  >  minE/1582771‑1582890
                                                           |                                                            
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:66529/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:64265/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:638846/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:2176737/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:216592/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:2152940/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:1980104/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:3426917/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:1912434/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:3151415/69‑1 (MQ=255)
aTGGCGCGAATAGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:1672/69‑1 (MQ=255)
   gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTg                                                  <  1:2351466/69‑1 (MQ=255)
   gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTg                                                  <  1:1784621/69‑1 (MQ=255)
   gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTg                                                  <  1:1501187/69‑1 (MQ=255)
   gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTg                                                  <  1:1254119/69‑1 (MQ=255)
   gcgcGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGACATCGCCACGGGTCGCGGCACTg                                                  <  1:464270/69‑1 (MQ=255)
           aGTACGCACATTAGACGTAATATCTTCCCCGGTGGTGCCATCGCCACGCGTCGCGGCa                                                     <  1:2837770/58‑1 (MQ=255)
                    aTTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTg                                                  <  1:1181614/52‑1 (MQ=255)
                          cGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:245243/69‑1 (MQ=255)
                          cGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:2920908/69‑1 (MQ=255)
                              atatCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATa                       >  1:2666041/1‑69 (MQ=255)
                              atatCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATa                       >  1:1431924/1‑69 (MQ=255)
                              atatCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATa                       >  1:2321115/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTGAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2052309/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTg                    >  1:2505117/1‑57 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCaa             >  1:3361556/1‑64 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2428719/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:1025868/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2745757/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3065531/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3194908/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3201362/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3286332/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3430094/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2427106/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:1724333/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:100997/1‑69 (MQ=255)
                                             gtgCCATCGCCACGCGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2222910/1‑69 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATGCAGAATACTGACGGCAAGACCa        >  1:1523961/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:1599731/48‑1 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           >  1:1772664/1‑48 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           >  1:3361236/1‑48 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:1383878/48‑1 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:1327514/48‑1 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           <  1:682526/48‑1 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAg                           >  1:2377276/1‑48 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:1174080/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3326762/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3314203/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:851031/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:963683/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:1654210/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:3160478/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:1990330/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2831245/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2818043/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2568902/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2491374/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCa        >  1:2398192/1‑67 (MQ=255)
                                               gCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGAACa        >  1:1079094/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:1751583/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:1304398/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:2737679/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:107048/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:2340336/1‑67 (MQ=255)
                                                     gCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGc  >  1:674260/1‑67 (MQ=255)
                                                           |                                                            
ATGGCGCGAATAGTACGCACATTAGACGTGATATCTTCCCCGGTGGTGCCATCGCCACGGGTCGCGGCACTGACTAAAACGCCATTTTCATACAGAATACTGACGGCAAGACCATCCAGC  >  minE/1582771‑1582890

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: