Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,927,949	A→C	R297R (CGT→CGG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,927,949	0	A	C	100.0%	178.1 / NA	51	R297R (CGT→CGG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/30);  total (21/30)

CTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927802‑1928091                                                                                                                                                    |                                                                                                                                               cTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCg                                                                                                                                              <  1:569820/150‑1 (MQ=255)       gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATaa                                                                                                                                        <  1:146184/150‑1 (MQ=255)              tGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGTCAAAg                                                                                                                                 >  4:314767/1‑150 (MQ=255)                gcgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCGGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtt                                                                                                                               >  1:485881/1‑150 (MQ=255)                  gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGtttt                                                                                                                             >  3:18095/1‑150 (MQ=255)                   cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGttttt                                                                                                                            <  1:143211/150‑1 (MQ=255)                     aaCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGc                                                                                                                          <  2:135615/150‑1 (MQ=255)                       cTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCg                                                                                                                        >  1:69548/1‑150 (MQ=255)                         ggTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTg                                                                                                                      >  2:398173/1‑150 (MQ=255)                           ttAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGc                                                                                                                    >  2:111716/1‑150 (MQ=255)                           ttAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGc                                                                                                                    <  2:467847/150‑1 (MQ=255)                                ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                                                               <  2:491097/150‑1 (MQ=255)                                   cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                                                            <  1:963/150‑1 (MQ=255)                                            tctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaa                                                                                                   >  1:394776/1‑150 (MQ=255)                                                  accaccAGCCAGTCGGTAAATAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                             <  4:130483/150‑1 (MQ=255)                                                  accaccAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                             >  2:3235/1‑150 (MQ=255)                                                    cgccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCACCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           <  2:423876/148‑1 (MQ=255)                                                    caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           >  3:248571/1‑150 (MQ=255)                                                    caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                                                           >  1:16864/1‑150 (MQ=255)                                                      ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGGATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                                                         >  3:93412/1‑150 (MQ=255)                                                       cagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                                                        <  1:90746/150‑1 (MQ=255)                                                          ccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTa                                                                                     >  2:159005/1‑150 (MQ=255)                                                            agTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGt                                                                                   <  1:153840/150‑1 (MQ=255)                                                                  tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGAAAAAAAAAAAAAAACGATGCTGGTAGTGATCaa                                                                             <  2:205017/150‑1 (MQ=255)                                                                   aaaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACGGTTCCCACTGTTGGTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                                                            <  2:244503/150‑1 (MQ=255)                                                                    aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                                                           >  2:517785/1‑150 (MQ=255)                                                                    aaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                                                           <  4:12609/150‑1 (MQ=255)                                                                     aCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGCTGCGGGTAGTGATCAAGGc                                                                          >  2:186600/1‑150 (MQ=255)                                                                     aCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGAAGATAAAAACGATGCTGGTAGTGATCAAGGc                                                                          <  2:317688/150‑1 (MQ=255)                                                                        gCGTAAACCCTTCCGGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGt                                                                       <  2:251825/150‑1 (MQ=255)                                                                            aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                                                   <  1:262232/150‑1 (MQ=255)                                                                                       gCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTCGAACGCAGACTAGCCCGACGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                                                        <  4:16945/150‑1 (MQ=255)                                                                                             tgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg                                                  <  3:239431/150‑1 (MQ=255)                                                                                               ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGGGATCTAGGCTGCATTCGTGATCTTTgcggcgcggc                                                >  1:281830/1‑150 (MQ=255)                                                                                                  tCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGa                                             <  2:340038/150‑1 (MQ=255)                                                                                                     cgTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaa                                          <  3:161659/150‑1 (MQ=255)                                                                                                       tGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                        <  3:325393/150‑1 (MQ=255)                                                                                                          cGGAATTAACTCTTCCAACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaa                                     >  3:342127/1‑150 (MQ=255)                                                                                                           ggAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaacaaac                                    <  1:398173/150‑1 (MQ=255)                                                                                                                 aaCTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCg                              <  1:462151/150‑1 (MQ=255)                                                                                                                  aCTCTTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCgg                             >  1:430056/1‑150 (MQ=255)                                                                                                                   ctctTCCCACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGc                            >  1:121276/1‑150 (MQ=255)                                                                                                                         ccACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCACAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGc                      <  4:341582/150‑1 (MQ=255)                                                                                                                          cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCa                                                                                                                                    >  2:110209/1‑39 (MQ=255)                                                                                                                          cACTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCa                                                                                                                                    <  1:110209/39‑1 (MQ=255)                                                                                                                            cTGTTGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCAt                   >  1:471421/1‑150 (MQ=255)                                                                                                                               ttGCTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTc                >  4:160402/1‑150 (MQ=255)                                                                                                                                  cTCATCACGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc             <  1:111716/150‑1 (MQ=255)                                                                                                                                    catcaCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGt           <  2:69548/150‑1 (MQ=255)                                                                                                                                        aCGCATACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcag       <  1:516973/150‑1 (MQ=255)                                                                                                                                             tACGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGc  <  2:16864/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               CTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAGGC  >  NC_000913_3_bme_pgi/1927802‑1928091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑