Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,392,890	G→A	*300* (TAG→TAA)	pgrR →	murein peptide degradation regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,392,890	0	G	A	100.0%	147.6 / NA	46	*300* (TAG→TAA)	pgrR	murein peptide degradation regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/22);  total (24/22)

GATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAA  >  NC_000913/1392764‑1393020                                                                                                                               |                                                                                                                                   gATACGTGTTTTGGATAAATTCACA‑CAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGa                                                                                                                       <  2:168132/139‑1 (MQ=255)  aTACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGa                                                                                                                       <  1:717699/139‑1 (MQ=255)                             gATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTtctc                                                                                                  >  2:634531/1‑133 (MQ=255)                             gATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTtctc                                                                                                  <  1:634531/133‑1 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCGTTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  <  1:235148/96‑1 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCGTTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  >  2:235148/1‑96 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                       >  1:437303/1‑107 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                       <  2:437303/107‑1 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                <  2:806129/114‑1 (MQ=255)                                                  tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                >  1:806129/1‑114 (MQ=255)                                                       tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                >  2:800106/1‑109 (MQ=255)                                                       tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                <  1:800106/109‑1 (MQ=255)                                                       tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAAc                                                                                             >  1:153133/1‑112 (MQ=255)                                                       tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAAc                                                                                             <  2:153133/112‑1 (MQ=255)                                                       tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTg                                                                  >  1:815429/1‑139 (MQ=255)                                                         cacaCCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  <  2:231510/89‑1 (MQ=255)                                                         cacaCCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  >  1:231510/1‑89 (MQ=255)                                                            acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTg                                                                                                                        <  2:567429/80‑1 (MQ=255)                                                            acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTg                                                                                                                        >  1:567429/1‑80 (MQ=255)                                                            acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAt                                                             >  2:14831/1‑139 (MQ=255)                                                            acCGTCGACATGCTGGCTCGGCATTCTCATTA‑TTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtt                                                            >  2:663396/1‑139 (MQ=255)                                                                 cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtttttt                                                        <  2:815429/139‑1 (MQ=255)                                                                  gACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                       >  2:125965/1‑91 (MQ=255)                                                                  gACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                       <  1:125965/91‑1 (MQ=255)                                                                   aCATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  <  2:58607/79‑1 (MQ=255)                                                                   aCATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                  >  1:58607/1‑79 (MQ=255)                                                                    cATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTAtt                                                     <  1:959857/139‑1 (MQ=255)                                                                      tGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTAttct                                                   <  1:14831/139‑1 (MQ=255)                                                                          ggCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                               <  2:341025/123‑1 (MQ=255)                                                                          ggCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                               >  1:341025/1‑123 (MQ=255)                                                                           gCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTa                                              >  2:697595/1‑139 (MQ=255)                                                                                 cATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATgg                                        >  1:870830/1‑139 (MQ=255)                                                                                              tataGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGa                           >  1:435310/1‑139 (MQ=255)                                                                                              tataGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGAGACAGATGCTGa                           >  1:407249/1‑139 (MQ=255)                                                                                                     agGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGAt                    >  1:795143/1‑139 (MQ=255)                                                                                                      ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGt                                                                <  1:473460/94‑1 (MQ=255)                                                                                                      ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGt                                                                >  2:473460/1‑94 (MQ=255)                                                                                                       gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACaa                                           >  1:857056/1‑114 (MQ=255)                                                                                                       gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACaa                                           <  2:857056/114‑1 (MQ=255)                                                                                                       gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATAc                  >  2:891524/1‑139 (MQ=255)                                                                                                        cTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATAcc                 <  2:873207/139‑1 (MQ=255)                                                                                                             gTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGaaaa            <  2:435310/139‑1 (MQ=255)                                                                                                                 aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACggg        >  1:791903/1‑139 (MQ=255)                                                                                                                       gCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGaa  <  2:791903/139‑1 (MQ=255)                                                                                                                          gTTTAACACTACTTGATGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                               <  1:602062/107‑1 (MQ=255)                                                                                                                          gTTTAACACTACTTGATGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                               >  2:602062/1‑107 (MQ=255)                                                                                                                               |                                                                                                                                   GATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAA  >  NC_000913/1392764‑1393020

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