Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,413 A→G I515V (ATT→GTT)  mrcA → penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,524,4130AG94.6% 175.6 / ‑0.5 56I515V (ATT→GTT) mrcApenicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (30/23);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAAC  >  NC_000913/3524320‑3524509
                                                                                             |                                                                                                
aGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCt                                                                                            <  2:65157/100‑1 (MQ=255)
    gctcttccGATCTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACg                                                                                        <  1:757346/93‑1 (MQ=255)
     gTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGt                                                                                       >  1:1391640/1‑100 (MQ=255)
              tgtgCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTg                                                                              <  2:697102/100‑1 (MQ=255)
               gtgCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTgg                                                                             <  1:168214/100‑1 (MQ=255)
                tgCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTggg                                                                            <  1:1003668/100‑1 (MQ=255)
                tgCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTggg                                                                            <  2:472626/100‑1 (MQ=255)
                 gCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTggga                                                                           >  2:1003662/1‑99 (MQ=255)
                  ccAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTc                                                                          <  2:1003800/100‑1 (MQ=255)
                    aaTTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAg                                                                        <  2:1139784/100‑1 (MQ=255)
                     aTTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGt                                                                       <  1:1500963/100‑1 (MQ=255)
                      tttCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTc                                                                      <  1:1198643/100‑1 (MQ=255)
                        tctcGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGa                                                                    <  1:1216118/100‑1 (MQ=255)
                           cGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGaaaa                                                                 >  2:1566046/1‑100 (MQ=255)
                            gCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCAAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGaaaaa                                                                >  2:606675/1‑100 (MQ=255)
                             cTGGGATGCAAGTGCCGGTTCTGACTGGCAGCAGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAAc                                                               >  1:1497569/1‑100 (MQ=255)
                                 gATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACgtgg                                                           <  1:838908/100‑1 (MQ=255)
                                 gATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACgtgg                                                           <  2:625487/100‑1 (MQ=255)
                                      tcttGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATg                                                      <  1:911043/97‑1 (MQ=255)
                                      aaGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATg                                                      >  1:1019008/1‑100 (MQ=255)
                                       aGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAACCGTGGTGATgg                                                     >  1:176246/1‑100 (MQ=255)
                                        gTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGt                                                    >  1:97729/1‑100 (MQ=255)
                                         tGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATgaga                                                   >  2:911038/1‑97 (MQ=255)
                                           ccGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTAcg                                                 >  1:160832/1‑100 (MQ=255)
                                               ttCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCa                                             <  2:1210715/100‑1 (MQ=255)
                                                  tGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATg                                          >  2:102322/1‑100 (MQ=255)
                                                   gACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGc                                         >  2:391360/1‑100 (MQ=255)
                                                      tGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCggg                                      <  1:750104/100‑1 (MQ=255)
                                                       ggCAGCCGAAGAAGTCACCACCGCCGTATGCTGATCCAGTCCTCTTACGTCAGGGCCTGAGCCAGTCGAAAAACGAAGACATGGTACGCGCAATGCGGGc                                     >  2:555914/1‑100 (MQ=255)
                                                       ggCAGCCGAAGAACTCACCACCGGAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCggg                                      >  2:335555/1‑99 (MQ=255)
                                                          aGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTAGGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGc                                     <  1:673204/97‑1 (MQ=255)
                                                            ccGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATgg                                <  2:357665/100‑1 (MQ=255)
                                                             cGAAGAACTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATggg                               <  1:102322/100‑1 (MQ=255)
                                                                agaaCTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGt                            >  1:1489804/1‑100 (MQ=255)
                                                                 gaaCTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTc                           <  1:1566056/100‑1 (MQ=255)
                                                                 gaaCTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGCCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTc                           >  2:273333/1‑100 (MQ=255)
                                                                   aCTCACCACCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGa                         >  1:740575/1‑100 (MQ=255)
                                                                      caccacCGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGTTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTa                      >  2:182985/1‑100 (MQ=255)
                                                                        ccaccGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACg                    >  2:490211/1‑100 (MQ=255)
                                                                         caccGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGc                   >  2:390897/1‑100 (MQ=255)
                                                                           ccGCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTg                 >  1:1484657/1‑100 (MQ=255)
                                                                             gCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCa               >  1:1505560/1‑100 (MQ=255)
                                                                             gCAGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCa               >  1:1510289/1‑100 (MQ=255)
                                                                               aGTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGa             <  2:336983/100‑1 (MQ=255)
                                                                                gTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGaa            >  1:889728/1‑100 (MQ=255)
                                                                                gTATGCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGaa            <  1:416148/100‑1 (MQ=255)
                                                                                    gCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAGAACGTGGTGATCGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAatat         >  2:315568/1‑99 (MQ=255)
                                                                                    gCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATc        >  2:636955/1‑100 (MQ=255)
                                                                                    gCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATc        <  1:104640/100‑1 (MQ=255)
                                                                                    gCTGGTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATc        <  1:1398398/100‑1 (MQ=255)
                                                                                       ggTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGc     >  1:1379447/1‑100 (MQ=255)
                                                                                       ggTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGc     >  1:1094758/1‑100 (MQ=255)
                                                                                       ggTCCAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGc     >  1:241252/1‑100 (MQ=255)
                                                                                        gTCCAGTTCGCTTACGTCAGGGGAAAAAAAAGAAAAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCa    <  2:97729/100‑1 (MQ=255)
                                                                                          ccAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAAc  >  2:502937/1‑100 (MQ=255)
                                                                                          ccAGTTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAAc  >  2:795196/1‑100 (MQ=255)
                                                                                             |                                                                                                
AGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAAC  >  NC_000913/3524320‑3524509

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