Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 647,808 T→G K364Q (AAA→CAA)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606647,8080TG100.0% 74.2 / NA 19K364Q (AAA→CAA) mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/10);  total (9/10)

CTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAAC  >  REL606/647738‑647879
                                                                      |                                                                       
cTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCtgt                                                                    >  1:83113/1‑76 (MQ=255)
            aGATCCCGAGCCACTGCTTACGTATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTg                                                        <  1:831568/76‑1 (MQ=255)
            aGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTg                                                        <  1:714544/76‑1 (MQ=255)
                 ccGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGAGTGACAtt                                                   <  1:759544/76‑1 (MQ=255)
                  cGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTa                                                  >  1:740761/1‑76 (MQ=255)
                    aGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACa                                                <  2:645655/76‑1 (MQ=255)
                      ccACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAAt                                              >  2:805462/1‑76 (MQ=255)
                       cACTGCTTACGTATTTCTCACCTCATACACAATTTGCGACAANANTTGCCAGNGGGTGTGGATTGACATTACAATc                                             <  1:1068406/76‑1 (MQ=255)
                          tGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCg                                          <  1:1078840/76‑1 (MQ=255)
                            cTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGaa                                        >  2:102279/1‑76 (MQ=255)
                               aCGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAAccc                                     <  2:362511/76‑1 (MQ=255)
                                cGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCa                                    <  2:617565/76‑1 (MQ=255)
                                     tGCGCACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCAt                               <  1:700837/76‑1 (MQ=255)
                                        gcACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAg                            <  2:251426/76‑1 (MQ=255)
                                        gcACCTCTTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAg                            >  2:397862/1‑76 (MQ=255)
                                              ctTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAAt                      >  2:1137984/1‑76 (MQ=255)
                                              ctTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAAt                      >  2:659350/1‑76 (MQ=255)
                                              ctTACACGCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAAt                      <  1:248947/76‑1 (MQ=255)
                                                     gCTTTTCGACAACATTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGcc               >  2:490657/1‑76 (MQ=255)
                                                                  aTTTGCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAAc  >  2:755454/1‑76 (MQ=255)
                                                                      |                                                                       
CTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAAC  >  REL606/647738‑647879

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