Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,696,690 A→G 31.0% A512A (GCT→GCC 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,6900AG31.0% 95.1 / 52.3 106A512A (GCT→GCCmrcAfused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (16/17);  total (52/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 = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CGCATTGCGCGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAA  >  minE/2696625‑2696753
                                                                 |                                                               
cGCATTGCGCGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCa                                                              >  1:2548637/1‑69 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:642536/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:3274381/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:3329325/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:3382186/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:1283268/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:2275795/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:759774/69‑1 (MQ=255)
      gcgcGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:10602/69‑1 (MQ=255)
              cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTg                                                         >  1:773690/1‑60 (MQ=255)
              cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                >  1:464153/1‑69 (MQ=255)
              cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                >  1:1320143/1‑69 (MQ=255)
              cATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                >  1:2043758/1‑69 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:1718317/68‑1 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:2314778/68‑1 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:1573564/68‑1 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:1677540/68‑1 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:3038695/68‑1 (MQ=255)
               aTCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGa                                                <  1:1429752/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:2939019/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                            >  1:653696/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:959390/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                            >  1:1032837/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                            >  1:1033436/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                            >  1:148810/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGttc                                            >  1:1270127/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                             >  1:284918/1‑68 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                             >  1:1610244/1‑68 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                             >  1:738884/1‑68 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGtt                                             >  1:567076/1‑68 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:1273780/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:2894077/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:2572140/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:2549437/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:2386963/1‑69 (MQ=255)
                  accacGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGttc                                            >  1:871955/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:1160910/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:781585/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:1712051/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:1393789/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:2959156/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:649411/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:285271/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:2735250/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:1623943/69‑1 (MQ=255)
                                     ccAGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:2285935/69‑1 (MQ=255)
                                      cAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGc                                                        <  1:2664369/37‑1 (MQ=255)
                                       aGCCCCTGACGTAAGCGAATTGGACCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:1233564/67‑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:274783/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:387500/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:207290/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:2062776/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:196335/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:1215833/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:90834/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:1839586/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAg                            >  1:2607465/1‑51 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc          >  1:3336147/1‑69 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc          >  1:902801/1‑69 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc          >  1:2230526/1‑69 (MQ=255)
                                                    aGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATccc          >  1:1444408/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:1442447/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:26990/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:1892986/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2955146/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2690803/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2998976/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:3003605/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:1316703/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2637162/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2215335/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:1287756/69‑1 (MQ=255)
                                                            ggACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:2414963/69‑1 (MQ=255)
                                                             gACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGaa  <  1:1486307/68‑1 (MQ=255)
                                                              aCCGGCATATTGCGGTGGTGAGTTCTTCGGTTGCCAGTCAGAAc                         <  1:1762111/44‑1 (MQ=38)
                                                                 |                                                               
CGCATTGCGCGTACCATCACCACGTTTTTCGACTGACCCAGCCCCTGACGTAAGCGAATTGGACCAGCATACTGCGGTGGTGAGTTCTTCGGCTGCCAGTCAGAACCGGCACTTGCATCCCAGCGAGAA  >  minE/2696625‑2696753

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: