Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 212,234 T→C 59.6% K182K (AAA→AAG mltD ← predicted membrane‑bound lytic murein transglycosylase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE212,2340TC59.6% 15.4 / 25.9 42K182K (AAA→AAGmltDpredicted membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (2/15);  new base C (9/16);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATA  >  minE/212166‑212302
                                                                    |                                                                    
gCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACAtt                                                                      >  1:907229/1‑69 (MQ=255)
gCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACAtt                                                                      >  1:2902857/1‑69 (MQ=255)
gCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACAtt                                                                      >  1:1677860/1‑69 (MQ=255)
gCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACAtt                                                                      >  1:17494/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGGCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:3029874/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTg                                                        <  1:1406659/69‑1 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTg                                                        <  1:1517111/69‑1 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTg                                                        <  1:28734/69‑1 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTg                                                        <  1:2613400/69‑1 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTg                                                        <  1:2313365/69‑1 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:982602/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:832907/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:797708/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:484135/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:3178284/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:2720800/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:1253726/1‑69 (MQ=255)
              aCCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTg                                                        >  1:1731751/1‑69 (MQ=255)
                cTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGc                                                       >  1:1577969/1‑68 (MQ=255)
                cTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGc                                                       >  1:489459/1‑68 (MQ=255)
                     ccgcTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGcatca                                                   <  1:2145362/67‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:1017316/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:3363164/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:2897187/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:1565286/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:3413894/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:2645580/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:608742/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:2414716/69‑1 (MQ=255)
                                   cGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAgt                                   <  1:2246541/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1176920/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1076429/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:978061/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1977611/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1801667/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1430697/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:3091989/69‑1 (MQ=255)
                                          gTCAGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGaa                            <  1:1266560/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:900060/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:1492274/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:2737483/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:2671511/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:251701/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:2334127/69‑1 (MQ=255)
                                             aGAAGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:1842316/69‑1 (MQ=255)
                                               aaGCCAGTCGCCGTCAAACATCTTGTTCAGACGCTGCATCATGTTCAGCGCGGCAGTTGTTGAAGca                         <  1:2054441/67‑1 (MQ=255)
                                                                   ttttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCAt   <  1:3190748/69‑1 (MQ=255)
                                                                   ttttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCAt   <  1:350830/69‑1 (MQ=255)
                                                                   ttttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCAt   <  1:960756/69‑1 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCAt   <  1:1733507/68‑1 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATa  >  1:3223766/1‑69 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATa  >  1:3130922/1‑69 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATa  >  1:2652368/1‑69 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATa  >  1:1198862/1‑69 (MQ=255)
                                                                    tttGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATa  >  1:2043446/1‑69 (MQ=255)
                                                                    |                                                                    
GCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAACATTTTGTTCAGACGCTGCATCATATTCAGCGCGGCAGTTGTTGAAGCAACAACATCGCGACGCGCGTCATA  >  minE/212166‑212302

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: