Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,496,207 G→A 45.2% N39N (AAC→AAT emrD ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,496,2070GA45.2% 20.0 / 76.4 95N39N (AAC→AATemrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (24/19);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

ACCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAC  >  minE/2496153‑2496275
                                                      |                                                                    
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:2983124/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:1439406/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:1582270/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:2668440/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:272096/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:2977487/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:3051119/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:3222731/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:3315637/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:447707/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:606716/69‑1 (MQ=255)
acCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:80271/69‑1 (MQ=255)
 ccGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:959097/68‑1 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGGCATTgag                                                                 >  1:1156126/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:2414938/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:2562684/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:2087375/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1791695/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1649974/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:989324/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1408006/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1346535/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1236707/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1175102/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1152618/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1146690/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1101212/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1362115/1‑58 (MQ=255)
  cGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTgag                                                                 >  1:1052413/1‑58 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCa                                              <  1:1412213/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:685801/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:3317014/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:1013667/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:819103/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:1154781/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:964935/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:10113/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:2341689/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:206808/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:1516071/69‑1 (MQ=255)
          agcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:1687975/69‑1 (MQ=255)
           gcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCa                                              <  1:3269293/68‑1 (MQ=255)
           gcagATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCa                                              <  1:2777667/68‑1 (MQ=255)
             agATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCgcgcgc                                                        <  1:487964/56‑1 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:2608014/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:513532/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:1297103/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:1704296/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:3005465/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:2143094/1‑69 (MQ=255)
              gATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAg                                          >  1:85437/1‑69 (MQ=255)
                      cccATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCa                                              <  1:1334129/57‑1 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2039110/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2987416/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:3326239/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1737846/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:375281/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1666045/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:568530/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1634273/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:66821/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:160881/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1384764/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1545018/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1012277/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1495383/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2945903/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2889096/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2027418/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2068112/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1343977/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2812083/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:277496/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2287864/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:1009205/1‑69 (MQ=255)
                            aCGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATg                            >  1:2515832/1‑69 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:2508149/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:204364/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:805997/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:1860529/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:1338903/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:2664729/68‑1 (MQ=255)
                               ctctGCACCGCCCCTTCACGGACATTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGt                          <  1:2757227/68‑1 (MQ=255)
                                 ctGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGc                                               >  1:81878/1‑45 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:2981926/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:3013556/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:3030895/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:831160/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:2187188/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:2283463/69‑1 (MQ=255)
                                      ccGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCa                  <  1:1791368/69‑1 (MQ=255)
                                        gCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCAt                 >  1:2360340/1‑68 (MQ=255)
                                               cACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGa         >  1:1687864/1‑69 (MQ=255)
                                               cACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATATGACCGa         >  1:1720608/1‑69 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:645881/69‑1 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:825493/69‑1 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:897036/69‑1 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:490791/69‑1 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:2878011/69‑1 (MQ=255)
                                                      gTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAc  <  1:2562386/69‑1 (MQ=255)
                                                      |                                                                    
ACCGTAAGTCAGCAGATAAGCGCCCATTACGCTCTGCACCGCCCCTTCACGGACGTTGAGATCGCGCGCCATATCGGCAATAGCTGGAATATAAATGGTTTGCGCCATCTGACCGACGGCCAC  >  minE/2496153‑2496275

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: