Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 531,622 G→A 34.0% N314N (AAC→AAT ybhR ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE531,6220GA34.0% 63.9 / 43.4 82N314N (AAC→AATybhRpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/39);  new base A (5/23);  total (20/62)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGA  >  minE/531554‑531676
                                                                    |                                                      
cTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAgg                                                        <  1:1095705/69‑1 (MQ=255)
          tCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGc                                               <  1:3374154/68‑1 (MQ=255)
          tCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGc                                               <  1:333625/68‑1 (MQ=255)
          tCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGc                                               <  1:18199/68‑1 (MQ=255)
          tCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGc                                               <  1:2518716/68‑1 (MQ=255)
          tCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGc                                               <  1:3099832/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAATGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:2803570/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:2638522/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:3384902/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:217876/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:241687/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:3306632/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:687778/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:1511180/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:798771/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:2677337/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:3063983/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:135230/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:283376/68‑1 (MQ=255)
             aataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCAt                                            <  1:1024842/68‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:2104561/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3084640/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3208846/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3080591/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3323978/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3374484/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3012169/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:173130/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:978665/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1104317/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1208409/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:91331/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1416025/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:800112/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1532245/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1689695/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:2264323/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:765286/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:621747/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:1976073/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3425476/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:3419932/63‑1 (MQ=255)
               taaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGATTTTGCAGc                                               <  1:2208103/63‑1 (MQ=255)
                         tGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGttt                                >  1:1572977/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCgg                        >  1:577780/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCgg                        >  1:101090/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCgg                        >  1:1773858/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCgg                        >  1:2428682/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCg                         >  1:611775/1‑67 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCg                         >  1:2471681/1‑67 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCAACgg                        >  1:3030446/1‑68 (MQ=255)
                                 tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCACACCGGCATGTTTTCCACCgg                        >  1:1455968/1‑68 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCAGCGGAGaa                    <  1:2934289/68‑1 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGaaa                   <  1:637611/69‑1 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGaaa                   <  1:2502863/69‑1 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGaaa                   <  1:2981648/69‑1 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGaa                    <  1:2275050/68‑1 (MQ=255)
                                     tAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGaa                    <  1:2044105/68‑1 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTaa              >  1:2484722/1‑69 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTaa              >  1:2581649/1‑69 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTaa              >  1:611131/1‑69 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTaa              >  1:2650752/1‑69 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTaa              >  1:1000678/1‑69 (MQ=255)
                                          tGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTa               >  1:814404/1‑68 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGATAAACGTAACCGGAAAGgag   <  1:3121150/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:1796088/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:533141/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:523373/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:503550/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:416336/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:1391523/69‑1 (MQ=255)
                                                     gTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:3266341/69‑1 (MQ=255)
                                                      ttAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:2578937/68‑1 (MQ=255)
                                                      ttAATCCACGTCAGATTTTGCAGCCACATTGGCATGTTTTCCACCGGCGAAACGTAACCGGAAAGgaga  >  1:984039/1‑69 (MQ=255)
                                                      ttAATCCACGTCAGATTTTGCAGCCACATTGGCATGTTTTCCACCGGCGAAACGTAACCGGAAAGgaga  >  1:615197/1‑69 (MQ=255)
                                                      ttAATCCACGTCAGATTTTGCAGCCACATTGGCATGTTTTCCACCGGCGAAACGTAACCGGAAAGgaga  >  1:2080562/1‑69 (MQ=255)
                                                      ttAATCCACGTCAGATTTTGCAGCCACATTGGCATGTTTTCCACCGGCGAAACGTAACCGGAAAGgaga  >  1:2794024/1‑69 (MQ=255)
                                                      ttAATCCACGTCAGATTTTGCAGCCACATTGGCATGTTTTCCACCGGCGAAACGTAACCGGAAAGgag   >  1:2651380/1‑68 (MQ=255)
                                                                  aGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:1826390/56‑1 (MQ=255)
                                                                  aGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:3428087/56‑1 (MQ=255)
                                                                  aGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:1425120/56‑1 (MQ=255)
                                                                  aGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:220230/56‑1 (MQ=255)
                                                                  aGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGgag   <  1:2402426/56‑1 (MQ=255)
                                                                    |                                                      
CTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGA  >  minE/531554‑531676

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: