Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,881	A→G	intergenic (‑271/+189)	cyoA ← / ← ampG	cytochrome o ubiquinol oxidase subunit II/muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	451,881	0	A	G	100.0%	159.6 / NA	42	intergenic (‑271/+189)	cyoA/ampG	cytochrome o ubiquinol oxidase subunit II/muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

TGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATTACTACAGCTAAATAATATTTACAGATTACGTCAGATGCGTTTTTCGTAGCGCCAGATAA  >  NC_000913/451653‑452103                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               tgtgGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtg                                                                                                                                                                                                                      >  2:333736/1‑239 (MQ=255)                               tGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATc                                                                                                                                                                                        >  2:532181/1‑239 (MQ=255)                                                    gTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTgg                                                                                                                                                                                                                         <  1:243913/185‑1 (MQ=255)                                                    gTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTgg                                                                                                                                                                                                                         >  2:243913/1‑185 (MQ=255)                                                              tttGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                                                                                                                                                         >  1:117362/1‑239 (MQ=255)                                                                      ccAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAtt                                                                                                                                                 >  2:435860/1‑239 (MQ=255)                                                                           tCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTcc                                                                                                                                              <  1:435860/237‑1 (MQ=255)                                                                                     tGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtga                                                                                                                                                                                                                     >  1:215792/1‑156 (MQ=255)                                                                                     tGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtga                                                                                                                                                                                                                     <  2:215792/156‑1 (MQ=255)                                                                                     tGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGcc                                                                                                                                  >  1:361031/1‑239 (MQ=255)                                                                                                    gTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACaaa                                                                                                                   >  1:241225/1‑239 (MQ=255)                                                                                                         ttttAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATta                                                                                                              <  1:532181/239‑1 (MQ=255)                                                                                                              agaagCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAg                                                                                                         <  2:241225/239‑1 (MQ=255)                                                                                                                      ttGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGaa                                                                                                                                                                                          <  2:616098/150‑1 (MQ=255)                                                                                                                      ttGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGaa                                                                                                                                                                                          >  1:616098/1‑150 (MQ=255)                                                                                                                        gCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGAt                                                                                               <  1:239614/239‑1 (MQ=255)                                                                                                                        gCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGAt                                                                                               <  1:16509/239‑1 (MQ=255)                                                                                                                                gCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGaa                                                                                                                                                                                          >  1:692151/1‑140 (MQ=255)                                                                                                                                gCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGaa                                                                                                                                                                                          <  2:692151/140‑1 (MQ=255)                                                                                                                                gCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGATGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGTCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGtttattt                                                                                       >  2:424925/1‑239 (MQ=255)                                                                                                                                              ggggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAg                                                                                                         <  1:111562/207‑1 (MQ=255)                                                                                                                                              ggggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAg                                                                                                         >  2:111562/1‑207 (MQ=255)                                                                                                                                                   gTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACg                                                                       >  2:38790/1‑236 (MQ=255)                                                                                                                                                   gTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACg                                                                       <  1:38790/236‑1 (MQ=255)                                                                                                                                                     gCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGa                                                                   <  2:117362/238‑1 (MQ=255)                                                                                                                                                     gCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGa                                                                   <  2:361031/238‑1 (MQ=255)                                                                                                                                                       tGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTg                                                                <  2:58743/239‑1 (MQ=255)                                                                                                                                                                cTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGATGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATtacta                                                         <  1:424925/237‑1 (MQ=255)                                                                                                                                                                       ggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAg                                                                                                         >  1:301892/1‑182 (MQ=255)                                                                                                                                                                       ggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAg                                                                                                         <  2:301892/182‑1 (MQ=255)                                                                                                                                                                              ggTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATTACTACAGCTAAATAATAttt                                         <  1:333736/239‑1 (MQ=255)                                                                                                                                                                                    aCCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGa                                                                   <  2:519750/207‑1 (MQ=255)                                                                                                                                                                                    aCCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGa                                                                   >  1:519750/1‑207 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATgag                                                                                                                                                                                                                 >  1:511212/1‑60 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATgag                                                                                                                                                                                                                 <  2:511212/60‑1 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATta                                                                                                              >  2:123358/1‑159 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATta                                                                                                              <  1:123358/159‑1 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGtttat                                                                                         >  1:28947/1‑180 (MQ=255)                                                                                                                                                                                         tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGtttat                                                                                         <  2:28947/180‑1 (MQ=255)                                                                                                                                                                                                                     aatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATTACTACAGCTAAATAATATTTACAGATTACGTCAGATGCGTTTTTCGTAGCGCCAGATaa  >  1:728800/1‑239 (MQ=255)                                                                                                                                                                                                                               agagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATtacta                                                         >  1:716472/1‑174 (MQ=255)                                                                                                                                                                                                                               agagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATtacta                                                         <  2:716472/174‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               TGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAAATTCCTGTTACTTGGATTTGGCTGCCAGGTCACATATGTTGCAAAACATATCAGGTTGATTTTAAGAAGCTGTTGCAAAAGGGCGTTATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCGCACAAAAGAACAAATATTAATAAGCGATACTGATGTTTATTTAGCCGTAATAATTACGGCGAGTGATTACTACAGCTAAATAATATTTACAGATTACGTCAGATGCGTTTTTCGTAGCGCCAGATAA  >  NC_000913/451653‑452103

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