Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,672	G→A	intergenic (‑1619/‑925)	topAI ← / → ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,535,672	0	G	A	100.0%	51.5 / NA	19	intergenic (‑1619/‑925)	topAI/ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/13);  total (6/13)

GTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGAGCAGCATGATGAATTTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGGGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAATCGTGTTCACAGCAAAATTGTGATTGAGGATGAAGAACTGCTGTGTGTCGGATC  >  NC_000913/4535461‑4535901                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      gTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGAGCAGCATGATGAATTTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTc                                                                                                                                                                                                        <  1:125284/243‑1 (MQ=255)                                                              cGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGAGCAGCATGATGAATTTTTGAACAAGACGCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGGTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCg                                                                                                                                                                                    <  1:186144/202‑1 (MQ=255)                                                                                                            ntgtgGAGCAGCATGATGAATTTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCACTGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAAt                                                                                     >  1:149685/2‑251 (MQ=255)                                                                                                                      gcaTGATGAATGTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGAGGATAATCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAAt                                                                                     <  1:120274/241‑1 (MQ=255)                                                                                                                          gatgaATTTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCTTCAGGAGAAACAACATCTTCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGc                                                                       <  1:126367/251‑1 (MQ=255)                                                                                                                           attaaTTTTTGAACAAGACTCTTGCCGGAGCGCAAAAGAAAATAACGATTAGTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGTCTTTGGCCCGTTCCCGGGGCATTGATATTAc                                                                                                                                                                                      <  1:182071/136‑1 (MQ=255)                                                                                                                           atgaATTTNTNNNCNNGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTNGTTATCCTGGCAAAANGTNGAGANAACNAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTAc                                                                                                                                                                                      >  2:182071/1‑139 (MQ=255)                                                                                                                                        cAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTCTCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCa                                                                                                                                                                                                                         >  1:141401/1‑91 (MQ=255)                                                                                                                                                  gCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGCTCCCGGGGCATTGATATTACCgt                                                                                                                                                                                   <  1:86368/119‑1 (MQ=255)                                                                                                                                                      gAGCTCAAAATAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATATCGCATGTCGATGAGGATAAGCGTCAGGAGAAACAACATCTGCTGAATTATGCTGGTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAATCGTGTTCACAGc                                           <  1:174413/251‑1 (MQ=255)                                                                                                                                                          tctaaaGAGAATAACGATTATTTCCCNTTGGTTNNCTNGGCAAAAAGTGGAGCAAACNAGNNNTNNNGNCTCAANGGNNTTGGCCGGTNCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATNTCGATGATGNNAAGCg                                                                                                                                   >  2:56244/4‑159 (MQ=255)                                                                                                                                                          tCAAAAGAAAATAACGATTATGTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCAGCAATGGCTTTGGCCCGTTCCCGGTGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCg                                                                                                                                   <  1:56244/159‑1 (MQ=255)                                                                                                                                                            aaaaGACAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCATTGGCTTTGGACCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGATGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAATAGTGTTCACAGCAAAAtt                                     >  1:227049/1‑251 (MQ=255)                                                                                                                                                                      aaCGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAAGCGTATTCACAGCAAAATTGTGATTGAg                            <  1:201153/250‑1 (MQ=255)                                                                                                                                                                                   ccttgtttttacTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGTTGACGGATAAAAACTGCAATATAGCGCATGTCTATGATGATAAGCGTCATGATAAACAACATCGGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAAGCGTGTTcaca                                             <  1:42964/209‑1 (MQ=255)                                                                                                                                                                                           aTCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATNNAAACTg                                                                                                                                                                 >  2:4434/1‑96 (MQ=255)                                                                                                                                                                                            tCCTGGCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTTCGAAATTGGTTAATCGTGTTCACAGCAAAAtt                                     <  1:90925/219‑1 (MQ=255)                                                                                                                                                                                                ggCAAAAAGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAAGCGTGTTCACAGCAAAATTGTGATTGAGGATGAAGAACTGCTGTGTGTCGGAtn  <  1:49072/250‑2 (MQ=255)                                                                                                                                                                                                    aaaaGTGGAGCAAACGAGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGTATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAatgatg                                                                                                      <  1:116246/146‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      GTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGAGCAGCATGATGAATTTTTGAACAAGACTCTTGCCGGAGCTCAAAAGAAAATAACGATTATTTCCCCTTGGTTATCCTGGCAAAAAGTGGAGCAAACGGGATTTCTGGCATCAATGGCTTTGGCCCGTTCCCGGGGCATTGATATTACCGTGGTGACGGATAAAAACTGCAATATAGCGCATGTCGATGATGATAAGCGTCAGGAGAAACAACATCTGCTGAATGATGCTGTTGAAAAGCTCAATAAGATGGGTATCGCTACGAAATTGGTTAATCGTGTTCACAGCAAAATTGTGATTGAGGATGAAGAACTGCTGTGTGTCGGATC  >  NC_000913/4535461‑4535901

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