Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,073	A→G	intergenic (‑1020/‑1524)	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,073	0	A	G	100.0%	89.0 / NA	29	intergenic (‑1020/‑1524)	topAI/ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/10);  total (19/10)

AGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTAC  >  NC_000913/4534847‑4535323                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           aGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAACGCGGAATGTATCTGCATGCGCATCGTCGCTGTTTCAATATCATCATCGACTACTGTAACTCACTTTGCTATCACGGAAAGTTGCGGACTAAAAGAGGGATGGAAAAAGGACGTTTTTCCCGGCACTGGGATATTTACATATTGTTGGTAGAGGCATGAAGCACAAA‑GGAGGG‑GCCGGACAAAGCCCTCTTAAAGCTGAAAAAAt                                                                                                                                                                                                                                       >  2:143375/1‑248 (MQ=255)    ccAGTTCAACTCTCACTACCAGTATGATCCCAATTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCGATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATAGTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATagcagc                                                                                                                                                                                                                                 <  2:310243/251‑1 (MQ=255)                                  cAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCAAGGCTGGTGGGCCCTAAGGATGATTTagg                                                                                                                                                                                                   >  2:266241/1‑250 (MQ=255)                                          aGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCGATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCAtt                                                                                                                                                                                           >  1:143003/1‑251 (MQ=255)                                                               aGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTgaacca                                                                                                                                                                        >  1:191256/1‑244 (MQ=255)                                                                  aTCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCGATTACGGGGAACCACTGAACAAAgttg                                                                                                                                                                   >  2:328591/1‑251 (MQ=255)                                                                     gtcgCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGACGGAAACAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CACATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGttgttg                                                                                                                                                                >  1:35635/1‑251 (MQ=255)                                                                        gCTGTTTCAATAACATCATCGACTCCTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGCAGCTGAAACAATAGCAGCATGGCTGGTGGCCCATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGgtg                                                                                                                                                             >  2:52437/1‑251 (MQ=255)                                                                           gTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGtgttg                                                                                                                                                          <  1:308428/251‑1 (MQ=255)                                                                                                   gTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATCTTTCCATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGCAGCTGAAACAATAGCAGCATGGCTGGTGGCTCATAAGGATGATATAGAGCGCCATTACGGGGAAACACTGTACAAAGTTGTTGGTGTTGTGACGCCCTTTTCGGCGTAGGtt                                                                                                                                   >  1:233147/1‑250 (MQ=255)                                                                                                   gTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑TAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGTGAACCACTGTACAAAGTTGTTGGTGTTTTGACGCCATTTTCGGCGTAGGtt                                                                                                                                   <  1:266241/250‑1 (MQ=255)                                                                                                     aattCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGTAAAAAGGACGATTTTCCCTGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCTTAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTTGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTTGTGTTGTGACGCCATTTTCGGCGTAGGTTAAt                                                                                                                                <  1:317445/248‑1 (MQ=255)                                                                                                                       ggTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTg                                                                                                              <  2:87226/251‑1 (MQ=255)                                                                                                                          acgttgccgcttacacgagGGATGGACTACGGACGCTTTTCCCCGCACTGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CCCTTGGCGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCCCATACGGATGATTTAGTGCGCCAATAAGGGGAGCCACTGTACAAAGTTGTTGGTGTAGTGACGCCATTTTCGGCGTGGGTTAAAGATATAAAAATGTCATTgggt                                                                                                           >  2:68285/17‑248 (MQ=255)                                                                                                                                           agGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAt                                                                                                                     <  1:148277/224‑1 (MQ=255)                                                                                                                                           agGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGAGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAt                                                                                                                     >  2:148277/1‑224 (MQ=255)                                                                                                                                                    aaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGa                                                                                 <  1:328591/251‑1 (MQ=255)                                                                                                                                                    aaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAAGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTaa                                                                                                                                 >  1:155694/1‑203 (MQ=255)                                                                                                                                                    aaaaGGACGATGTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTaa                                                                                                                                 <  2:155694/203‑1 (MQ=255)                                                                                                                                                        ggACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGc                                                                               >  1:148046/1‑249 (MQ=255)                                                                                                                                                               attCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTa                                                                       <  2:375418/249‑1 (MQ=255)                                                                                                                                                                        aaTGGGATATTTACCTATTGATGGTAGAGGCATGAAGC‑CACATGGAGGGAGTCG‑TCATAGCCCTCTTGAAGCTGAAAAAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGTACCCCTGGACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTAACTGACGGTgg                                                             >  2:121863/1‑251 (MQ=255)                                                                                                                                                                            ggATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCAATTACGGGGAACCACTGTACAACGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGGGTAAGCTGGAAATTAATGGCAAAGACGTGACGGGTTTACTGACATTggggcc                                                         >  2:222504/1‑249 (MQ=255)                                                                                                                                                                              atatTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGCTATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACtg                                                       >  1:298664/1‑251 (MQ=255)                                                                                                                                                                                 tttACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                                >  2:417654/1‑191 (MQ=255)                                                                                                                                                                                 tttACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCAt                                                                                                                <  1:417654/191‑1 (MQ=255)                                                                                                                                                                                      atatTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTcc                                               >  1:340619/1‑251 (MQ=255)                                                                                                                                                                                                       aTGAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGaa                                >  2:217262/1‑249 (MQ=255)                                                                                                                                                                                                         gAAGC‑CAAATGGAGGGAGCCG‑TCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAggg                            >  2:455177/1‑251 (MQ=255)                                                                                                                                                                                                                                     aCCCTCTTGAAGCTGAAACACTAACAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCCGGTTAATGCTATCAAACTGTCATTGCGTAAGCTGGAAATTAATGGCCAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTAc  >  2:379207/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           AGCCCAGTTCAACTCTCACTACCAGTATGATCCCAAGTTCGAGCGCGGAATGTATCTGTATGAGCATCGTCGCTGTTTCAATAACATCATCGACTACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGC‑CAAATGGAGGGAGCCG‑TCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTAC  >  NC_000913/4534847‑4535323

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