Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,339,849 G→A 36.5% T276T (ACC→ACT wcaC ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,339,8490GA36.5% 44.6 / 41.6 71T276T (ACC→ACTwcaCpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (9/36);  new base A (9/17);  total (18/53)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TTTGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTT  >  minE/1339786‑1339905
                                                               |                                                        
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1926941/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2616364/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2690455/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2588202/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2403325/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2382556/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2820324/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2927724/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2274767/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:3359774/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:3372159/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:3385653/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1845661/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1825013/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1727068/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:76818/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:882292/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:780535/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1484908/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1718429/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:730524/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1515948/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:604868/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:1669200/69‑1 (MQ=255)
tttGTCGGTTTCAAAGCCGTGATTAACCACGTTGACAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:2872067/69‑1 (MQ=255)
         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)
                        aaCCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATg                                                     <  1:3009481/45‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:2885540/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:1194619/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:2771869/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:2693189/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:811581/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:1413812/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:1492647/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:2368874/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:232310/69‑1 (MQ=255)
                                      cGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTa               <  1:1649917/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:983302/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:1011143/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:1615416/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:446969/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:3372693/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2125855/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2157958/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:215852/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2600901/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2966518/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2342415/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:3162967/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2409273/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2744914/68‑1 (MQ=255)
                                             cGGCGAGAACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGtt         <  1:2526104/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:566260/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:968607/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:87173/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:1255864/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:502891/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:2090680/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:3000634/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCACCAGTTGCTGCtt  >  1:3065182/1‑69 (MQ=255)
                                                   gaACTTACCGAAAGTATGCAGTTCGATTTTGTCGCTCAGCGCCATCATCTCGCGCAACAGTTGCTGCtt  >  1:2306874/1‑69 (MQ=255)
                                                               |                                                        
TTTGTCGGTTTCAAAGCCGTGATTAACCACGTTGCCAGCGGTGAACGGCGAGAACTTACCAAAGGTATGCAGTTCAATTTTGTCGCCCAGCGCCATCATCTCGCGTACCAGTTGCTGGTT  >  minE/1339786‑1339905

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: