Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,445,728 A→G 7.4% D219D (GAT→GAC)  mepA ← murein DD‑endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,445,7280AG 7.4% 157.1 / 3.4 54D219D (GAT→GAC) mepAmurein DD‑endopeptidase
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base G (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGT  >  NC_000913/2445590‑2445870
                                                                                                                                          |                                                                                                                                              
gCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAgg                                                                                                                                        <  5:171537/147‑1 (MQ=255)
           gcAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCg                                                                                                                            >  3:89997/1‑148 (MQ=255)
           gcAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCg                                                                                                                            >  1:221154/1‑148 (MQ=255)
             aaCGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTGTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTAc                                                                                                                         <  6:104442/149‑1 (MQ=255)
               cGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTGCCCGGTTTTGGAGGGTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACAt                                                                                                                       <  1:222750/149‑1 (MQ=255)
                      ggTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATg                                                                                                                >  8:211774/1‑149 (MQ=255)
                            ttcttctCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATgcgcgc                                                                                                           >  8:53094/1‑148 (MQ=255)
                                    aGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTgg                                                                                                   >  5:119375/1‑148 (MQ=255)
                                       cTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAAc                                                                                                >  3:259561/1‑148 (MQ=255)
                                              gttCCCGGATTTGGAGGTTCAAACCAGCTTACCAGTTCAGACCCGCAACCATCGCCTGATGGCGGTAAAGCTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCg                                                                                        >  8:218939/1‑149 (MQ=255)
                                                  ccGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACt                                                                                    >  2:127512/1‑149 (MQ=255)
                                                   cGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACtt                                                                                   >  6:113604/1‑149 (MQ=255)
                                                        ttGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTgcg                                                                               <  1:56070/148‑1 (MQ=255)
                                                           gAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTGTCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACttt                                                                                  >  3:287494/1‑142 (MQ=255)
                                                           gAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTGTCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACttt                                                                                  >  3:292038/1‑142 (MQ=255)
                                                           gAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTGTCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACttt                                                                                  <  4:292038/142‑1 (MQ=255)
                                                           gAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTGTCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACttt                                                                                  <  4:287494/142‑1 (MQ=255)
                                                             ggTTCAAACCAGCTTTGCAGTTCTGCCGCGCATCCTTCGCCTTATGTCGGTAAAGGTTGATCTTCTCTCTGCAGAGTATCGGCATGGCAACGTAATCGTACATGCATATGGGCTCGATTCTGGAACCAGGGTCGAACTTTgtgg                                                                              <  7:60762/144‑4 (MQ=255)
                                                             ggTTCAAACAAGCATTGCAGTTCTCCCCCGCAAACATCGCCAGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTAAAGGAAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTTGAACCAGGGTCGCACtttcagc                                                                              >  8:60762/1‑140 (MQ=255)
                                                               ttCACACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATTTGAGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAGCCAGt                                                                       >  4:62963/1‑149 (MQ=255)
                                                               ttCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGt                                                                       >  3:186771/1‑149 (MQ=255)
                                                                   aaCCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTcgc                                                                    >  8:42449/1‑148 (MQ=255)
                                                                     ccAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGAt                                                                 >  2:148668/1‑149 (MQ=255)
                                                                       aGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCg                                                               >  6:127522/1‑149 (MQ=255)
                                                                         cTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGt                                                             <  8:206693/149‑1 (MQ=255)
                                                                         cTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGt                                                             <  2:221154/149‑1 (MQ=255)
                                                                         cTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGt                                                             <  4:89997/149‑1 (MQ=255)
                                                                          tttGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGATGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTg                                                            >  5:83763/1‑149 (MQ=255)
                                                                                 ttCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCCGTCGCGATCGGTGCCCGCAt                                                     >  6:32927/1‑149 (MQ=255)
                                                                                      ccccGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAgg                                                >  3:140164/1‑149 (MQ=255)
                                                                                      ccccGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAg                                                 <  7:53094/148‑1 (MQ=255)
                                                                                            aaCCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgt                                          <  4:140164/149‑1 (MQ=255)
                                                                                             aCCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgtt                                         <  3:306440/149‑1 (MQ=255)
                                                                                              ccATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgttg                                        <  5:32927/149‑1 (MQ=255)
                                                                                                   gCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgttgttgtt                                   <  3:301130/149‑1 (MQ=255)
                                                                                                       gATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAAt                               >  1:134983/1‑149 (MQ=255)
                                                                                                        aTGGGGGTAAAGGTTGATCTTCACCCTCCAGACTAGCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATGGGGGCCCGCATCAAGGCAAAGTTGTTGTTTAATc                              <  3:48993/149‑1 (MQ=255)
                                                                                                           gCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGc                            <  3:275574/148‑1 (MQ=255)
                                                                                                            cGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAAc                                                                           >  8:116875/1‑100 (MQ=255)
                                                                                                            cGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAAc                                                                           <  7:116875/100‑1 (MQ=255)
                                                                                                                   ggTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTaa                    >  7:146372/1‑148 (MQ=255)
                                                                                                                    gTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACCTGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAAc                   >  6:99655/1‑148 (MQ=255)
                                                                                                                    gTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAAc                   >  7:29764/1‑148 (MQ=255)
                                                                                                                     ttGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACaa                 <  7:211774/149‑1 (MQ=255)
                                                                                                                     ttGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACa                  >  1:311431/1‑148 (MQ=255)
                                                                                                                        aTCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTTGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGGAAAGTTGTTGTTTAATCGCCGGATTAACaaaa               <  4:259561/148‑1 (MQ=255)
                                                                                                                         tCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACaaaaa              <  5:241319/148‑1 (MQ=255)
                                                                                                                             cacaCTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACaaaaa              <  1:127512/144‑1 (MQ=255)
                                                                                                                              acacTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGGACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTATCAAAAATgcgc         <  6:83763/148‑1 (MQ=255)
                                                                                                                                  tCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGAc    >  4:229112/1‑149 (MQ=255)
                                                                                                                                  tCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGAc    >  8:26838/1‑149 (MQ=255)
                                                                                                                                    cAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgttgtt                                      <  4:206092/113‑1 (MQ=255)
                                                                                                                                    cAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAgttgtt                                      >  3:206092/1‑113 (MQ=255)
                                                                                                                                    cAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGt  >  2:136684/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTCCAGACTATCGGCAGGGCAACGTAATCGTACATGCATATGCGCGCGATGCTGGAACCAGGGTCGCACTTTGCGCAACCAGTCGCGATCGGTGCCCGCATCAAGGCAAAGTTGTTGTTTAATCGCCGGATTAACAAAAATGCGCGTGACGT  >  NC_000913/2445590‑2445870

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