Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,060,522	C→T	P26S (CCG→TCG)	yihL →	putative DNA‑binding transcriptional regulator YihL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,522	0	C	T	100.0%	123.7 / NA	36	P26S (CCG→TCG)	yihL	putative DNA‑binding transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/16);  total (20/16)

AGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGCCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTAC  >  NC_000913/4060307‑4060763                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTggg                                                                                                                                                                                                                 >  1:184530/1‑250 (MQ=255)         cattATAATTCTGTGGTTACGTCAATCGAGAGCTTTGTTACAATTATCTGCAAATTCATACGCCGTTAATTGCTTTCTTTTTTTGCTTAAGTGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGTAGAATCAATCCACCGTAGAAAAGGCAAAAGAGAAACTGGATCGGTGGTTGAAAGACGGCATCCCCACGTCGGGTGTAAAACCTCCTTCAGAAAGAGAGCTGGGCGAActgct                                                                                                                                                                                                        <  2:213509/249‑1 (MQ=255)            tataAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAActgctg                                                                                                                                                                                                       <  2:18668/249‑1 (MQ=255)                aGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCAtt                                                                                                                                                                                                 >  2:130264/1‑251 (MQ=255)                           aCGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGa                                                                                                                                                                                       <  2:396631/250‑1 (MQ=255)                                    gagagCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGcc                                                                                                                                                                              >  1:379544/1‑250 (MQ=255)                                             gTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCgttgtt                                                                                                                                                                    <  2:184530/251‑1 (MQ=255)                                             gTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCgttgtt                                                                                                                                                                    <  1:625735/251‑1 (MQ=255)                                               canaATTATCCGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCGGGGCATTAAACGTATGACGCTGCGCCAGGCCttgttg                                                                                                                                                                   <  2:602420/250‑1 (MQ=255)                                                               gTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGactc                                                                                                                                                  >  2:249815/1‑248 (MQ=255)                                                                         gTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTc                                                                                                                                        <  1:29113/251‑1 (MQ=255)                                                                             aTTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTa                                                                                                                                    <  1:284135/251‑1 (MQ=255)                                                                                 ctttctttTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGa                                                                                                                                                                                                                             <  2:542701/158‑1 (MQ=255)                                                                                 ctttctttTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGa                                                                                                                                                                                                                             >  1:542701/1‑158 (MQ=255)                                                                                       tttttGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCg                                                                                                                             <  2:460989/248‑1 (MQ=255)                                                                                        ttttGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTaa                                                                                                                          <  2:409583/250‑1 (MQ=255)                                                                                          ttGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAgg                                                                                                                        <  1:165603/250‑1 (MQ=255)                                                                                                 aagcgtaagATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGGTCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCCCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTc                                                                                                                >  2:55721/1‑251 (MQ=255)                                                                                                            gCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGc                                                                                                      <  2:8177/250‑1 (MQ=255)                                                                                                               tCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCcgcg                                                                                                  >  1:240501/1‑251 (MQ=255)                                                                                                                 aTCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGCAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGCATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGAt                                                                                                >  1:32186/1‑251 (MQ=255)                                                                                                                     ggTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTa                                                                                             <  1:199226/250‑1 (MQ=255)                                                                                                                       tttAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGAGATTTAAtt                                                                                          >  2:395150/1‑251 (MQ=255)                                                                                                                        ttAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAAtt                                                                                          >  2:422006/1‑250 (MQ=255)                                                                                                                          aaaCCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTAc                                                                                        >  1:531584/1‑250 (MQ=255)                                                                                                                                        aaaCAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGGGGAAAACTACCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGACTCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGt                                                                           >  2:254190/1‑249 (MQ=255)                                                                                                                                                                         aTGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCg                                                                         >  1:358141/1‑218 (MQ=255)                                                                                                                                                                         aTGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAggg                                        >  2:128179/1‑251 (MQ=255)                                                                                                                                                                         aTGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGCGGAAACCTCCCTTCAGAGAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGATTTCAGCGGGCCGCAATTGAGCAAggg                                        >  1:29184/1‑251 (MQ=255)                                                                                                                                                                         aTGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCAGCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCg                                                                         <  2:358141/218‑1 (MQ=255)                                                                                                                                                                             aaaaGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCg                                      >  2:103282/1‑249 (MQ=255)                                                                                                                                                                                 gagaAACTGGATCGGTGGTTGAAAGATGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGcc                                >  1:41424/1‑251 (MQ=255)                                                                                                                                                                                               gtggtTGAAAGCTGGCATCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTCTGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGtttt                  >  2:534325/1‑250 (MQ=255)                                                                                                                                                                                                               aTCACCACGTCGGGCGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTAc  >  2:69348/1‑251 (MQ=255)                                                                                                                                                                                                                tCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTaa                                                                                                                          <  2:234228/130‑1 (MQ=255)                                                                                                                                                                                                                tCACCACGTCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAAGATCGTaa                                                                                                                          >  1:234228/1‑130 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  AGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGCCGGGTGGAAAACTCCCTTCAGAAAGAGAGCTGGGAGAACTGCTGGGCATTAAACGTATGACGCTGCGCCAGGCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTAC  >  NC_000913/4060307‑4060763

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