Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,448,304	C→T	intergenic (+69/‑143)	ytfP → / → chpS	gamma‑glutamylamine cyclotransferase family protein YtfP/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,448,304	0	C	T	100.0%	87.8 / NA	30	intergenic (+69/‑143)	ytfP/chpS	gamma‑glutamylamine cyclotransferase family protein YtfP/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/19);  total (11/19)

GTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTG  >  NC_000913/4448063‑4448542                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               gTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCtt                                                                                                                                                                                                                                       <  1:140805/251‑1 (MQ=255)          aCGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGAGTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGGAATTTCCTCACCTTTTGCTTTtc                                                                                                                                                                                                                              <  2:160243/251‑1 (MQ=255)             ccctgCTGGCCGAACTTGATGCCTTGCTCACCAGGGGCGGGGAATACTCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCTCCTTCGGTTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTtctct                                                                                                                                                                                                                           <  2:172445/247‑1 (MQ=255)                         aaCTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTTCTTTTCTCTCCGAGCCGCttt                                                                                                                                                                                                               >  1:114551/1‑251 (MQ=255)                         aaCTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCttt                                                                                                                                                                                                               <  1:108606/251‑1 (MQ=255)                          aCTTGATGCCTTGCGCACCAGGGGCAGTGGATACGCGCGCCGGTTGAGTCAGACGCCGCACGGGGGAGCATGGTTGTACGTTTATCGGCGTCACGTCAATGGATTTCAGCTAATTGAAAGCGGCGACTGGATAGACCGGGATCAGAAACCGTATGCAGACGCCACCTTCGGGGGGCGTTGTTTTTTGCGAGAAGCAGCACATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTc                                                                                                                                                                                                              >  1:143228/1‑251 (MQ=255)                              gATGCCTGGCGCCCCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCat                                                                                                                                                                                                           <  2:195675/250‑1 (MQ=255)                                       cgcACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCCTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTACTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAAc                                                                                                                                                                                                 >  2:255907/1‑251 (MQ=255)                                        gcACCAGGGGCGGTGAATACGCGCGCCAGTTTATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGAGCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGTGATAAGTAACCATATGCAGACGTCACCTTCTGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACg                                                                                                                                                                                                <  1:239386/251‑1 (MQ=255)                                                        aTACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAACACTctgct                                                                                                                                                                                >  2:152951/1‑251 (MQ=255)                                                                                              cATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACGGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGGTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGGTTCCTGTa                                                                                                                                          <  1:223943/251‑1 (MQ=255)                                                                                               aTGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAg                                                                                                                                         <  1:48128/251‑1 (MQ=255)                                                                                                      tACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTg                                                                                                                                  >  1:149469/1‑251 (MQ=255)                                                                                                           ttATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCGCACCTTTTGCTGTTCTCTCCGAGCCGCTTTCCATAGCTAGTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGTGGTAAAACGTTTCCTTTAGCACCGTGAGtt                                                                                                                              <  2:241138/250‑1 (MQ=255)                                                                                                                    aCCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGTTAAGTAACCATATGCATACGCCACCTTCGGGAGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGAAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGt                                                                                                                     >  1:241925/1‑250 (MQ=255)                                                                                                                               gATTAAAGCTCATAGAAAGCGGAGACTGTTTACACAGGGATACGAAACCATATGCATACGCCCCCTTCGGGGGGCGTTGTTTTTTGCGAGACGACTCGCACTCTGTTTTGTAATTTCCTAACCTTTTGCTTTTCTCTCAGCGCCGCTTTCCATATCTATTAACGCATCTATAACTCTGCTGTCAGTCACAAATGCGCAGGAGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAagg                                                                                                          >  2:51712/1‑250 (MQ=255)                                                                                                                                     aGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATGTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAACTGCGCAGTGTGAAAACGATTCCTGTAGCACCGTGAGTTATACTTTGAAGAACTTAAGGAGGTGc                                                                                                    >  2:8469/1‑250 (MQ=255)                                                                                                                                          aTTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGc                                                                                              >  1:54122/1‑251 (MQ=255)                                                                                                                                                                                     cATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCAt                                                   <  1:60753/251‑1 (MQ=255)                                                                                                                                                                                        aCGCCACCTTCGGGTGGCGTTGTTTTTGGCGAGACGACGCGCAGTCTGTTTTGTAATTTCCTCTCCTTTTGCGTTTCTCTCCGAGCCGCTTTCCATATCTATGAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                 <  2:114551/250‑1 (MQ=255)                                                                                                                                                                                                       ggCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGaa                                 <  1:37745/251‑1 (MQ=255)                                                                                                                                                                                                             gttTTTTGAGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaact                           >  2:221387/1‑251 (MQ=255)                                                                                                                                                                                                               tttttGCGAGACGACTCACATTCTTTTTTGTAATTTCCTCACCTTTGGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCAGAAAAAACTCTGCTGGCATTCACAAATGCGCAGTGGTAAAACGTTGCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGAGGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCAGTCCCAATATCGGAATGAATGaactt                          <  2:38512/250‑1 (MQ=255)                                                                                                                                                                                                                ttttGCGAGACGACTCGCATTCTGTTTGTTAATTTCCTCCCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTACTACCATAAAAAGATGGGGGAACAGTGCAGGTATGTTCATTCCCAATATCGTAATGGAAGaactta                         <  2:249725/250‑1 (MQ=255)                                                                                                                                                                                                                    gCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAATACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAacttaactta                    <  1:152951/251‑1 (MQ=255)                                                                                                                                                                                                                         acgacTCGCATTCTGTTTTGTAATTTCCTCGCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATAGCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGGATAACTTAAGGAGGTGCAGATGCGGATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGcc               <  1:123007/251‑1 (MQ=255)                                                                                                                                                                                                                         acgacTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGcc               <  2:149469/251‑1 (MQ=255)                                                                                                                                                                                                                               cGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTta                                                                                                                             >  1:69315/1‑135 (MQ=255)                                                                                                                                                                                                                               cGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTta                                                                                                                             <  2:69315/135‑1 (MQ=255)                                                                                                                                                                                                                                       gTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTg  <  1:178792/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               GTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTG  >  NC_000913/4448063‑4448542

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