Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,696,696 C→T 31.6% Q510Q (CAG→CAA mrcA ← fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,696,6960CT31.6% 95.5 / 50.9 104Q510Q (CAG→CAAmrcAfused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base T (16/17);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 113 total aligned reads displayed.

GCGCGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATC  >  minE/2696631‑2696764
                                                                 |                                                                    
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:10602/69‑1 (MQ=255)
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:759774/69‑1 (MQ=255)
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:3274381/69‑1 (MQ=255)
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:642536/69‑1 (MQ=255)
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:3329325/69‑1 (MQ=255)
gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:1283268/69‑1 (MQ=255)
        cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           >  1:2043758/1‑69 (MQ=255)
        cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           >  1:1320143/1‑69 (MQ=255)
        cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           >  1:464153/1‑69 (MQ=255)
         aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           <  1:1677540/68‑1 (MQ=255)
         aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           <  1:2314778/68‑1 (MQ=255)
         aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           <  1:624870/68‑1 (MQ=255)
         aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           <  1:885888/68‑1 (MQ=255)
         aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                           <  1:1022318/68‑1 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:653696/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:2939019/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:956835/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1962345/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1862687/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1780084/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1032837/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1033436/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:1270127/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                                       >  1:148810/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                                        >  1:284918/1‑68 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                                        >  1:738884/1‑68 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                                        >  1:1610244/1‑68 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:871955/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:2386963/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:2549437/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:2572140/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:1440318/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:1273780/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                                       >  1:1600370/1‑69 (MQ=255)
            accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGAACGGCATATTGCGGTGGTGAGttc                                                       >  1:1707903/1‑69 (MQ=255)
                        gACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCa                                            <  1:651599/68‑1 (MQ=255)
                          gtGACCCAGCGCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:148984/42‑1 (MQ=255)
                           tGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGt                                                         >  1:1713489/1‑52 (MQ=255)
                           tGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCt                                                >  1:3129449/1‑61 (MQ=255)
                           tGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTc                                         >  1:3350454/1‑68 (MQ=255)
                           tGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTc                                         >  1:3271788/1‑68 (MQ=255)
                            gACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       <  1:408389/69‑1 (MQ=255)
                             aCCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:1075196/40‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:649411/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:781585/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:467922/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:2735250/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:2959156/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:285271/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1160910/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1623943/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:184962/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1712051/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:2285935/69‑1 (MQ=255)
                               ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1393789/69‑1 (MQ=255)
                                cAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                                   <  1:2664369/37‑1 (MQ=255)
                                        gACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTg                           >  1:425310/1‑69 (MQ=255)
                                        gACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTg                           >  1:1962084/1‑69 (MQ=255)
                                        gACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTg                           >  1:1438410/1‑69 (MQ=255)
                                         aCGTAAGCGAATTGGACCGGCATATTGCGGTCGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1895685/59‑1 (MQ=255)
                                          cGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:2948878/58‑1 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:1839586/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:1215833/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:2062776/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:90834/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:387500/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:196335/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:207290/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                                       >  1:274783/1‑51 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc                     >  1:1444408/1‑69 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc                     >  1:3336147/1‑69 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc                     >  1:902801/1‑69 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc                     >  1:2230526/1‑69 (MQ=255)
                                              aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATcc                      >  1:2911393/1‑68 (MQ=255)
                                              aGCGAATTGGACCAGCATACCGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc                     >  1:314609/1‑69 (MQ=255)
                                                   aTTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:34653/49‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:249415/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2637162/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2690803/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:26990/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:1442447/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:1892986/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2215335/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2955146/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:1287756/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:1316703/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2958366/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2414963/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:3003605/69‑1 (MQ=255)
                                                      ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:2998976/69‑1 (MQ=255)
                                                       gACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa             <  1:1486307/68‑1 (MQ=255)
                                                        aCCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                                    <  1:1762111/44‑1 (MQ=38)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:1566586/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:760742/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:543449/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:42357/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:103479/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:1379660/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:1880947/1‑69 (MQ=255)
                                                                 cTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATc  >  1:950620/1‑69 (MQ=255)
                                                                 |                                                                    
GCGCGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAAATTGGCACATC  >  minE/2696631‑2696764

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: