Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,499,661	T→C	intergenic (+161/‑1599)	insD6 → / → yjgZ	IS2 insertion element protein InsB/uncharacterized protein YjgZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,499,661	0	T	C	100.0%	113.4 / NA	33	intergenic (+161/‑1599)	insD6/yjgZ	IS2 insertion element protein InsB/uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/17);  total (16/17)

GCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACA  >  NC_000913/4499417‑4499894                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCGCTCTCCCTCAAAAAc                                                                                                                                                                                                                                     <  1:135473/251‑1 (MQ=255)     tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCAt                                                                                                                                                                                                                                  >  1:316703/1‑250 (MQ=255)                 cgcCACGGGAATATCTGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGgttagt                                                                                                                                                                                                                       <  2:414260/249‑1 (MQ=255)                       gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCg                                                                                                                                                                                                               >  1:147711/1‑251 (MQ=255)                            tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGt                                                                                                                                                                                                          >  2:116966/1‑251 (MQ=255)                                  cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTTAGCgg                                                                                                                                                                                                    <  1:103434/251‑1 (MQ=255)                                              tGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCTCCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCg                                                                                                                                                                                         <  2:191726/250‑1 (MQ=255)                                                             gATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTAc                                                                                                                                                                          <  1:419493/250‑1 (MQ=255)                                                                   aGATGTCTGGAAATATAGGGGCAAATCCCCCTTGTGCTGATATGAATACTGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATATCAGGTTAGTTAACTTCGTCTGTCAGCGGATTGAATGGCGTTATATTTTCAGTACAAATAcc                                                                                                                                                                   <  2:255930/251‑1 (MQ=255)                                                                    gATGTCTGGAAATATAGGGGCATATCCACCTGGTGCTGATAGGAATTCGGATATTTTATACAGAATACTATTCAGTTCAAGCGACTGTATTACAGTCTGTCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGGCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCg                                                                                                                                                                  <  2:369330/251‑1 (MQ=255)                                                                    gATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCg                                                                                                                                                                  <  1:295842/251‑1 (MQ=255)                                                                      tGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGt                                                                                                                                                                >  1:415702/1‑251 (MQ=255)                                                                                    ggggAAAATCCACCTTGTGCTGATATGAATACGGAGATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGaaaaa                                                                                                                                                  <  2:291792/251‑1 (MQ=255)                                                                                     gggCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAATACAAATACCCGTCATCCACTTaaaaa                                                                                                                                                  >  2:9828/1‑250 (MQ=255)                                                                                                ccTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAg                                                                                                                                                                                                       >  1:319995/1‑186 (MQ=255)                                                                                                ccTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTGCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAg                                                                                                                                                                                                       <  2:319995/186‑1 (MQ=255)                                                                                                  ttGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCCTCCCATGAAAAATACCTACGAAAtt                                                                                                                                     >  2:212239/1‑249 (MQ=255)                                                                                                   tgtgCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATg                                                                                                                                     <  2:147711/249‑1 (MQ=255)                                                                                                                      atatTTTATACAGAATACTTTTCAGGACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTAAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAACAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGACTGGCGGTATATTTTCAGAACAAATACCGGTCATCCACTCAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACt                                                                                                                >  1:424577/1‑251 (MQ=255)                                                                                                                         ttttATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGAc                                                                                                             >  2:437805/1‑251 (MQ=255)                                                                                                                                             aGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGt                                                                                         <  1:323182/251‑1 (MQ=255)                                                                                                                                                                cAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCGGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAg                                                                      <  2:218577/251‑1 (MQ=255)                                                                                                                                                                cAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATAGTCCTCTCTCCCTCAAAAACCATAACATGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCAGCCACAGAAAAATACCTGCGAAATGACGGGCACTGATTTTGCTACTGACGGGCTGATGAGACGTTATGTCACATGCGCTAATAATCAg                                                                      <  2:318872/251‑1 (MQ=255)                                                                                                                                                                    cTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTCTATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCt                                                                                                                                                                                                            >  2:370752/1‑113 (MQ=255)                                                                                                                                                                    cTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCt                                                                                                                                                                                                            <  1:370752/113‑1 (MQ=255)                                                                                                                                                                      gCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGccc                                                                >  1:448095/1‑249 (MQ=255)                                                                                                                                                                              tGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCAACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGAGGGTATGTCAGTGGCGGTAATTATCAGAGGGACACTGTAgg                                                        >  2:210739/1‑250 (MQ=255)                                                                                                                                                                                 ttACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGCTTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGGAATTAACTGTGGGACACTGtgtgccc                                                     >  2:75305/1‑245 (MQ=255)                                                                                                                                                                                                   aaTAAGTGGATGAGGTTGCAAACTGGCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCTGTATATTTTCAGTACAAATACCGGTCATCCACTGGAAAATGCCTGCGAAATGACTGGCTCTGATTTTTCTACTGACTGGCTGATGAGACGGTGGGTCGCTGGCGGTAATGGTTTTTGGCACACTGTAGCCTCCCTGCACATGACCATaa                                    <  1:55703/250‑1 (MQ=255)                                                                                                                                                                                                       aGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAAATATCAGTGGCAAACTGGAGGCCCCCCGCACATGACCATAaagatg                               >  1:428927/1‑251 (MQ=255)                                                                                                                                                                                                          tgCTGAGTTTTCTAACTGTCTGTCGGTGTATTCCTCTTTCCCTCAAAAACCATAACTGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGTTATATTTTCAGTACAAATACCGGTCATCCTCTGAAAAATACCTGCGAAATGACGGGCACTGATTTCTCTGCTGACGGGCTTATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGGACACGGTAGCCTCCCTGCACATGATCATAATGATGaa                             <  1:69578/249‑1 (MQ=255)                                                                                                                                                                                                                  tttCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCTTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCAACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGCCGGGCTGATGAGACGTTATGTCACTGGCGGTTATAATCAGAGGGACACATTAGCCTCCCTGCACATGACCATAatgaatgacaagatatg                    >  1:540501/1‑238 (MQ=255)                                                                                                                                                                                                                                     atatTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGGAAATATCAGTGGCGCACTGTAGCCTCCCCGCACCAGACCATAATGATGCAAAGGATTTGCACCCTCGCTGACCcaca  >  2:377224/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          GCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACA  >  NC_000913/4499417‑4499894

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