Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,940,584 C→A 100% R351L (CGT→CTT)  mepM ← peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,940,5840CA97.7% 100.8 / ‑6.4 38R351L (CGT→CTT) mepMpeptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (18/17);  total (20/18)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGATTAG  >  NC_000913/1940446‑1940725
                                                                                                                                          |                                                                                                                                             
tCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAAt                                                                                                                                           <  1:780662/143‑1 (MQ=255)
    actGCAGACGCCGCCCGGTTGAACGCCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGATTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGcc                                                                                                                                       <  1:177365/141‑1 (MQ=255)
    aGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGCAGAATAGcc                                                                                                                                       >  1:447375/1‑143 (MQ=255)
     gTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGACGCCCCGTTTCACCTTCTGTCCCGGTTTCACCAGAATCATGCGCAAGTGCCTATAACGCGTGGCGCCGCTGCGACCATGCAGAACAGCca                                                                                                                                      >  1:524561/1‑143 (MQ=255)
             gCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACAtaataa                                                                                                                              <  1:537707/143‑1 (MQ=255)
              cGGCCCGGGTGAACGCCCGGTATCACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCAGGAAGAATAGCCACataagatc                                                                                                                             <  2:787867/143‑5 (MQ=255)
                gCCCGGTTGAACGGCCGGTCTTACCGGAAAGCGCGAGACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCGCATAATAACct                                                                                                                           <  1:558194/143‑1 (MQ=255)
                 cccGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACctg                                                                                                                          <  2:120479/143‑1 (MQ=255)
                   cGGTTGAACGTCCGGTATTACCGGAAAGCGCGCTACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAGTCTTGCGCAAGTGCATATAACGCGTGGCGTAGCTGCGACCATGACGAATAGCCACACAATAGCcggct                                                                                                                        >  2:527393/1‑143 (MQ=255)
                    ggTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACctgctg                                                                                                                       <  1:446361/143‑1 (MQ=255)
                        gAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGCTTCACCAGAATCTTGCGCAAGTGCATATAACTCGTGGTGTAGCTGCGACCATGAAGAATAGCACCAGAATCACCTGCCACGcc                                                                                                                   >  2:846627/1‑143 (MQ=255)
                         aaCGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATATTTCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCa                                                                                                                  <  2:824160/143‑1 (MQ=255)
                                 gaaTTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATACAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGttt                                                                                                          <  2:145239/141‑1 (MQ=255)
                                    ttACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGc                                                                                                       <  1:616810/143‑1 (MQ=255)
                                      aCCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGCTGTAGCTGCGACCCTGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCaa                                                                                                     >  1:479276/1‑143 (MQ=255)
                                                 cgttACGTTCGCCACGGTTCACCTTCTGGCCCGTTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCCCTTCa                                                                                          <  2:743864/140‑1 (MQ=255)
                                                            ccACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCAc                                                                               <  1:196364/143‑1 (MQ=255)
                                                                  ttCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATACCGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTCGGCAACCACCACTTCACCGTCACCCACTGCAAg                                                                         >  2:50270/1‑143 (MQ=255)
                                                                  ttCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGCAAg                                                                         >  2:14457/1‑143 (MQ=255)
                                                                           tGTCCGGGTTTCCCCAGAAGCTTGCGCAAGTGCATATAGCGGGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGt                                                                <  2:835256/143‑1 (MQ=255)
                                                                             tCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGCCCATGACGAATAGCCACATCATAACCTGCTGCGCCACTGCGTTTGGCGACCACCGCTTCACCGTCACCCACGGACAGACCTGGCGTAc                                                              >  1:152030/1‑143 (MQ=255)
                                                                                cGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGCCGTACCCt                                                           >  2:793598/1‑143 (MQ=255)
                                                                                      cACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCa                                                     >  1:22787/1‑143 (MQ=255)
                                                                                              tCTTGCGCGAGGGCGTATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGaaa                                             <  2:447375/143‑1 (MQ=255)
                                                                                                ttGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATc                                           >  1:79179/1‑143 (MQ=255)
                                                                                                  gcgcAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCaa                                         <  2:744284/143‑1 (MQ=255)
                                                                                                     cAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAacac                                      <  1:422754/143‑1 (MQ=255)
                                                                                                     cAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACAGTAACCCACTGATAGCACTGGCGTACCCTGCGGCATGGCGAAATCAacac                                      <  1:803900/143‑1 (MQ=255)
                                                                                                        gTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGCAATCAACACctc                                   >  1:80876/1‑143 (MQ=255)
                                                                                                             tataACGCGTGGTGTAGCTGCGACAATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGCTTGGCAACCACCACTACAACGTCACCCACTGAAAGCACTGGCCTACCCCGCGGCACGGCGCAATCACCACCTAcgtgt                              >  2:352807/1‑143 (MQ=255)
                                                                                                               taACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTgtgg                            >  1:728001/1‑143 (MQ=255)
                                                                                                                                   aCCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCCACGCGACCGGTCACCg        >  2:796845/1‑143 (MQ=255)
                                                                                                                                   aCCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGGACCCTGCGGCATGACGAACTCACCCCCTCACCGTGGTCGAACGCGACCGGTCCCCg        >  2:346925/1‑143 (MQ=255)
                                                                                                                                    ccATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCCCTTCACCGTCACCCACTGACAGCACTGGCGTACCCCCCGCCATGGCGCACCCAACACCTCTGTGCCGTGCAAAGCCACCGGTCAcctg       >  1:806138/1‑141 (MQ=255)
                                                                                                                                    ccATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCCCTGCAAGCACTGGCGTACCCTGCGGCCTGGCGAAATCAACAACACACCGGCGCGCAACGCGACCGGTCCCCgg       >  1:751947/1‑143 (MQ=255)
                                                                                                                                      aTGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAAACACCACTTCACAGACACCCACTGCAAGCACTGGCGTACCCTGCAGCATGAGGAAGTCAAAGCGGCCGCGTGGGGCAACGCGCCCGAGCACCggct     >  2:390542/1‑141 (MQ=255)
                                                                                                                                       tGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGACAACCACCACTGCACCGTCACCCACTGAAAGCACTGGCGTACCCCGCGGCATGGCGAAATCAACACCCCCGTGCGGTGCAACGCGACCGGGCAcacgcac    >  2:284756/1‑137 (MQ=255)
                                                                                                                                         aCGAATACCTACCTAACCTCCAGCTGCTCCACTTCGTTTGGCGACCACCACTTCACCGTCAGCCGGTGAAAGCGCTGCCTGACCCTGCGGCATGGCGAAATGAACGCCTCTGTGTGGTGCAACGCGACCGGTCACCGGATTAg  <  2:337217/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
TCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGATTAG  >  NC_000913/1940446‑1940725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: