Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,996,372	C→T	*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	2,996,372	0	C	T	100.0%	172.8 / NA	55	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/27);  total (28/27)

CTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTG  >  NC_000913/2996237‑2996496                                                                                                                                        |                                                                                                                             cTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTtat                                                                                                                           <  1:630984/139‑1 (MQ=255)    ataaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAtt                                                                                                                        >  1:751077/1‑139 (MQ=255)      aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                      >  2:692696/1‑139 (MQ=255)      aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                      >  2:343695/1‑139 (MQ=255)        ataatCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                    <  2:517706/139‑1 (MQ=255)           atCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGAc                                                                                                                 >  2:735436/1‑139 (MQ=255)              ccTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATc                                                                                                              >  1:737022/1‑139 (MQ=255)                  aGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                          >  1:652694/1‑139 (MQ=255)                        tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACAt                                                                                                               <  1:261662/128‑1 (MQ=255)                        tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACAt                                                                                                               >  2:261662/1‑128 (MQ=255)                        tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  2:737022/139‑1 (MQ=255)                        tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  2:652694/139‑1 (MQ=255)                         cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTta                                                                                                                            >  2:19750/1‑114 (MQ=255)                         cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                   >  2:606504/1‑139 (MQ=255)                         cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                   >  1:503917/1‑139 (MQ=255)                           ggAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTaa                                                                                                 <  2:706945/139‑1 (MQ=255)                            gAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAGCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                                >  2:431973/1‑139 (MQ=255)                             aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                               <  1:62415/139‑1 (MQ=255)                              aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGTTTTGCCCTTATATTTCCAGa                                                                                                                  >  1:506355/1‑119 (MQ=255)                              aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGa                                                                                                                  <  2:506355/119‑1 (MQ=255)                              aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    >  1:445842/1‑133 (MQ=255)                              aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  2:445842/133‑1 (MQ=255)                              aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCC‑‑ACATCTGTTATCACTTAACCCAt                                                                                            >  2:480802/1‑139 (MQ=255)                                 ttATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                           >  2:254075/1‑139 (MQ=255)                                          aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                  >  1:451065/1‑139 (MQ=255)                                              cgcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    <  2:91718/117‑1 (MQ=255)                                              cgcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                    >  1:91718/1‑117 (MQ=255)                                                               tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         >  1:70986/1‑95 (MQ=255)                                                               tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                         <  2:70986/95‑1 (MQ=255)                                                                 tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  1:735436/139‑1 (MQ=255)                                                                 tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  1:581254/139‑1 (MQ=255)                                                                 tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  1:343695/139‑1 (MQ=255)                                                                      aaaaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  2:522894/134‑1 (MQ=255)                                                                      aaaaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           >  1:522894/1‑134 (MQ=255)                                                                        aaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                    <  2:293340/139‑1 (MQ=255)                                                                        aaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                    <  1:431973/139‑1 (MQ=255)                                                                                aCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                            <  2:503917/139‑1 (MQ=255)                                                                                  cACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                          <  2:451065/139‑1 (MQ=255)                                                                                        gTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGc                                    <  2:137588/139‑1 (MQ=255)                                                                                           aaCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCAAGCGATAACCCAGCGCACTATGCGGa                                 >  1:162343/1‑139 (MQ=255)                                                                                            aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgag                                                         <  2:380578/114‑1 (MQ=255)                                                                                            aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgag                                                         >  1:380578/1‑114 (MQ=255)                                                                                               gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  2:382385/109‑1 (MQ=255)                                                                                               gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           >  1:382385/1‑109 (MQ=255)                                                                                               gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                             <  2:751077/139‑1 (MQ=255)                                                                                                aTCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCC‑TTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                           >  1:20969/1‑139 (MQ=255)                                                                                                 taaaTGGAAATAAATTCCAACTGTTTATCGGATTTGCCCGTATATTTCCAGACATCTGTTATCACTTAACCCAATACAAGCCCACTGCCGCAGATCTGCCCGTGGCCAGCGATAACCCATCGCACTATGCCGATGCCAt                           <  1:650973/137‑1 (MQ=14)                                                                                                 tCAATGGAAGTAAATTCCCATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           >  1:446804/1‑107 (MQ=255)                                                                                                 tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                           <  2:446804/107‑1 (MQ=255)                                                                                                   aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                    >  1:81639/1‑48 (MQ=255)                                                                                                   aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                    <  2:81639/48‑1 (MQ=255)                                                                                                     tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcga                                                      <  1:782958/108‑1 (MQ=255)                                                                                                     tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcga                                                      >  2:782958/1‑108 (MQ=255)                                                                                                     tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt                       >  2:132325/1‑139 (MQ=255)                                                                                                                       gTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct     <  1:254075/139‑1 (MQ=18)                                                                                                                          tATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTg  >  1:232000/1‑139 (MQ=11)                                                                                                                                        |                                                                                                                             CTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTG  >  NC_000913/2996237‑2996496

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