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%	182.5 / NA	60	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/31);  total (29/31)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624                                                                                                                              |                                                                                                                             gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  2:750423/1‑139 (MQ=11) gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:740290/139‑1 (MQ=11)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGa                                                                                                        >  2:1041079/1‑139 (MQ=38)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:127004/139‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  1:1002672/1‑111 (MQ=11)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:1002672/111‑1 (MQ=11)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  1:1062473/116‑1 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  2:1062473/1‑116 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTta                                                                                     >  2:1003049/1‑139 (MQ=255)                                        aGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:167079/1‑120 (MQ=255)                                        aGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:167079/120‑1 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  1:925591/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  1:826697/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  2:344485/1‑139 (MQ=255)                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:416348/1‑114 (MQ=255)                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:416348/114‑1 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:355606/1‑113 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:355606/113‑1 (MQ=255)                                                     aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  1:750423/139‑1 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  2:201850/139‑1 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:424449/139‑1 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  2:870175/139‑1 (MQ=255)                                                         ctcgcCACGGTAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         <  2:978700/139‑1 (MQ=255)                                                             ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTAt                                                     >  2:57850/1‑139 (MQ=255)                                                              cACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTAtt                                                    <  1:1003049/139‑1 (MQ=255)                                                               aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                   <  2:925591/139‑1 (MQ=255)                                                                cGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTAc                                                  <  2:435626/139‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCATATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                            <  1:354846/139‑1 (MQ=255)                                                                           gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCata                                       <  1:57850/139‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:818201/68‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:818201/1‑68 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                    <  1:716796/71‑1 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:234413/78‑1 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:234413/1‑78 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                 <  2:847701/122‑1 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                 >  1:847701/1‑122 (MQ=255)                                                                                  cGGGCTTGTAATGGGGTAAGTGATAACAGATGTCTGGAAATATAAGGTCAAATCCACCTTGTGCTGATATg                                                                                                    >  2:716796/1‑71 (MQ=17)                                                                                      cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCt                            <  2:826697/139‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:68450/55‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       >  2:68450/1‑55 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           >  2:371983/1‑94 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:371983/94‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        >  2:791959/1‑113 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        >  2:509975/1‑113 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        <  1:509975/113‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                        <  1:791959/113‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGgaagaa              >  2:240196/1‑139 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGgaagaa              >  2:240195/1‑139 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               >  2:57192/1‑100 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               <  1:57192/100‑1 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  1:911700/85‑1 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             >  2:911700/1‑85 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg    >  2:413329/1‑139 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:915411/1‑49 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:915411/49‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:884931/1‑48 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:884931/48‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg  >  2:788484/1‑139 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                          >  1:318616/1‑98 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                          <  2:318616/98‑1 (MQ=255)                                                                                                                              |                                                                                                                             GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624

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: ‑