Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,499,500	G→A	*302* (TAG→TAA)	insD1 →	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,499,500	0	G	A	100.0%	200.4 / NA	64	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/34);  total (30/34)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAA  >  NC_000913/4499375‑4499632                                                                                                                              |                                                                                                                                     gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                         <  2:331886/139‑1 (MQ=11) gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                         <  1:1245712/139‑1 (MQ=11)     agGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgt                                                                                                                     >  1:1374928/1‑139 (MQ=18)      ggCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                    <  2:635735/139‑1 (MQ=21)      ggCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                    <  1:1508872/139‑1 (MQ=21)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGa                                                                                                                >  2:1509660/1‑139 (MQ=38)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     >  1:1114528/1‑139 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  2:638814/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  1:1509660/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  1:1599079/139‑1 (MQ=255)                            tGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATAtttt                                                                                              >  1:1837982/1‑139 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  1:817191/131‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     >  2:817191/1‑131 (MQ=255)                               cATCCGCACAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATATCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                           <  2:1114528/139‑1 (MQ=255)                                           gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                               <  2:1374928/139‑1 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  1:1099694/116‑1 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                              >  2:1471093/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGAAGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     >  2:1099694/1‑116 (MQ=255)                                             gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                    >  2:1074900/1‑132 (MQ=255)                                             gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                    <  1:1074900/132‑1 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     >  1:1863254/1‑113 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                     <  2:1863254/113‑1 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                               >  2:382986/1‑98 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                               <  1:382986/98‑1 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                      >  1:1427062/1‑139 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                      >  2:681950/1‑139 (MQ=255)                                                      tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCg                                                                    >  1:653453/1‑139 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                   <  1:1788786/139‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                                       <  1:744579/83‑1 (MQ=18)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCtt                                                                                                                       >  2:744579/1‑83 (MQ=18)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                     >  2:152920/1‑133 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                     <  1:152920/133‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                    <  1:681950/139‑1 (MQ=255)                                                                        tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                  <  2:653453/139‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                       >  1:350732/1‑129 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                       <  2:350732/129‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                    >  1:927461/1‑132 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                    <  2:927461/132‑1 (MQ=255)                                                                               cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                     >  2:527734/1‑113 (MQ=255)                                                                               cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                     <  1:527734/113‑1 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                           >  2:1242345/1‑120 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                           <  1:1242345/120‑1 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                         >  2:751304/1‑138 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                         <  1:751304/138‑1 (MQ=255)                                                                                      cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCGTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCt                                    >  1:1024199/1‑139 (MQ=255)                                                                                             tGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCa                             <  2:1837982/139‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATATGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAgg                           >  2:482349/1‑139 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                    <  2:1372977/82‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                    >  1:1372977/1‑82 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCa                             >  1:1563568/1‑132 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCa                             <  2:1563568/132‑1 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                         >  1:1252619/1‑66 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATCAATACGGATATTTTATACa                                                                                         <  2:1252619/66‑1 (MQ=37)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                               >  1:1456792/1‑71 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                               <  2:1456792/71‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                   <  1:1713012/83‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                   >  2:1713012/1‑83 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg          >  2:399060/1‑139 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATAAAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCa                                                 <  2:264459/98‑1 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATAAAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCa                                                 >  1:264459/1‑98 (MQ=255)                                                                                                                    tGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                           <  2:1043606/86‑1 (MQ=255)                                                                                                                    tGGAAATATAAGGGCAAAGCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                           >  1:1043606/1‑86 (MQ=255)                                                                                                                      gAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCa    <  2:1427062/139‑1 (MQ=255)                                                                                                                        aaTATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCaaa  <  2:1024199/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                     GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAA  >  NC_000913/4499375‑4499632

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