Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 515,063 G→A 45.2% D236D (GAC→GAT mepM ← Murein DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440515,0630GA45.2% 8.0 / 53.6 42D236D (GAC→GATmepMMurein DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base A (6/13);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

AGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATT  >  KT2440/514933‑515192
                                                                                                                                  |                                                                                                                                 
aGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCg                                                                                                                                >  4:186747/1‑134 (MQ=255)
       gCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCa                                                                                                                         >  3:175619/1‑134 (MQ=255)
       gCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAg                                                                                                                        >  1:75971/1‑135 (MQ=255)
        ccGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAg                                                                                                                        <  7:59323/134‑1 (MQ=255)
          gcgGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCAc                                                                                                                     >  6:197829/1‑135 (MQ=255)
             gTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCAcgcg                                                                                                                  <  5:219855/135‑1 (MQ=255)
              ttGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCgg                                                                                                                 <  1:122145/135‑1 (MQ=255)
              ttGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCgg                                                                                                                 <  4:175619/135‑1 (MQ=255)
              gtGGGGAAGCGCGAAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGTACTCGTCCCCAGTGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCgg                                                                                                                 <  6:152623/134‑1 (MQ=255)
               tGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCggg                                                                                                                <  7:110573/135‑1 (MQ=255)
                ggTGAAGCGCGCAGATAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACACCGCACGCGGGc                                                                                                               >  8:142741/1‑135 (MQ=255)
                                   aTGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCtg                                                                                            >  5:208349/1‑135 (MQ=255)
                                           ggTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTATACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAg                                                                                    >  1:215589/1‑135 (MQ=255)
                                           ggTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAg                                                                                    >  1:230506/1‑135 (MQ=255)
                                           ggTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAg                                                                                    >  3:193856/1‑135 (MQ=255)
                                               ccGACGACCTTGCCATCCATCACCTTCTGTTCGTAGCCCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCCGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACCGGCCg                                                                                >  6:199712/1‑135 (MQ=255)
                                                 gacgacCTTGCCATCCATCTCCTTCTGTTAGTAGACCACGTCGAACTCGTCCCCCGGTCGAATATCATGTGCGACGTCTATATCGTCACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGgc                                                                              <  8:173226/135‑1 (MQ=255)
                                                      ccTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCt                                                                         <  7:293419/135‑1 (MQ=255)
                                                            catccatcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCgg                                                                   <  4:145818/135‑1 (MQ=255)
                                                            catccatcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCg                                                                    >  4:47556/1‑134 (MQ=255)
                                                               ccatcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAg                                                                <  7:204685/135‑1 (MQ=255)
                                                                 atcACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCGGTGGGACAGGCCGGCGCGCTGGGCGGAAGc                                                               <  8:30887/134‑1 (MQ=255)
                                                                     ccTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACa                                                          <  5:139025/135‑1 (MQ=255)
                                                                     ccTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACa                                                          <  7:53222/135‑1 (MQ=255)
                                                                          tGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGa                                                      <  7:27322/134‑1 (MQ=255)
                                                                          tGTTCGTAGACCACGACGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGGGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGa                                                      <  6:97633/134‑1 (MQ=255)
                                                                           gTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGCAGCAGACAGCGAGc                                                    >  7:60647/1‑135 (MQ=255)
                                                                           gTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGc                                                    >  3:122341/1‑135 (MQ=255)
                                                                                     cACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAc                                          <  8:197890/135‑1 (MQ=255)
                                                                                      ccGTCGACCTCTTCCCCAGGTCGAATATCCTGTGCGACGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCaac                                         <  8:34353/134‑3 (MQ=255)
                                                                                      aCGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcc                                         <  7:228598/135‑1 (MQ=255)
                                                                                      aCGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCAGGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcc                                         <  7:261280/135‑1 (MQ=255)
                                                                                      aCGGCGACCTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGCCCATGGTCATGCTGTGGGACTGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcc                                         <  2:115996/135‑1 (MQ=255)
                                                                                              cTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTggg                                  <  2:215589/134‑1 (MQ=255)
                                                                                                     ccAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATgg                                                                                                    <  3:123784/61‑1 (MQ=255)
                                                                                                     ccAGGGCGAATATCCTGTGCGAAGTCGATATCGTAACCCAGCACGCGGGCCATGTCCATgg                                                                                                    >  4:123784/1‑61 (MQ=255)
                                                                                                              atatCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCTGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcccc                                       <  1:264753/113‑1 (MQ=255)
                                                                                                              atatCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcccc                                       >  2:264753/1‑113 (MQ=255)
                                                                                                                      tgCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGTGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcc                                         <  1:299583/103‑1 (MQ=255)
                                                                                                                      tgCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGTGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCAcc                                         >  2:299583/1‑103 (MQ=255)
                                                                                                                       gCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTg        >  6:155037/1‑135 (MQ=255)
                                                                                                                             tCGATGTCGTAACCCAGCACGCGGGACATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGAtt  <  4:74225/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
AGGTTTTGCCGCGGTTGGTGAAGCGCGCAGAGAGGATGTTACCGGTGCCGACGACCTTGCCATCCATCACCTTCTGTTCGTAGACCACGTCGAACTCGTCCCCAGGGCGAATATCCTGTGCGAAGTCGATGTCGTAACCCAGCACGCGGGCCATGTCCATGGTCATGCTGTGGGACAGGCCGGCGCGCTGGGCGGAAGCAGACAGCGAGCTCTTGATCACCCCGTGGGCATAGGCGGTGCGTATCACCGGCTTGCTGATT  >  KT2440/514933‑515192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: