Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,735,191 T→C Y253H (TAC→CAC)  mepH → murein DD‑endopeptidase, space‑maker hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,735,1910TC98.0% 185.2 / NA 50Y253H (TAC→CAC) mepHmurein DD‑endopeptidase, space‑maker hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/23);  minor base A (0/1);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAG  >  NC_000913/1735047‑1735336
                                                                                                                                                |                                                                                                                                                 
gAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGc                                                                                                                                             >  2:805536/1‑151 (MQ=255)
 agcTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTAACCGCGCACAGGTCAGGAAATTCAAAGCACTTCTCTCAGTGAAGACCACTGGc                                                                                                                                             <  2:279443/148‑1 (MQ=255)
 aaCTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCa                                                                                                                                            <  1:451906/151‑1 (MQ=255)
       aaaaCGGCGACCTGGTCCTTTTCGGTACTCAGGGACGCGGCCCAGCGGATCATGTCGGCGTGTATGTCGGCATCGGCAAATTTTTTCAGTCACCGCGCACAGGTCCGGAAATTCAAATCTCTTCTCTCAGTGAAGACCACTTGCAgc                                                                                                                                          <  2:873568/147‑1 (MQ=255)
                   tGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGgcgc                                                                                                                          >  2:878185/1‑151 (MQ=255)
                    ggTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCt                                                                                                                         >  2:21358/1‑151 (MQ=255)
                      tCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcg                                                                                                                       >  2:182549/1‑151 (MQ=255)
                      tCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcg                                                                                                                       >  2:553744/1‑151 (MQ=255)
                          tttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGTGCACAGGTTAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCgg                                                                                                                   >  1:316933/1‑151 (MQ=255)
                            tCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGt                                                                                                                 >  1:599724/1‑151 (MQ=255)
                            tCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGt                                                                                                                 >  2:516815/1‑151 (MQ=255)
                             ccGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTa                                                                                                                <  2:625046/151‑1 (MQ=255)
                                   tCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTAAGGAAATTCAAATCACTTCTCTCAGTGAAGACCGCTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTa                                                                                                                <  2:253754/145‑1 (MQ=255)
                                    cAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGAccc                                                                                                         >  2:867389/1‑151 (MQ=255)
                                          cgcgTCACACCCGCTAATGTCGGCGTGTCTGTCGGCAGCGGTAAACTTATTCAGTCACCGCGCGTAGGTCAGGAAATTCAAACCTCTTCTCTCAGTGAAGACAACTTTCAGCGCCGCAAAGTTGGCGCTCGTCGGGTAAt                                                                                                              <  2:159879/140‑1 (MQ=255)
                                           gcgGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaaaa                                                                                                   >  2:734360/1‑150 (MQ=255)
                                                 caGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCg                                                                                            <  2:858025/151‑1 (MQ=255)
                                                   gCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCt                                                                                          >  1:676193/1‑151 (MQ=255)
                                                       aTCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTaaaa                                                                                      <  1:517331/151‑1 (MQ=255)
                                                        tCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAAc                                                                                     >  2:27618/1‑151 (MQ=255)
                                                          aTGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACtt                                                                                   >  1:810054/1‑151 (MQ=255)
                                                             tCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAc                                                                                >  1:651968/1‑151 (MQ=255)
                                                               ggCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAccc                                                                              <  1:734401/151‑1 (MQ=255)
                                                               ggCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAccc                                                                              >  1:495615/1‑151 (MQ=255)
                                                               ggCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAcca                                                                              <  2:947270/151‑2 (MQ=255)
                                                                                cggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAg                                                             <  1:516842/151‑1 (MQ=255)
                                                                                 ggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAgg                                                            <  1:287042/151‑1 (MQ=255)
                                                                                 ggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAgg                                                            >  1:639333/1‑151 (MQ=255)
                                                                                    aaaTTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTa                                                         <  2:27061/151‑1 (MQ=255)
                                                                                      aTTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTGCCTAAAATTTTACCTTGTTGTTACGCCAACAgggaa                                                         >  1:874566/1‑147 (MQ=255)
                                                                                      aTTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTa                                                         >  2:756844/1‑149 (MQ=255)
                                                                                        ttATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGt                                                      <  1:357861/150‑1 (MQ=255)
                                                                                          aTTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTc                                                    <  1:59711/150‑1 (MQ=255)
                                                                                            tCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATc                                                 >  1:773723/1‑151 (MQ=255)
                                                                                             cAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATca                                                <  2:323547/151‑2 (MQ=255)
                                                                                               gTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCttt                                              <  2:160027/151‑1 (MQ=255)
                                                                                                tCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCtttt                                             >  1:863363/1‑151 (MQ=255)
                                                                                                   ccGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtc                                          >  1:926670/1‑151 (MQ=255)
                                                                                                         acaGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCt                                    <  2:810008/151‑1 (MQ=255)
                                                                                                          caGGCCAGGAAATTCAAATCACTTCTCTCCGCGAAGACCACTGTCCGCGCCACTATGTAGGAGCTCGTCGGGTAATCACCCCAAAAACACGTCGCTAACACTGTACCCTGGTGTGACGGCAATAGGGTAACTTCATCTTTTGTCTCACCtt                                   >  2:1050457/1‑151 (MQ=255)
                                                                                                           aGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCttt                                  <  2:651933/151‑1 (MQ=255)
                                                                                                           aGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCttt                                  >  1:606632/1‑151 (MQ=255)
                                                                                                                    aaTTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCt                         <  1:602959/151‑1 (MQ=255)
                                                                                                                      ttCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTa                        <  2:19913/150‑1 (MQ=255)
                                                                                                                       tCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTAcc                      >  2:135500/1‑151 (MQ=255)
                                                                                                                       tCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTAc                       <  2:916119/150‑1 (MQ=255)
                                                                                                                                 tctctcAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATAc            >  1:609880/1‑151 (MQ=255)
                                                                                                                                    ctcAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATAc            <  1:706317/148‑1 (MQ=255)
                                                                                                                                    ctcAGTGAAGAC‑ACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACg           <  1:391040/148‑1 (MQ=255)
                                                                                                                                           aaGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAg  <  1:234418/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAG  >  NC_000913/1735047‑1735336

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑