Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	650,773	T→C	Y10C (TAT→TGT)	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
*	REL606	650,773	0	T	C	100.0%	144.4 / NA	35	Y10C (TAT→TGT)	mrdA	transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/22);  total (13/22)

CACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGG  >  REL606/650700‑650841                                                                          |                                                                     cACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAg                                                                    >  1:1023239/1‑76 (MQ=255)  aCGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGt                                                                   <  1:1796168/76‑1 (MQ=255)        gTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGaa                                                             <  1:1000929/76‑1 (MQ=255)         tCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGaaa                                                            <  1:1590127/76‑1 (MQ=255)                 aaaaTCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCt                                                    >  1:516350/1‑76 (MQ=255)                 aaaaTCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCt                                                    >  1:594075/1‑76 (MQ=255)                    aTCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTa                                                 <  1:709528/76‑1 (MQ=255)                     cccccAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAg                                                <  2:759131/75‑1 (MQ=255)                      ccccAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGt                                               <  2:580889/76‑1 (MQ=255)                        ccAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGttt                                             <  2:1158784/76‑1 (MQ=255)                        ccAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGttt                                             <  2:667785/76‑1 (MQ=255)                          aaaaaGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCa                                           >  1:1180988/1‑76 (MQ=255)                          aaaaaGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCa                                           <  2:631181/76‑1 (MQ=255)                          aaaaaGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCa                                           <  2:647802/76‑1 (MQ=255)                             aaGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATcc                                        <  2:632266/76‑1 (MQ=255)                              aGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCg                                       <  1:1774101/76‑1 (MQ=255)                                gCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCt                                     <  1:760472/76‑1 (MQ=255)                                gCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCt                                     >  1:247421/1‑76 (MQ=255)                                    ccAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGt                                 >  1:1058236/1‑76 (MQ=255)                                     cAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGtt                                >  2:638513/1‑76 (MQ=255)                                       gcgcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGtttt                              <  2:1139486/76‑1 (MQ=255)                                       gcgcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGtttt                              <  2:910456/76‑1 (MQ=255)                                        cgcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTc                             <  2:1864145/76‑1 (MQ=255)                                        cgcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGCTTTc                             >  2:996431/1‑76 (MQ=255)                                         gcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCt                            <  2:1786188/76‑1 (MQ=255)                                         gcCCGGCGCACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCt                            >  2:72400/1‑76 (MQ=255)                                                gcACAAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCaaa                     <  1:1846044/76‑1 (MQ=255)                                                  acaAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGc                   >  1:1423158/1‑76 (MQ=255)                                                   caAACAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCt                  <  2:1040118/76‑1 (MQ=255)                                                    aaaCAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTc                 >  1:340011/1‑76 (MQ=255)                                                     aaCAGCGCGGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTcc                <  2:400694/76‑1 (MQ=255)                                                            cgGACTCTGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTAtca         <  1:1869618/76‑1 (MQ=255)                                                                   tGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACgg  >  2:338334/1‑76 (MQ=255)                                                                   tGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACgg  <  2:1057656/76‑1 (MQ=255)                                                                   tGCCGTACAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCACTATTACTCACgg  >  1:791246/1‑76 (MQ=255)                                                                          |                                                                     CACGCCGGTCAGCAGCAAAATCCCCAAAAAGGCGACCAGCGCCCGGCGCACAAACAGCGCGGACTCTGCCGTATAGTCGCGAAAAGAGTTCTGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGG  >  REL606/650700‑650841

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