Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,829 C→A V325F (GTT→TTT)  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,8290CA97.5% 119.6 / ‑5.5 40V325F (GTT→TTT) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (15/24);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAT  >  REL606/649757‑649901
                                                                        |                                                                        
tGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAaccc                                                                       <  1:1314567/76‑1 (MQ=255)
  cgcgTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCt                                                                     <  1:1357987/76‑1 (MQ=255)
  cgcgTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCt                                                                     <  1:185487/76‑1 (MQ=255)
     gTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCtgtg                                                                  >  1:148294/1‑76 (MQ=255)
         tCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAAACCCTGTGTGgc                                                              >  2:270304/1‑76 (MQ=255)
          cACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGCTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGgcg                                                             >  2:695798/1‑76 (MQ=255)
                   gcTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCAcc                                                    <  1:259393/76‑1 (MQ=255)
                   gcTAAAGGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCAcc                                                    <  1:844560/76‑1 (MQ=255)
                     tCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAg                                                  <  2:1235349/76‑1 (MQ=255)
                       aaggCCAAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCg                                                <  1:1012164/73‑1 (MQ=255)
                       aaTGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCg                                                <  2:26812/76‑1 (MQ=255)
                         tGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgt                                              >  2:29727/1‑76 (MQ=255)
                           ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                            <  2:169554/76‑1 (MQ=255)
                           ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                            >  1:8515/1‑76 (MQ=255)
                           ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                            >  1:723993/1‑76 (MQ=255)
                           ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                            >  1:406567/1‑76 (MQ=255)
                           cagAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt                                            <  2:314738/74‑1 (MQ=255)
                            cGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                                           >  1:176871/1‑76 (MQ=255)
                            cGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTa                                           >  1:1501138/1‑76 (MQ=255)
                               aaCCGTCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                        <  1:1347560/76‑1 (MQ=255)
                               aaCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc                                        >  2:785903/1‑76 (MQ=255)
                               aaCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTCCACCAGCGGTGTATTc                                        >  2:231728/1‑76 (MQ=255)
                                     cacaTAGTGTTTAACTGTAGGCGCGGGAGGATAAATCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCg                                  <  2:362282/76‑1 (MQ=255)
                                     cacaTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCg                                  >  2:554766/1‑76 (MQ=255)
                                       caTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGtt                                <  2:329237/76‑1 (MQ=255)
                                        aTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTc                               >  1:891868/1‑76 (MQ=255)
                                          aGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTcaa                             >  2:842134/1‑76 (MQ=255)
                                            ggTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTcaaca                           <  2:1014974/76‑1 (MQ=255)
                                                tAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGc                       >  1:940529/1‑76 (MQ=255)
                                                 aaCTGGAGACGCGGGAGGATAAAGCCGCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCg                      <  2:1437601/76‑1 (MQ=255)
                                                  aCTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCgg                     <  1:978167/76‑1 (MQ=255)
                                                  aCTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCgg                     <  1:696856/76‑1 (MQ=255)
                                                   cTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGa                    <  1:819675/76‑1 (MQ=255)
                                                    tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa                   <  1:966035/76‑1 (MQ=255)
                                                    tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa                   <  1:774059/76‑1 (MQ=255)
                                                    tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa                   <  1:376733/76‑1 (MQ=255)
                                                    tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa                   <  1:248901/76‑1 (MQ=255)
                                                    tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa                   <  1:1293955/76‑1 (MQ=255)
                                                     gTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaat                  <  1:774330/76‑1 (MQ=255)
                                                     gTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaat                  <  1:1326375/76‑1 (MQ=255)
                                                      tAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaata                 <  1:1473472/76‑1 (MQ=255)
                                                           gcgGGAGGATAACCCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCtt            <  2:1424112/76‑1 (MQ=255)
                                                                 ggaTAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTgg      <  1:1365706/76‑1 (MQ=255)
                                                                   aTAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGgag    <  1:583364/76‑1 (MQ=255)
                                                                     aaaaCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt  <  1:689084/76‑1 (MQ=255)
                                                                        |                                                                        
TGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAT  >  REL606/649757‑649901

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