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%	188.2 / NA	61	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/29);  total (32/29)

ATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/2996240‑2996493                                                                                                                                      |                                                                                                                          ataaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATAtt                                                                                                                     >  1:229563/1‑139 (MQ=255)   aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTc                                                                                                                   >  1:1014053/1‑139 (MQ=255)   aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTc                                                                                                                   >  2:1128675/1‑139 (MQ=255)   aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTc                                                                                                                   >  2:167195/1‑139 (MQ=255)   aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTc                                                                                                                   >  1:479215/1‑139 (MQ=255)       aatCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGa                                                                                                               <  2:445226/139‑1 (MQ=255)              tAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                        >  1:594411/1‑139 (MQ=255)               aGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAATCCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                       <  1:767475/139‑1 (MQ=255)                   cATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCa                                                                                                   >  1:132427/1‑139 (MQ=255)                     tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 >  2:127861/1‑139 (MQ=255)                     tCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 <  2:10541/139‑1 (MQ=255)                          aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  1:167195/139‑1 (MQ=255)                              ttATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                        >  2:543220/1‑139 (MQ=255)                                      aaaCCCGCGTTCTGCAAACCTTTCGCGATAAAAGTGACAACCGCTGAGTAACCAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGccc                                                                                <  2:165330/139‑1 (MQ=255)                                       aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               <  2:944384/139‑1 (MQ=255)                                             cgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                         >  2:319532/1‑139 (MQ=255)                                              gTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                        <  2:930627/139‑1 (MQ=255)                                                  tGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTAATATTGCATGAAATTGGTTCTAATTTAACCCTTACCAAGCCGGCTGCCGCAAat                                                                    >  1:755603/1‑139 (MQ=255)                                                            tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                          >  2:1179079/1‑139 (MQ=255)                                                             ctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                         <  1:543220/139‑1 (MQ=255)                                                              tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        <  1:319532/139‑1 (MQ=255)                                                              tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        <  1:283032/139‑1 (MQ=255)                                                              tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        <  2:1014053/139‑1 (MQ=255)                                                              tcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        <  1:1128675/139‑1 (MQ=255)                                                                          aCAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                            <  1:127861/139‑1 (MQ=255)                                                                                    aGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTaa                                                                                              <  1:1026458/79‑1 (MQ=255)                                                                                    aGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTaa                                                                                              >  2:1026458/1‑79 (MQ=255)                                                                                         aCAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 <  2:1091188/71‑1 (MQ=255)                                                                                         aCAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 >  1:1091188/1‑71 (MQ=255)                                                                                         aCAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                             >  1:226601/1‑139 (MQ=255)                                                                                            gATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  2:894026/139‑1 (MQ=255)                                                                                            gATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  2:132427/139‑1 (MQ=255)                                                                                            gATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          >  1:894026/1‑139 (MQ=255)                                                                                            gATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  2:594411/139‑1 (MQ=255)                                                                                            gATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          >  2:893828/1‑139 (MQ=255)                                                                                             aTCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCTAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATACGGATGCCa                         >  2:542059/1‑139 (MQ=255)                                                                                             aTCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                         <  1:1179079/139‑1 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  2:948536/71‑1 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  1:948536/1‑71 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        >  1:688536/1‑107 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        <  2:688536/107‑1 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                        >  1:929105/1‑139 (MQ=255)                                                                                              tCAATGGAAGTAAATTCCAATTGTGTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAAAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAATATGCGGATGCCAt                        >  2:745141/1‑139 (MQ=255)                                                                                                aattggaGTAAATTCCAATTTTTTATTTGGATTT‑CCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGAGAACCCAGCGCACTATGCGGATGCCATTc                      <  2:755603/133‑1 (MQ=14)                                                                                                aaTGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  2:898814/69‑1 (MQ=255)                                                                                                aaTGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  1:898814/1‑69 (MQ=255)                                                                                                aaTGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  2:1135315/135‑1 (MQ=255)                                                                                                aaTGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          >  1:1135315/1‑135 (MQ=255)                                                                                                aaTGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                      >  1:740509/1‑139 (MQ=255)                                                                                                  tGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt                    >  2:122884/1‑139 (MQ=255)                                                                                                    gAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGa                              >  1:839211/1‑127 (MQ=255)                                                                                                    gAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGa                              <  2:839211/127‑1 (MQ=255)                                                                                                       gTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTCTCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa               >  1:9773/1‑139 (MQ=255)                                                                                                       gTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           <  2:269953/79‑1 (MQ=255)                                                                                                       gTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           >  2:1024997/1‑79 (MQ=255)                                                                                                       gTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           >  1:269953/1‑79 (MQ=255)                                                                                                       gTAAATTCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           <  1:1024997/79‑1 (MQ=255)                                                                                                            ttCCAATTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTc          >  2:347912/1‑139 (MQ=255)                                                                                                                 aTTGTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACg     >  2:704370/1‑139 (MQ=25)                                                                                                                    gTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct  <  2:479215/139‑1 (MQ=18)                                                                                                                    gTTTA‑TTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct  <  2:929105/139‑1 (MQ=18)                                                                                                                                      |                                                                                                                          ATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTA‑TTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/2996240‑2996493

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