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%	108.7 / NA	30	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 (15/15);  total (15/15)

ATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGA  >  NC_000913/4499427‑4499907                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       atccctcgcCACGGGAATAACTGCGGCAGCGGGTTTGGAATGGGTTAAGTGATAATAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGgt                                                                                                                                                                                                                                        <  2:189476/247‑1 (MQ=255)       cgcCACGGGAATATCTGCGGCAGCGGGATTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGgttagt                                                                                                                                                                                                                                    <  1:62738/249‑1 (MQ=255)           aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCtt                                                                                                                                                                                                                              >  1:63194/1‑251 (MQ=255)                 atatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTACACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCGCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTg                                                                                                                                                                                                                        <  1:93166/251‑1 (MQ=255)                       gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGc                                                                                                                                                                                                                   <  2:63194/250‑1 (MQ=255)                             gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCTATATATGATTACCTGTTATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCAGAACAGGTTAGTTATCTTCGTCTGTCAGCGTATTGa                                                                                                                                                                                                            <  2:230608/251‑1 (MQ=255)                              cGGGCTTGTCATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGa                                                                                                                                                                                                            >  2:152121/1‑250 (MQ=255)                              cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTTGAAAGATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATTATTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATGAGGAAGAATAATTGGGTGAGGTTTCAAACTGTCTGCCGTTATATTCCTCTCGCCTTCAAAGACCATTTCATGTTAGTTTTCTGCATCGGTGAGCGGATTGaa                                                                                                                                                                                                           <  2:251017/251‑1 (MQ=255)                                                         aGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATAcc                                                                                                                                                                                <  2:49477/251‑1 (MQ=255)                                                           aTGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTGCGTCGGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCgg                                                                                                                                                                              <  2:156118/251‑1 (MQ=255)                                                                      tataGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAACTACCGGTCATCCACTg                                                                                                                                                                    >  2:68717/1‑250 (MQ=255)                                                                           gggCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGaaaaa                                                                                                                                                               >  2:228220/1‑250 (MQ=255)                                                                                         tgtgCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGAc                                                                                                                                                >  2:211407/1‑251 (MQ=255)                                                                                                  atGAACACGGATATTTTATACCGAATACTTTTCAGTACAAGCGAATGTATTAAAGTCAGCCATATATGATAACATGTGATCAGGAAGAGTAAGTGGCTGAGGTTTCAAACTGTCAGCCGGTATATTACTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACACATACCGGTCATCCCCTGAAAAATACATGCGAAATGACGGACACtt                                                                                                                                        >  2:220041/1‑249 (MQ=255)                                                                                                                 ttATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACg                                                                                                                         <  1:121613/250‑1 (MQ=255)                                                                                                                  tataCAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATCCCTGCGAAATGACGGGCACTGATTTTTCTACTGACgg                                                                                                                        <  1:124897/250‑1 (MQ=255)                                                                                                                         aaTACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACTTGTTATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGCATGGCGGTATATTTTCAGTACAAATACCGGTCATCCCCTGAAAAATACCTTCGAAATAAATGGAACTGATTTTTCTACTAACGGACtgatg                                                                                                                 >  2:86674/1‑250 (MQ=255)                                                                                                                            aCTTTTCAGTCCAAGCGACTGTCTTACAGTATGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCGGCGGATTGAATGGCGGTATATTTTCAGTACAACTACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGAc                                                                                                             >  2:161615/1‑251 (MQ=255)                                                                                                                              ttttCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGt                                                                                                           >  2:40314/1‑251 (MQ=255)                                                                                                                                                   ttACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGCTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGCAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGtcaaaat                                                                                      >  1:133760/1‑245 (MQ=255)                                                                                                                                                           tGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGgaattattagggggc                                                                              >  1:54063/1‑240 (MQ=255)                                                                                                                                                            gCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGc                                                                               <  1:49535/249‑1 (MQ=255)                                                                                                                                                                        aaCCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAACACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGGAATAATCACTAGCacaccgtgtcccac                                                                 >  2:151233/2‑241 (MQ=255)                                                                                                                                                                          ccTGTGCTCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAAAAATCAGTGGCACCATGTTGCCTcccc                                                               >  1:225317/1‑250 (MQ=255)                                                                                                                                                                                    ggAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCCGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTGATAAACAATGGCCCCCTGTAGTCTCCCTGcaccatacc                                                     >  1:96248/1‑245 (MQ=255)                                                                                                                                                                                                                 gTCTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTc                          <  1:202008/249‑1 (MQ=255)                                                                                                                                                                                                                   cTGCCGGTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTTTAGCCTCCCCGCACCTGACCCTAATGATGAACAGGATTGGCACTGTCGCTg                      >  1:162803/1‑251 (MQ=255)                                                                                                                                                                                                                        ggTATATTCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGa                 <  1:152121/251‑1 (MQ=255)                                                                                                                                                                                                                               tCCTCTCTCCCTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGccc           <  2:162803/250‑1 (MQ=255)                                                                                                                                                                                                                                       cccTCAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGa  <  1:228220/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       ATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGA  >  NC_000913/4499427‑4499907

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