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%	187.9 / NA	58	*300* (TAG→TAA)	pgrR	murein peptide degradation regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/26);  total (32/26)

AAAAACTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGT  >  NC_000913/1392757‑1393023                                                                                                                                      |                                                                                                                                      aaaaaCTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACacta                                                                                                                                  <  2:161049/139‑1 (MQ=255)                              acacCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTtctc                                                                                                     <  1:371627/139‑1 (MQ=255)                                     aTTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaatt                                                                                                                >  1:68944/1‑121 (MQ=255)                                     aTTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaatt                                                                                                                <  2:68944/121‑1 (MQ=255)                                               ttATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaa                                                                                                                  <  2:468429/109‑1 (MQ=255)                                               ttATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaa                                                                                                                  >  1:468429/1‑109 (MQ=255)                                                  tCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCt                                                                                 >  2:470553/1‑139 (MQ=255)                                                             aTCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                          >  2:642807/1‑103 (MQ=255)                                                             aTCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                          <  1:642807/103‑1 (MQ=255)                                                              tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   <  1:657905/109‑1 (MQ=255)                                                              tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   >  2:657905/1‑109 (MQ=255)                                                              tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTg                                                                     >  1:476110/1‑139 (MQ=255)                                                                cacaCCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGt                                                                   >  1:566287/1‑139 (MQ=255)                                                                 acacCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTCATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  <  2:115078/139‑1 (MQ=255)                                                                   acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   <  1:1009018/104‑1 (MQ=255)                                                                   acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   >  2:1009018/1‑104 (MQ=255)                                                                   acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAt                                                                >  1:105152/1‑139 (MQ=255)                                                                    ccGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtt                                                               <  2:105152/139‑1 (MQ=255)                                                                      gTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACactac                                                                                                                                 <  1:1048706/71‑1 (MQ=255)                                                                      gTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACactac                                                                                                                                 >  2:1048706/1‑71 (MQ=255)                                                                       tCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACacta                                                                                                                                  <  2:437973/69‑1 (MQ=255)                                                                       tCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACacta                                                                                                                                  >  1:437973/1‑69 (MQ=255)                                                                        cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGa                                                                                                                          <  1:511664/76‑1 (MQ=255)                                                                        cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGa                                                                                                                          >  2:511664/1‑76 (MQ=255)                                                                        cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      >  2:595820/1‑112 (MQ=255)                                                                        cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      <  1:595820/112‑1 (MQ=255)                                                                                        cATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  >  1:438470/1‑116 (MQ=255)                                                                                        cATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  <  2:438470/116‑1 (MQ=255)                                                                                             tcATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACa                                                                                                                                     >  1:382702/1‑44 (MQ=255)                                                                                                ttatttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGAt                                   >  2:938885/1‑139 (MQ=255)                                                                                                 tatttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  <  2:476110/139‑1 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttctt                                                                                                        >  1:405144/1‑57 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttctt                                                                                                        <  2:405144/57‑1 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  <  1:373947/95‑1 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  >  2:373947/1‑95 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGa                                        <  1:448241/121‑1 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGa                                        >  2:448241/1‑121 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATa                      >  2:1020412/1‑139 (MQ=255)                                                                                                             ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATa                      >  2:830533/1‑139 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttctt                                                                                                        <  2:741567/56‑1 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttctt                                                                                                        >  1:741567/1‑56 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAAc                                                                                                <  2:290823/64‑1 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAAc                                                                                                >  1:290823/1‑64 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  <  1:1031661/94‑1 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTc                                                                  >  2:1031661/1‑94 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATAc                     <  1:938885/139‑1 (MQ=255)                                                                                                              gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATAc                     >  1:388949/1‑139 (MQ=255)                                                                                                                       taAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGt                                                                                              >  2:2842/1‑57 (MQ=255)                                                                                                                       taAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGt                                                                                              <  1:2842/57‑1 (MQ=255)                                                                                                                        aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  >  1:988032/1‑116 (MQ=255)                                                                                                                        aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  <  2:988032/116‑1 (MQ=255)                                                                                                                        aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACggg           >  1:534633/1‑139 (MQ=255)                                                                                                                            gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGt                                                                                              >  2:474609/1‑52 (MQ=255)                                                                                                                            gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGt                                                                                              <  1:474609/52‑1 (MQ=255)                                                                                                                               cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  <  1:790546/109‑1 (MQ=255)                                                                                                                               cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  >  2:790546/1‑109 (MQ=255)                                                                                                                               cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAc    >  1:894903/1‑139 (MQ=255)                                                                                                                                 gTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTAAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGt  >  2:790235/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                      AAAAACTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGT  >  NC_000913/1392757‑1393023

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