Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,519,014	G→A	pseudogene (294/504 nt)	insB1 →	IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,519,014	0	G	A	100.0%	165.9 / NA	62	pseudogene (294/504 nt)	insB1	IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/24);  total (38/24)

GTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518884‑4519144                                                                                                                                   |                                                                                                                                   gTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTACAAATGACGGa                                                                                                                            <  2:707309/139‑1 (MQ=255)      cgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       <  1:586090/139‑1 (MQ=255)           gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:274011/139‑1 (MQ=255)           gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  1:287394/1‑139 (MQ=255)           gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:846835/139‑1 (MQ=255)           gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:875128/139‑1 (MQ=255)                tttACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTa                                                                                                             <  1:454230/139‑1 (MQ=255)                         aTGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                                   >  1:218951/1‑139 (MQ=11)                         aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                                    >  2:597301/1‑139 (MQ=255)                           gACAGGCTCCGGAAGACGGTCGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCt                                                                                                 >  1:181639/1‑139 (MQ=11)                                gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                                            >  1:341556/5‑139 (MQ=11)                                 aTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                                     <  2:594956/114‑1 (MQ=255)                                 aTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                                     >  1:594956/1‑114 (MQ=255)                                   ccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaa                                                                                         >  2:997019/2‑139 (MQ=11)                                   aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                                          >  2:460156/1‑139 (MQ=255)                                     ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                                       >  1:465987/5‑139 (MQ=14)                                     ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                                       >  1:888494/1‑139 (MQ=14)                                       aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGc                                                                                      >  1:965920/1‑139 (MQ=255)                                          aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                                   >  1:117343/1‑139 (MQ=255)                                          aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                                   >  1:583896/1‑139 (MQ=255)                                            ggTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                                >  1:445884/1‑139 (MQ=14)                                            ggTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                                >  2:987271/1‑139 (MQ=14)                                                gtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTa                                                                             >  2:369210/1‑139 (MQ=255)                                                    cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  1:973450/139‑1 (MQ=14)                                                    cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        >  2:973450/1‑139 (MQ=14)                                                    cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  2:181639/139‑1 (MQ=14)                                                     gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                        <  2:965920/139‑1 (MQ=255)                                                      cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                                       >  2:594948/1‑139 (MQ=255)                                                                gTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTacac                                                             <  1:460156/139‑1 (MQ=255)                                                                   aaCGCACTATGGCGACGCTGGAGCGTCTTCTGCGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    >  2:628590/1‑80 (MQ=18)                                                                   aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  1:628590/80‑1 (MQ=18)                                                                      gCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc                                                                    >  1:133321/1‑125 (MQ=17)                                                                      gCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc                                                                    <  2:133321/125‑1 (MQ=14)                                                                           ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTg                                                  >  2:992542/1‑139 (MQ=255)                                                                           atGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                 <  1:891650/138‑1 (MQ=11)                                                                            tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  1:235115/139‑1 (MQ=11)                                                                            tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  2:465987/139‑1 (MQ=11)                                                                            tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  1:658400/139‑1 (MQ=11)                                                                                  cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  <  2:32/97‑1 (MQ=21)                                                                                  cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  >  1:32/3‑99 (MQ=21)                                                                                  cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                                          >  1:811915/3‑135 (MQ=11)                                                                                    cTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATa                                         >  2:880329/1‑139 (MQ=255)                                                                                    cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        >  1:669375/1‑107 (MQ=18)                                                                                    cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                                        >  2:928528/1‑133 (MQ=14)                                                                                    cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCATTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  2:669375/107‑1 (MQ=17)                                                                                        aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                     >  1:802655/1‑139 (MQ=255)                                                                                            tcttctGAGCCTGCTGTCGACCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                                           <  2:526159/81‑1 (MQ=255)                                                                                            tcttctGAGCCTGCTGTCGACCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                                           >  1:526159/1‑81 (MQ=255)                                                                                             cttctGAGCCTGCTGTCGGCCTTTGAGGACGTGGTATAAATGACTGATGGCTGGCCGCTGTATGAATCACGCCTGACGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatc                                >  1:564506/1‑139 (MQ=255)                                                                                               tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                              >  2:192289/1‑139 (MQ=255)                                                                                               tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                              >  2:723288/1‑139 (MQ=255)                                                                                                 tGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTgag                            >  1:204826/1‑139 (MQ=255)                                                                                                      ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCa                                                                          >  1:921977/1‑88 (MQ=255)                                                                                                      ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCa                                                                          <  2:921977/88‑1 (MQ=255)                                                                                                         gctgtggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATatacgc                                                           >  1:117148/6‑95 (MQ=14)                                                                                                             cGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGc                >  2:960770/1‑139 (MQ=255)                                                                                                               gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                                      <  2:800632/83‑1 (MQ=255)                                                                                                               gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                                      >  1:800632/1‑83 (MQ=255)                                                                                                                         tCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg    <  1:597301/139‑1 (MQ=255)                                                                                                                           gTGGTATAAATGACGGATCGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGa  >  2:587609/1‑139 (MQ=255)                                                                                                                            tGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                                <  2:834194/60‑1 (MQ=255)                                                                                                                            tGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                                >  1:834194/1‑60 (MQ=255)                                                                                                                                   |                                                                                                                                   GTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518884‑4519144

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