Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	647,840	A→G	F353S (TTC→TCC)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	647,840	0	A	G	100.0%	142.9 / NA	37	F353S (TTC→TCC)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/16);  total (21/16)

TACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACC  >  REL606/647768‑647908                                                                         |                                                                     tACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGAcc                                                                   <  2:565809/76‑1 (MQ=255) tACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGAcc                                                                   >  1:657176/1‑76 (MQ=255) tACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGAcc                                                                   >  1:728509/1‑76 (MQ=255)   cGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCa                                                                 <  2:55235/76‑1 (MQ=255)     cATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGc                                                               >  1:1062697/1‑76 (MQ=255)        tGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCAt                                                            <  1:870479/76‑1 (MQ=255)         gcgcACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATc                                                           <  2:1221986/76‑1 (MQ=255)            cACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGc                                                        <  2:1017910/76‑1 (MQ=255)            cACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGc                                                        >  2:1222226/1‑76 (MQ=255)            cACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGc                                                        >  2:94467/1‑76 (MQ=255)             aCCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGca                                                       <  2:1112627/76‑1 (MQ=255)                 ctTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAAt                                                   >  1:744758/1‑76 (MQ=255)                 ctTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAAt                                                   >  1:744763/1‑76 (MQ=255)                    acacGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAg                                                <  1:955287/76‑1 (MQ=255)                     cacGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAgc                                               >  2:318188/1‑76 (MQ=255)                      acGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAgcg                                              >  2:65574/1‑76 (MQ=255)                        gCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGcc                                            <  1:799674/76‑1 (MQ=255)                         cTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCg                                           >  1:1157834/1‑76 (MQ=255)                          ttttCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGa                                          >  1:110061/1‑76 (MQ=255)                              cGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGAtcct                                      <  1:853025/76‑1 (MQ=255)                                acaacaTTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGAtcctcc                                    <  2:975313/76‑1 (MQ=255)                                acaacaTTTTCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGAtcctcc                                    >  2:127150/1‑76 (MQ=255)                                       tttCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTg                             >  2:644474/1‑76 (MQ=255)                                        ttCCTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCCATCCCTCATAAATGa                            >  1:315675/1‑76 (MQ=255)                                          ccTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGatc                          >  2:340091/1‑74 (MQ=255)                                          ccTGTGGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGAcc                          >  2:340100/1‑76 (MQ=255)                                             gtgGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGt                       >  2:259918/1‑76 (MQ=255)                                              tgGGTGTGGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTg                      >  2:657339/1‑76 (MQ=255)                                                  tgtgGATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAg                  <  1:956655/76‑1 (MQ=255)                                                      gATTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGa              <  2:896716/76‑1 (MQ=255)                                                       aTTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGaa             >  1:651265/1‑76 (MQ=255)                                                       aTTGACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGaa             <  2:1192151/76‑1 (MQ=255)                                                          gACATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAAccc          >  1:229227/1‑76 (MQ=255)                                                            cATTACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCt        >  2:1007197/1‑76 (MQ=255)                                                               tACAATCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACa     <  1:165206/76‑1 (MQ=255)                                                                  aaTCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAAcc  <  2:845502/76‑1 (MQ=255)                                                                  aaTCCCGGACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAAcc  <  1:828110/76‑1 (MQ=255)                                                                         |                                                                     TACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACC  >  REL606/647768‑647908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑