Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,601 T→C I401V (ATC→GTC)  mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606649,6010TC97.8% 182.2 / ‑3.8 46I401V (ATC→GTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG  >  REL606/649529‑649673
                                                                        |                                                                        
gTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGAcccc                                                                       <  2:1280535/76‑1 (MQ=255)
  tcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCa                                                                     >  1:1553111/1‑76 (MQ=255)
   cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCat                                                                    <  1:1093402/76‑1 (MQ=255)
   cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCat                                                                    <  2:2137844/76‑1 (MQ=255)
      ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATc                                                                 <  2:1721604/76‑1 (MQ=255)
      ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATc                                                                 >  1:461254/1‑76 (MQ=255)
       cGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCg                                                                <  1:492810/76‑1 (MQ=255)
         ccAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTa                                                              >  1:543269/1‑76 (MQ=255)
          cAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAg                                                             >  2:61194/1‑76 (MQ=255)
          cAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAg                                                             <  2:2126707/76‑1 (MQ=255)
           aGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAgg                                                            >  2:596623/1‑76 (MQ=255)
             gTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGcc                                                          <  1:1963246/76‑1 (MQ=255)
              tCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTATGccg                                                         >  1:2074395/1‑75 (MQ=255)
              tCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCa                                                         <  2:1006094/76‑1 (MQ=255)
               cGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCAc                                                        >  2:1947342/1‑76 (MQ=255)
                  tACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTg                                                     >  1:1486258/1‑76 (MQ=255)
                   aCCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGa                                                    <  1:1930659/76‑1 (MQ=255)
                   aCCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGa                                                    >  1:248587/1‑76 (MQ=255)
                    ccGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGAt                                                   <  2:876909/76‑1 (MQ=255)
                       gtgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATaga                                                >  1:2107312/1‑76 (MQ=255)
                       gtgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATaga                                                <  1:197071/76‑1 (MQ=255)
                          tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagacga                                             >  1:455899/1‑76 (MQ=255)
                          tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagaaga                                             >  1:1067419/1‑76 (MQ=255)
                           aaTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAgaagaa                                            >  2:1549007/1‑76 (MQ=255)
                           aaTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAgaagaa                                            >  1:965082/1‑76 (MQ=255)
                            aTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCCGATAGaagaag                                           >  2:729859/1‑76 (MQ=255)
                              gACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGt                                         <  2:1375493/76‑1 (MQ=255)
                                ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt                                       <  1:1728701/76‑1 (MQ=255)
                                ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt                                       <  2:374279/76‑1 (MQ=255)
                                ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt                                       <  1:1539275/76‑1 (MQ=255)
                                   tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg                                    <  2:764102/76‑1 (MQ=255)
                                   tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg                                    <  2:1198512/76‑1 (MQ=255)
                                              aCCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTcc                         >  2:854334/1‑76 (MQ=255)
                                                 catccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGc                      <  2:1426044/76‑1 (MQ=255)
                                                  atccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCg                     <  1:578985/76‑1 (MQ=255)
                                                     catTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATc                  >  2:2184407/1‑76 (MQ=255)
                                                      atTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGACGGTATCCGCAGATTCTTCCAGCGATCt                 >  2:1196926/1‑76 (MQ=255)
                                                      atTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCt                 >  2:163612/1‑76 (MQ=255)
                                                       ttCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCtt                <  2:1825796/76‑1 (MQ=255)
                                                        tCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtg               <  1:1628642/76‑1 (MQ=255)
                                                        tCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtg               <  1:1744246/76‑1 (MQ=255)
                                                         cGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgt              >  1:550812/1‑76 (MQ=255)
                                                            agagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGAc           >  2:605599/1‑76 (MQ=255)
                                                             gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa          <  1:118647/76‑1 (MQ=255)
                                                             gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa          <  2:1775602/76‑1 (MQ=255)
                                                             gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa          <  1:2026177/76‑1 (MQ=255)
                                                              agGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACAt         >  2:953869/1‑76 (MQ=255)
                                                                 cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa      >  2:1531067/1‑76 (MQ=255)
                                                                 cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa      >  1:719690/1‑76 (MQ=255)
                                                                 cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa      <  1:256874/76‑1 (MQ=255)
                                                                  gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg     <  1:1432489/76‑1 (MQ=255)
                                                                  gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg     >  1:1342631/1‑76 (MQ=255)
                                                                  gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg     >  2:850706/1‑76 (MQ=255)
                                                                   atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa    <  2:1938707/76‑1 (MQ=255)
                                                                   atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa    <  2:1389325/76‑1 (MQ=255)
                                                                     cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg  <  2:1749323/76‑1 (MQ=255)
                                                                     cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg  >  2:1600299/1‑76 (MQ=255)
                                                                     cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg  <  1:1631208/76‑1 (MQ=255)
                                                                     cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg  >  1:2057728/1‑76 (MQ=255)
                                                                        |                                                                        
GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG  >  REL606/649529‑649673

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