Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,645,013	C→T	*771* (TAG→TAA)	pbpC ←	penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,645,013	0	C	T	100.0%	131.7 / NA	41	*771* (TAG→TAA)	pbpC	penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/18);  total (23/18)

TGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644886‑2645148                                                                                                                                |                                                                                                                                        tgttgttgtCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa                                                                                                                                  >  2:455541/1‑135 (MQ=255) tgttgttgtCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa                                                                                                                                  <  1:455541/135‑1 (MQ=255)        gtCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTc                                                                                                                       >  1:946640/1‑139 (MQ=255)           tGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTAACt                                                                                                                    >  1:85860/1‑139 (MQ=255)                   cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt                                                                                                            >  1:970176/1‑139 (MQ=255)                    gTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAtt                                                                                                           <  2:85860/139‑1 (MQ=255)                        gagaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  2:970176/139‑1 (MQ=255)                             ggCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  >  2:811290/1‑139 (MQ=255)                               cgcgGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTc                                                                                                >  1:437831/1‑139 (MQ=255)                                        taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCAc                                                                                       <  2:946640/139‑1 (MQ=255)                                              gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCAt                                                                                                                                 >  2:871394/1‑91 (MQ=255)                                              gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCAt                                                                                                                                 <  1:871394/91‑1 (MQ=255)                                              gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       >  1:515943/1‑117 (MQ=255)                                              gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  2:515943/117‑1 (MQ=255)                                               cAACATGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCACGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  1:419228/139‑1 (MQ=255)                                               cAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  2:437831/139‑1 (MQ=255)                                                       ccGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  2:290017/131‑1 (MQ=255)                                                       ccGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                >  1:290017/1‑131 (MQ=255)                                                               tATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTAt                                                                >  2:403643/1‑139 (MQ=255)                                                                 ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATcc                                                              >  2:792835/1‑139 (MQ=255)                                                                 ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATcc                                                              <  2:740983/139‑1 (MQ=255)                                                                 ttGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATcc                                                              >  1:740983/1‑139 (MQ=255)                                                                    atgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTc                                                           >  1:710357/1‑139 (MQ=255)                                                                       aaGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCaaa                                                        >  2:990584/1‑139 (MQ=255)                                                                         ggTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCa                                                                                                                      >  2:714139/1‑75 (MQ=255)                                                                         ggTAACATTTTTTTAAAATAAAATGAGTTTCAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCa                                                                                                                      <  1:714139/75‑1 (MQ=255)                                                                           tAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCa                                                    <  2:710357/139‑1 (MQ=255)                                                                            aaCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAg                                                   >  2:504582/1‑139 (MQ=255)                                                                               aTTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGt                                                <  1:811290/139‑1 (MQ=255)                                                                                              aatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCcgcg                                 <  1:504582/139‑1 (MQ=255)                                                                                                       tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTaa                        >  2:745631/1‑139 (MQ=255)                                                                                                            acaaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTAt                                                                >  1:23289/1‑94 (MQ=255)                                                                                                            acaaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTAt                                                                <  2:23289/94‑1 (MQ=255)                                                                                                             caaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTc                  >  1:619857/1‑139 (MQ=255)                                                                                                              aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  <  2:375540/58‑1 (MQ=255)                                                                                                              aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  >  1:375540/1‑58 (MQ=255)                                                                                                              aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCg                 >  2:52992/1‑139 (MQ=255)                                                                                                                       cAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGa        >  2:809540/1‑139 (MQ=255)                                                                                                                         aaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAAt                                                                         >  1:924077/1‑72 (MQ=255)                                                                                                                         aaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAAt                                                                         <  2:924077/72‑1 (MQ=255)                                                                                                                             aGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAaccac  <  1:403643/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                        TGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644886‑2645148

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