Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	280,709	C→T	W9* (TGG→TGA)	yagB ←	orphan antitoxin YagB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	280,709	0	C	T	100.0%	78.2 / NA	27	W9* (TGG→TGA)	yagB	orphan antitoxin YagB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/12);  total (15/12)

GGCTTCACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCC  >  NC_000913/280466‑280938                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      ggCTTCACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGttt                                                                                                                                                                                                                                >  1:84787/1‑251 (MQ=255)      cACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGAACCTGGGATTCAGATCCCAGGCGAGCAGCATCAGCTCCAGATGAGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGCATGAGCCGGTGAATCCAGCCCGGTCAGCAAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGACACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAg                                                                                                                                                                                                                             >  1:373636/1‑249 (MQ=255)        ccAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTAAGCTCACAGGAGAGCCGCATCAGCTCCCGATGTGCGACAAAATGCGGGACGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTACGCCCTCAAGTTTGAGt                                                                                                                                                                                                                            >  1:359640/1‑248 (MQ=255)                                                                                             cGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCAAGCCCGGTCAGCCAGAAAATGCAGCCGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGACGACAATCAGGTGGTc                                                                                                                                    >  1:349941/1‑251 (MQ=255)                                                                                                      gCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTAACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCttt                                                                                                                                                                                                                   >  1:166030/1‑163 (MQ=255)                                                                                                              gCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGt                                                                                                                         <  1:228958/245‑1 (MQ=255)                                                                                                                                          cACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATc                                                                                       <  1:471751/251‑1 (MQ=255)                                                                                                                                                    gCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTc                                                                                                                                                                                                                      >  1:212785/1‑114 (MQ=255)                                                                                                                                                    gCCGGTGAATCCAACCCTGTCAGCCAGAAAATGCAGACTATGGCCTTCCTACATCAGTCGTGCTCCGAAACACGGTGTGAGATCTCGTTTAAGCCCTCAAGTTTGAGTag                                                                                                                                                                                                                          <  1:280069/110‑1 (MQ=255)                                                                                                                                                     ccGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCggn                                                                                                                                                                           <  1:58809/156‑2 (MQ=255)                                                                                                                                                      cGGTGAATCCAGCCCGGTCAGCCAGAAAAGGCAGACGATGGCCTTCCTGCACCAGTCGTGCTCCGAAACTCGGTGTGATATCTCGTTTGAGCCCTCAAGTTTGAGTAGTTTCGATATGCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTtgcctg                                                                                                                                                        <  1:479048/174‑1 (MQ=255)                                                                                                                                                                           ccAGAAAATGCAGACGATGACCTTCCTGCTCCAGTCGAGCTCCGAAACACGGTGTGATATCCCTTTTAAGCCCTCAAGTGTGAGTAGAGTCTTTATTCAtca                                                                                                                                                                                                           <  1:477853/102‑3 (MQ=255)                                                                                                                                                                                   tGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGt                                                                                                                                                                                                                            <  1:38894/77‑1 (MQ=255)                                                                                                                                                                                   tGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGt                                                                                                                                                                                                                            >  2:38894/1‑77 (MQ=255)                                                                                                                                                                                    gCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGCAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTAGGGGACATGGTGAACACTTTTTAACATCCTTTACCCCTGGTGATCGACttttt                                                                               >  1:328495/1‑217 (MQ=255)                                                                                                                                                                                     naGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAggtggt                                                                                                                                     >  1:298480/2‑162 (MQ=255)                                                                                                                                                                                       gACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTaccac                                          <  1:478473/251‑1 (MQ=255)                                                                                                                                                                                               ccTTCNNGCNNNAGNNGTGCTCCGAAACACGGTGTGATANCTCGTTTAAGCCCTCAAGTTTGAGTAGAGNCNTNATTCNNAATNTTATTCCTGTNAAAGAAGACCATGt                                                                                                                                                                                >  2:371530/1‑109 (MQ=255)                                                                                                                                                                                                    cntcNCCAGTNGTGCNCCGAAACACGGTGTGNTATCTCGTATAAGCCCTCANNTNTGAGTNGAGTc                                                                                                                                                                                                                      <  2:212785/63‑1 (MQ=255)                                                                                                                                                                                                      gCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTTCCCATGGTGATCGGCGTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGc                             <  1:75621/249‑1 (MQ=255)                                                                                                                                                                                                       cANNNNNCGTNCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAggtggt                                                                                                                                     <  2:298480/144‑1 (MQ=255)                                                                                                                                                                                                                  ntCCGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATAAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGc               >  1:95307/2‑251 (MQ=255)                                                                                                                                                                                                                   nccGAAACACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCAccccc                                                            >  1:306898/2‑205 (MQ=255)                                                                                                                                                                                                                    ccGAAACACGGAGTGATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCAt             >  1:14362/1‑251 (MQ=255)                                                                                                                                                                                                                      naaaCACGGTGTGATATCTCGTTTAAGCCCTCAAGTTTGAGAAGAGTCTTTATTCATAATTTTATTCCTGTTAACGAAGACCATGTTCGGCGAATGAGCATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCg                                                                     >  1:101817/2‑193 (MQ=255)                                                                                                                                                                                                                      gAAACACGGTGTGATCTCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACACTCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCAtct           >  1:415250/1‑251 (MQ=255)                                                                                                                                                                                                                                gtgATATCTCGTTTAAGCCCTCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGcc  <  1:201542/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      GGCTTCACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTCGGCGAATGAGAATATTTGCCTGCCGCCGACAATCAGGTGGTCGGGGACATGGTGAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCC  >  NC_000913/280466‑280938

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