Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 280,674 C→T 69.9% T447T (ACG→ACA)  ampG ← muropeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE280,6740CT69.9% 78.1 / 40.3 86T447T (ACG→ACA) ampGmuropeptide transporter
Reads supporting (aligned to +/- strand):  ref base C (7/19);  new base T (7/53);  total (14/72)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTC  >  minE/280609‑280725
                                                                 |                                                   
aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTc                                                  <  1:436405/69‑1 (MQ=255)
aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTc                                                  <  1:743080/69‑1 (MQ=255)
aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTc                                                  <  1:270779/69‑1 (MQ=255)
    aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAAt                                               >  1:1045771/1‑68 (MQ=255)
    aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAAt                                               >  1:3363924/1‑68 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTc                                                  >  1:215375/1‑63 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:3160655/1‑64 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:543842/1‑64 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:533714/1‑64 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:3356141/1‑64 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:121446/1‑64 (MQ=255)
      aCCCCGACTTCCAGCAGCGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATTGTCa                                                 >  1:1131309/1‑64 (MQ=255)
       ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                           <  1:3219696/69‑1 (MQ=255)
       ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                           <  1:2110344/69‑1 (MQ=255)
       ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                           <  1:3399772/69‑1 (MQ=255)
       ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                           <  1:2603145/69‑1 (MQ=255)
       ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                           <  1:1619579/69‑1 (MQ=255)
                         caggcaggAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                         >  1:1177924/1‑69 (MQ=255)
                         caggcaggAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                         >  1:1954699/1‑69 (MQ=255)
                         caggcaggAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                         >  1:285426/1‑69 (MQ=255)
                            gcaggAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTc                                                  <  1:944258/41‑1 (MQ=255)
                               ggAAAGAGAATTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:57510/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:3002621/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:3213644/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:3062367/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:2234009/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:2227097/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:277031/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:308673/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:214470/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1966774/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1971161/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1112671/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:115546/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1211647/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1229140/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:77560/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:763681/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:729104/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:723682/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:1780771/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:663665/69‑1 (MQ=255)
                               ggAAAGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCAAACGGCCAACAGGCTGa                   <  1:1798271/69‑1 (MQ=255)
                                 aaaGAGAAGTGCGCCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGCTGa                   <  1:2823559/67‑1 (MQ=255)
                                   agagAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCAGACGGCCAACAGGCTGAcgcc               <  1:2094732/69‑1 (MQ=255)
                                   agagAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc               <  1:1802814/69‑1 (MQ=255)
                                   agagAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc               <  1:1315470/69‑1 (MQ=255)
                                   agagAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc               <  1:3095193/69‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:715337/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:707248/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:81290/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:3295798/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:617016/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:839281/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:976102/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1015302/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2105298/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:104493/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1158559/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1199878/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1488257/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1713918/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1723942/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1743315/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1754500/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:1963780/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:3145986/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2137565/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2183877/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2967349/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2317350/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2890700/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:287670/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2329001/53‑1 (MQ=255)
                                            cGTCAAATCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:2366968/53‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:460054/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:1421227/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:1586174/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:1634126/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:333246/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  <  1:2723114/69‑1 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  >  1:183328/1‑69 (MQ=255)
                                                aaaTCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc  >  1:3436153/1‑69 (MQ=255)
                                                  aTCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCGAACAGGCTGa                   <  1:5050/50‑1 (MQ=255)
                                                  aTCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:3106103/47‑1 (MQ=39)
                                                   tCCAGCGCGTCCATTGTCAACAGCAGTAACCACACGGCCAACAGGc                      <  1:641346/46‑1 (MQ=39)
                                                                 |                                                   
ACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTC  >  minE/280609‑280725

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: ‑