Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,459,063	A→G	T369A (ACC→GCC)	pmbA →	metalloprotease PmbA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,459,063	0	A	G	100.0%	181.1 / NA	51	T369A (ACC→GCC)	pmbA	metalloprotease PmbA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/28);  total (23/28)

TGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCACCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCAGGTA  >  NC_000913/4458929‑4459198                                                                                                                                       |                                                                                                                                        tGACTCAGTGGTTGCTGACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCggg                                                                                                                                    <  2:935776/140‑1 (MQ=255)       aGTGGCTGCTGACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCtggtg                                                                                                                              >  2:536529/1‑140 (MQ=255)             tgctgACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACg                                                                                                                        >  1:733927/1‑140 (MQ=255)               ctgACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGa                                                                                                                      >  2:122901/1‑140 (MQ=255)               ctgACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGa                                                                                                                      <  1:265662/140‑1 (MQ=255)                 gACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAAt                                                                                                                    >  2:29374/1‑140 (MQ=255)                 gACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGACGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAAt                                                                                                                    >  1:699966/1‑140 (MQ=255)                     aGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGAt                                                                                                                <  2:813566/140‑1 (MQ=255)                     aGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGAt                                                                                                                <  2:311148/140‑1 (MQ=255)                         aCTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGc                                                                                                            >  1:162247/1‑140 (MQ=255)                         aCTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGc                                                                                                            <  1:536529/140‑1 (MQ=255)                               cgcgGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGc                                                                                                      >  2:842199/1‑140 (MQ=255)                                  ggAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTg                                                                                                   <  2:720140/140‑1 (MQ=255)                                  ggAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTg                                                                                                   <  1:734820/140‑1 (MQ=255)                                     aaCTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGt                                                                                                <  2:263404/140‑1 (MQ=255)                                              tGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTAcc                                                                                       <  1:842199/140‑1 (MQ=255)                                               gAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCg                                                                                      <  1:122901/140‑1 (MQ=255)                                                                 gggcggTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCgtggtg                                                                    <  2:162247/140‑1 (MQ=255)                                                                      gTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCg                                                               <  1:29374/140‑1 (MQ=255)                                                                      gTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCg                                                               <  1:816375/140‑1 (MQ=255)                                                                        aTTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCggg                                                             <  1:247639/140‑1 (MQ=255)                                                                        aTTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCggg                                                             <  1:942550/140‑1 (MQ=255)                                                                           cacaACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCtt                                                          <  2:733927/140‑1 (MQ=255)                                                                               aCTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTgg                                                      <  2:172066/140‑1 (MQ=255)                                                                                 tGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGt                                                    >  1:958273/1‑140 (MQ=255)                                                                                  ggCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTa                                                   <  1:336968/140‑1 (MQ=255)                                                                                  ggCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTa                                                   <  2:568954/140‑1 (MQ=255)                                                                                   gCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTag                                                  >  2:874405/1‑140 (MQ=255)                                                                                        ttGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAAt                                             >  1:529131/1‑140 (MQ=255)                                                                                        ttGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAAt                                             >  2:370667/1‑140 (MQ=255)                                                                                            cGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCg                                         <  2:958273/140‑1 (MQ=255)                                                                                            cGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCCGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCg                                         <  2:18684/140‑1 (MQ=255)                                                                                              gACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGaa                                       <  2:529131/140‑1 (MQ=255)                                                                                                  aGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTc                                   <  2:699966/140‑1 (MQ=255)                                                                                                           cTTCGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCgg                          <  1:370667/140‑1 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGGGAAATTCAGTATCCGGTGa                       >  2:231636/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  2:135745/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  1:933719/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  2:73117/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  1:946644/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  1:668623/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  2:45788/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       >  2:285880/1‑140 (MQ=255)                                                                                                              cGAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGa                       <  2:316825/140‑1 (MQ=255)                                                                                                               gAGCAGATGCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAg                      >  2:173632/1‑140 (MQ=255)                                                                                                                       gCTCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCa              <  1:671235/140‑1 (MQ=255)                                                                                                                         tCAAAGAGATGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCAcc            <  1:173632/140‑1 (MQ=255)                                                                                                                               agaTGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCa      >  2:454763/1‑140 (MQ=255)                                                                                                                                 aTGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCAgg    <  2:668623/140‑1 (MQ=255)                                                                                                                                 aTGGGCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCAgg    >  1:532642/1‑140 (MQ=255)                                                                                                                                   gggCGCCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCAGGTa  >  2:717017/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                        TGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGAAACTGGGGCTGAAAAGCACCGGACATGCGGGCGGTATTCACAACTGGCGGATTGCCGGACAAGGTCTAAGCTTCGAGCAGATGCTCAAAGAGATGGGCACCGGGCTGGTGGTGACGGAATTGATGGGCCAGGGCGTGAGTGCCATTACCGGTGATTATTCCCGTGGTGCAGCGGGCTTCTGGGTAGAGAATGGCGAAATTCAGTATCCGGTGAGCGAAATCACCATCGCAGGTA  >  NC_000913/4458929‑4459198

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