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

ACTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAC  >  NC_000913/1392761‑1393021                                                                                                                                  |                                                                                                                                    aCTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTg                                                                                                                            <  1:1250952/139‑1 (MQ=255)                     aaTTCACACCAGATGTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGTCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                        <  2:593825/139‑1 (MQ=255)                               aGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTatc                                                                                              <  1:933998/139‑3 (MQ=255)                                          gTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                   >  1:128228/1‑107 (MQ=255)                                          gTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                   <  2:128228/107‑1 (MQ=255)                                          gTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAAGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATtaatttaatt                                                                                                         >  2:1502788/1‑117 (MQ=255)                                          gTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAAGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATtaatttaatt                                                                                                         <  1:1502788/117‑1 (MQ=255)                                                     tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTGACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGt                                                                                   >  2:1098454/1‑128 (MQ=255)                                                     tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                 >  1:248362/1‑114 (MQ=255)                                                     tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                 <  2:248362/114‑1 (MQ=255)                                                     tACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGt                                                                                   <  1:1098454/128‑1 (MQ=255)                                                          tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACa                                                                                                                                   <  1:40429/75‑1 (MQ=255)                                                          tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTg                                                                   >  1:611262/1‑139 (MQ=255)                                                          tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTg                                                                   >  1:768173/1‑139 (MQ=255)                                                                         tGCTGGCTGGGTATTCTCTTTAGATATAGATAGGCTGAAGTATAAAGGTACTGTTTAACACTACTTGGTGATACATTAATTTAATTCTTCTCTTAACATATTCTCAGTGCCTTTCAACGTTTTGGTCATTTTTTAtgct                                                    <  1:1436277/139‑1 (MQ=255)                                                                             ggCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGt                                                >  2:1061977/1‑139 (MQ=255)                                                                             ggCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGt                                                <  1:495973/139‑1 (MQ=255)                                                                                         tcATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACa                                    <  2:768173/139‑1 (MQ=255)                                                                                              atttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtttttt                                                         >  1:663614/1‑113 (MQ=255)                                                                                              atttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtttttt                                                         <  2:663614/113‑1 (MQ=255)                                                                                               tttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCt                              <  2:1178249/139‑1 (MQ=255)                                                                                                  atagatagGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                >  2:1339703/1‑134 (MQ=255)                                                                                                  atagatagGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                <  1:1339703/134‑1 (MQ=255)                                                                                                       tagGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaa                                                                                                                <  1:233022/49‑1 (MQ=255)                                                                                                       tagGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAttaa                                                                                                                >  2:233022/1‑49 (MQ=255)                                                                                                         ggCTGTAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCGTCTCTTAACGTATTCTGAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATa                    >  1:782302/1‑139 (MQ=255)                                                                                                         ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCa                                                                           <  1:344532/84‑1 (MQ=255)                                                                                                         ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCa                                                                           >  2:344532/1‑84 (MQ=255)                                                                                                         ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                <  1:312102/127‑1 (MQ=255)                                                                                                         ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                >  2:312102/1‑127 (MQ=255)                                                                                                         ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATa                    >  1:1380989/1‑139 (MQ=255)                                                                                                          gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGAc                                     >  1:1414208/1‑121 (MQ=255)                                                                                                          gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGAc                                     <  2:1414208/121‑1 (MQ=255)                                                                                                                    aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGAc                                     >  1:839173/1‑111 (MQ=255)                                                                                                                    aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGAc                                     <  2:839173/111‑1 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGt                                                                                   >  2:1331837/1‑58 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGt                                                                                   <  1:1331837/58‑1 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGGAAACGGGTTTGAAc  >  1:330908/1‑139 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAc  >  2:330842/1‑139 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAc  >  2:1378871/1‑139 (MQ=255)                                                                                                                           cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTCTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAc  >  2:878784/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                    ACTGATACGTGTTTTGGATAAATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAAC  >  NC_000913/1392761‑1393021

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