Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,339,861 A→G 30.5% I272I (ATT→ATC wcaC ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,339,8610AG30.5% 77.1 / 40.6 85I272I (ATT→ATCwcaCpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (23/36);  new base G (9/17);  total (32/53)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

TTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAT  >  minE/1339795‑1339928
                                                                  |                                                                   
ttCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCa                                                                     >  1:326623/1‑67 (MQ=255)
ttCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAAt                                                                   >  1:1278472/1‑69 (MQ=255)
            aTTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCa                                                        >  1:1224751/1‑68 (MQ=255)
            aTTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAg                                                       >  1:1586361/1‑69 (MQ=255)
            aTTAACAACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAg                                                       >  1:443180/1‑69 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:2693189/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:2885540/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:2771869/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:1649917/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:1492647/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:1413812/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:2368874/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:811581/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:1194619/69‑1 (MQ=255)
                             cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa                                      <  1:232310/69‑1 (MQ=255)
                                 gAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAg                                  >  1:2710428/1‑69 (MQ=255)
                                 gAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAg                                  >  1:1514855/1‑69 (MQ=255)
                                 gAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAg                                  >  1:3328843/1‑69 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2966518/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2744914/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:3162967/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:1011143/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:3372693/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2600901/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2526104/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2409273/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:446969/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2342415/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:983302/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:215852/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2157958/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2125855/68‑1 (MQ=255)
                                    cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:1615416/68‑1 (MQ=255)
                                     ggCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt                                <  1:2865821/67‑1 (MQ=255)
                                      gCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACGAGtt                                <  1:2686486/66‑1 (MQ=255)
                                      gCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGt                                       >  1:2027649/1‑59 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:3065182/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:502891/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:1255864/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:566260/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:3000634/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:87173/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:968607/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt                         >  1:2090680/1‑69 (MQ=255)
                                          gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCAACAGTTGCTGCtt                         >  1:2306874/1‑69 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:3026953/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:3023344/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:2269770/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:360453/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:1488525/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:1610238/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:652391/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:1718869/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:1730500/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcgtcgt           <  1:1879665/69‑1 (MQ=255)
                                                        tATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCcggcgt           <  1:1315875/69‑1 (MQ=255)
                                                          tGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAg                                  >  1:728492/1‑44 (MQ=255)
                                                          tGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:118026/1‑69 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCTTCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:1410680/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCTTCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3116243/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:1494856/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:682053/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:729092/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:836992/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:2933027/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:1597128/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3093686/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3371830/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3358056/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3294567/1‑67 (MQ=255)
                                                            cAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTaa         >  1:3241930/1‑67 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:1680129/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:858127/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:237388/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:1221119/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:739413/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:559980/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:1671549/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:2773301/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:588905/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:2861006/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:3160841/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:2956296/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:1344771/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAt  <  1:3356250/69‑1 (MQ=255)
                                                                 cAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCCTAACGCAGAt  <  1:2308262/69‑1 (MQ=255)
                                                                  |                                                                   
TTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTTAGTTTTGCCGTCGTAACGCAGAT  >  minE/1339795‑1339928

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: