Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,616,535	A→G	D177D (GAT→GAC)	yjjI ←	DUF3029 domain‑containing protein YjjI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,616,535	0	A	G	100.0%	103.9 / NA	33	D177D (GAT→GAC)	yjjI	DUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

GCTCGGCAATGGCTTTCAGGTTAAGGCGTACCAGCGTG‑CTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGT  >  NC_000913/4616288‑4616772                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              gcTCGGCAATGGCTTGCACGTTAAGGCGTACCA‑TGTGGCTGCCACCCCCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCGTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAAGGGTGATTTCAGGATCGTAGATAAAGGCCAGGTTCGGGGAAACCTGCTTCATCtctgcgtct                                                                                                                                                                                                                                            <  2:47183/251‑1 (MQ=255)                    gTTAAGGCGTACCAGCGTG‑CTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCggg                                                                                                                                                                                                                         <  1:22631/251‑1 (MQ=255)                                    gTG‑CTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                                                                                                                                                                                                            >  1:235474/1‑248 (MQ=255)                                          gcCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACg                                                                                                                                                                                                       <  1:236525/248‑1 (MQ=255)                                                           gCAGTGAGTTGTAACAGCTCACAATCCCGTGGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATa                                                                                                                                                                                   <  2:16237/251‑1 (MQ=255)                                                                  gTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGt                                                                                                                                                                            <  2:60197/251‑1 (MQ=255)                                                                       aaCAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCTCAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGc                                                                                                                                                                       <  2:56397/251‑1 (MQ=255)                                                                           gCTCACAATCCCGTAGCCCCATTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTAATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGCCTACACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATCTTGGCGTGCATAACGGCGTATGGCAggg                                                                                                                                                                   >  2:102444/1‑251 (MQ=255)                                                                                         aGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTa                                                                                                                                                      >  1:218578/1‑250 (MQ=255)                                                                                                    gtgAAAATTTTAACATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGAAGTTCTTCGAAACTTCCCGCAGCAGGTCATCAGGGGTGATTTCAGGATCGGAGATAAAGGTCAGGTTCGGGGAAACCTGCTTAAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGACGATATTGGCGTGCATAAAGGCGTCTGGTAGGGGTTTGTCGCGGTAACGcccgacacg                                                                                                                                          >  2:228290/1‑244 (MQ=255)                                                                                                            tttATCATGCACCGAACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAACGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGAAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGATTTATTc                                                                                                                                  >  2:104930/1‑251 (MQ=255)                                                                                                             ttATCATGCACCGGACCGTTGGCGTTGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCg                                                                                                                                 >  2:49641/1‑251 (MQ=255)                                                                                                                     cACCGAACAGTTGGCGATGAGAGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCCTCAGGGGTGATTTCAGGGTCGTAGATCAAGGTCAGGTTCGGGGAAACCTGATTAAACTCTGAGTCTGACCGTAAGATCGCGCGGGGAATGGGCGAATAAGACGGGCCGATATTGGCGTGCATCAAGGCGTCTGGCAGGGTGCTGTCGAGGTAACGCCAAAAACGTTTTATTCGAATATCg                                                                                                                          >  2:185897/1‑250 (MQ=255)                                                                                                                           aCCGTTGGCCATGAGCGGTATGCTACATTTACCGATGTTCTTCGACACTTCCGGCAGCAGGCCATAAGGGGTGCTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTACCGCCAGAAACGTTTTATTCGAATATCGATCTCGt                                                                                                                   >  2:155431/1‑251 (MQ=255)                                                                                                                             cGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCt                                                                                                                 >  1:18141/1‑251 (MQ=255)                                                                                                                                     aTGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGtt                                                                                                           >  1:85858/1‑249 (MQ=255)                                                                                                                                         gCGGTTTGCTACATGCACGGAGGTTCTTCGCCACTTCCAGCAGCAGTGCATCAGGGTTGATTTCAGGATCGTGGATAAAGGTCAGGTTCGGGGAAACCTGCGTCATCTCTGCGTCTGCACGTAAGATCTCGCGGGTAAGGGGCGAATCAGGCGTTCCGGTAGTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAAt                                                                                                      <  1:139449/250‑1 (MQ=255)                                                                                                                                             tttGCTACATTCACAGATGTTATTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTAGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATAAGACGGGCCGATATTGGAGAGCATAAAGGCGTCTGGCAGGGTTCTGTTGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTa                                                                                                  >  2:91261/1‑250 (MQ=255)                                                                                                                                                  tccattcccggATGTTCTTAGCCACTTCCAGCAGCCGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGCAACATGCTTCAACTCTGCGTCTGACCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGCGGTAACGCCCGAAACGTTTTATTCGAATATGGATCTCGTATTGTGTTAGAATTCTAACATat                                                                                            >  2:20057/11‑250 (MQ=255)                                                                                                                                                        cacaGATGTTCTTCGCCACTTCCCGCAGAAGGTCATCAGGGGTGATTTCAGGCTCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGc                                                                                      >  2:233310/1‑251 (MQ=255)                                                                                                                                                                         aCTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTg                                                                     <  2:260820/251‑1 (MQ=255)                                                                                                                                                                                      ggTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                                                        >  1:85516/1‑251 (MQ=255)                                                                                                                                                                                                ggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg                                              <  2:85858/251‑1 (MQ=255)                                                                                                                                                                                                 gTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGAGCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                             <  1:253880/251‑1 (MQ=255)                                                                                                                                                                                                 gTGATTTCAGGATCGTAGATAAAGGGCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCGGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGGCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTGGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                                             <  1:242243/251‑1 (MQ=255)                                                                                                                                                                                                    aTTTCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                                          >  1:156406/1‑251 (MQ=255)                                                                                                                                                                                                       tCAGGATCGTAGATAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGATCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGTTCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTCGACCGGCATCGATGTGACCg                                        >  1:247264/1‑250 (MQ=255)                                                                                                                                                                                                                    tAAAGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                          >  1:95073/1‑251 (MQ=255)                                                                                                                                                                                                                     aaaGGTCAGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg                          <  2:85516/250‑1 (MQ=255)                                                                                                                                                                                                                            aGGTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt                  <  1:104930/251‑1 (MQ=255)                                                                                                                                                                                                                              gTTCGGGGAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt                  >  1:16105/1‑249 (MQ=255)                                                                                                                                                                                                                                   gggAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                             <  1:82273/199‑1 (MQ=255)                                                                                                                                                                                                                                   gggAAACCTGCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAATATCGATCTCGTCTTGTGTTAGAATTCTAACATATGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                                             >  2:82273/3‑201 (MQ=255)                                                                                                                                                                                                                                            gCTTCAACTCTGCGTCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATAGTGGCGTGCATAAAGGCGTCTGGCATGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAAGATCGATCGCGTCTTGTGTTAGAATTCTAACATATGGTGGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGt  <  2:105625/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              GCTCGGCAATGGCTTTCAGGTTAAGGCGTACCAGCGTG‑CTGCCACCACCCGCCAGCGGCAGTGAGTTGTAACAGCTCACAATCCCGTAGCCCCCTTTTGTGAAAATTTTATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGT  >  NC_000913/4616288‑4616772

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