Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000964	2,581,732	(T)6→7	intergenic (‑115/+39)	pstS ← / ← pbpA	phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000964	2,581,732	1	.	T	100.0%	86.0 / NA	23	intergenic (‑115/+39)	pstS/pbpA	phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (15/8);  total (15/8)

CCTTTTGCAGAATCGCTATTCGATTTTTTACTTTCTCCTGCATTTCCGCATGCAGCTGCAATCATCATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTT‑GTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATCAATGTAAGTCCAAGGAACTACGACAGAAATAGCAATCTCAGGATTGTC  >  NC_000964/2581505‑2581977                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ccTTTTGCAGAATCGCTATTCGATTTTTTACTTTCTCCTGCATTTCCGCATGCAGCTGCAATCATCATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCAcat                                                                                                                                                                                                                                  <  1:27317/250‑1 (MQ=255)            aTCGCTATTCGATTTTTTACTTTCTCCTGCATTTCCGCATGCAGCTGCAATCATCATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGtt                                                                                                                                                                                                                      >  1:35582/1‑251 (MQ=255)                                                    gCAGCTGCAATCATCATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCAtttt                                                                                                                                                                              <  1:298463/251‑1 (MQ=255)                                                             atcatcATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt                                                                                                                                                                     <  1:270984/251‑1 (MQ=255)                                                                atcatAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGt                                                                                                                                                                  <  1:105647/251‑1 (MQ=255)                                                                     aaaaGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCttt                                                                                                                                                              <  1:687/250‑1 (MQ=255)                                                                      aacgcagcCCTCAGAAGCATAAGAACAAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAATGTTCTATGTAAAAAGGGATTGAAACGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTCTCAGAAGACGATGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTAt                                                                                                                                                            >  2:25669/4‑251 (MQ=255)                                                                                   gAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTT‑AGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGt                                                                                                                                               >  2:30895/1‑250 (MQ=255)                                                                                       cATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACct                                                                                                                                           >  2:61934/1‑251 (MQ=255)                                                                                          aaGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCACTTTTGTCTTTATTCTTTTGCTGAGTACcttct                                                                                                                                        >  2:204831/1‑251 (MQ=255)                                                                                                cAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGAATAGAGTAAAGTTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCCCATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCATTTTTGTCTTTATTCTTTTGCTGCGTACCTTCTTTGTcc                                                                                                                                  >  1:120252/1‑251 (MQ=255)                                                                                                     tGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCtt                                                                                                                             >  1:43507/1‑251 (MQ=255)                                                                                                                       gAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAat                                                                                                            >  1:7233/1‑250 (MQ=255)                                                                                                                       gAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAat                                                                                                            >  1:137666/1‑250 (MQ=255)                                                                                                                            cccccTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGCAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGTTTCGATTTCAGCTCAAAATATTTAt                                                                                                      >  1:85049/1‑251 (MQ=255)                                                                                                                                  gCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTa                                                                                                <  1:261395/251‑1 (MQ=255)                                                                                                                                                               ataGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAg                                                                   <  2:338529/251‑1 (MQ=255)                                                                                                                                                                      aaaGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTATTTATCAGAAGACGTTGTGTTTTCTGCATTTTTTTCCTTTTTTTCTTTATTCTTTTGCTGGGTACTTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCATTTTCGTTTGTAATTGAGTATCttt                                                            >  2:105233/1‑251 (MQ=255)                                                                                                                                                                               aTGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAAt                                                   >  1:166818/1‑251 (MQ=255)                                                                                                                                                                                           tttGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGTTTTCGTTTGTATTTGAGTATCTTTGATTGTAATCAATGTAAGTcc                                       >  2:107514/1‑251 (MQ=255)                                                                                                                                                                                           tttGAACCGAATGTTAAGGGTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATCAATGTAAGTcc                                       >  2:324378/1‑251 (MQ=255)                                                                                                                                                                                              gAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATCAATGTAAGTCCa                                      >  1:109174/1‑249 (MQ=255)                                                                                                                                                                                                       gTTAAGGTTTTGTAAATCAGTGCGTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATCAATGTAAGTCCAAGGAACTACg                            <  2:35582/250‑1 (MQ=255)                                                                                                                                                                                                                                ttttt‑GTAAAAAAAAAACGCCCACATCATGTGAGCTTTTTTTTTTAGTTATCAGAAGACGTTTTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTTGTTTGCTGAGTTCCTTCTATGTCCTTCTTCTATTTCAGCTCAAAATATTTTTCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATGAATGTAAGTCCAACGAACTACGACAGAAATAGCAATCTCAGGATTGTc  <  2:27380/250‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CCTTTTGCAGAATCGCTATTCGATTTTTTACTTTCTCCTGCATTTCCGCATGCAGCTGCAATCATCATAAAAGCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGTTTTTT‑GTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAATATTTATCCATTACTTCCCTGCCGATTTCGTTTGTAATTGAGTATCTTTGATTGTAATCAATGTAAGTCCAAGGAACTACGACAGAAATAGCAATCTCAGGATTGTC  >  NC_000964/2581505‑2581977

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