Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,187,415	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	3,187,415	0	G	A	100.0%	217.8 / NA	67	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/29);  total (38/29)

AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCA  >  NC_000913/3187292‑3187550                                                                                                                            |                                                                                                                                        agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                          <  1:1124985/139‑1 (MQ=17) agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                          <  1:1407861/139‑1 (MQ=17) agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                          <  2:12934/139‑1 (MQ=17)        ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATcatcat                                                                                                                   >  2:1615569/1‑139 (MQ=38)           gAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTgg                                                                                                                <  2:884650/139‑1 (MQ=255)                  ataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                         <  1:1346225/139‑1 (MQ=255)                   taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                        <  1:581178/139‑1 (MQ=255)                   taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTAGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                        >  1:1110490/1‑139 (MQ=255)                        aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTtg                                                                                                             >  2:1231272/1‑129 (MQ=255)                        aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTtg                                                                                                             <  1:1231272/129‑1 (MQ=255)                           ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                                >  2:1602595/1‑139 (MQ=255)                           ggCATCCGCATAGTGCGCTGGGTTATCGATCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                                >  1:1041848/1‑139 (MQ=255)                                   cATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGACAGCGGGCTGGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGa                                                                                        >  1:993676/1‑139 (MQ=255)                                        tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATg                                                                                               >  2:53513/1‑127 (MQ=255)                                        tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATg                                                                                               <  1:53513/127‑1 (MQ=255)                                          cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                 >  2:1398072/1‑139 (MQ=255)                                          cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                 >  2:1615705/1‑139 (MQ=255)                                             tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTtg                                                                                                             <  2:326996/108‑1 (MQ=255)                                             tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTtg                                                                                                             >  1:326996/1‑108 (MQ=255)                                                     cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                      <  1:960740/139‑1 (MQ=255)                                                     cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                      <  1:1602595/139‑1 (MQ=255)                                                     cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                      <  1:1398072/139‑1 (MQ=255)                                                       ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    <  2:993676/139‑1 (MQ=255)                                                       ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    >  2:1361777/1‑139 (MQ=255)                                                         cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGaa                                                                  >  2:1240358/1‑139 (MQ=255)                                                         cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGC‑AATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAg                                                                 >  1:1138078/1‑139 (MQ=255)                                                           ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                        >  2:1303067/1‑99 (MQ=255)                                                           ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                        <  1:1303067/99‑1 (MQ=255)                                                                    tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGCTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGAt                                                                           <  1:1020200/119‑1 (MQ=255)                                                                    tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGAt                                                                           >  2:1020200/1‑119 (MQ=255)                                                                           ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAAc                                                >  2:1302567/1‑139 (MQ=255)                                                                            gCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    <  1:1142091/118‑1 (MQ=255)                                                                            gCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    >  2:1142091/1‑118 (MQ=255)                                                                                       gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                    <  2:589034/139‑1 (MQ=255)                                                                                       gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                    <  1:1615569/139‑1 (MQ=255)                                                                                        tAATGGGTTAAGTGATAATAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTccc                                   >  1:895712/1‑139 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTg                                                     >  2:461256/1‑116 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTg                                                     <  1:461256/116‑1 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                              >  2:482600/1‑139 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                              >  2:1028102/1‑139 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                              >  2:145167/1‑139 (MQ=255)                                                                                               ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAAt                                                                                      <  1:1539317/81‑1 (MQ=255)                                                                                               ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAAt                                                                                      >  2:1539317/1‑81 (MQ=255)                                                                                               ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGa                            >  2:678621/1‑139 (MQ=255)                                                                                                 aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                 <  1:258401/84‑1 (MQ=255)                                                                                                 aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                 >  2:258401/1‑84 (MQ=255)                                                                                                 aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGaaa                          >  2:841036/1‑139 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                         >  1:198010/1‑139 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                         >  1:495422/1‑139 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                         >  1:1374042/1‑139 (MQ=255)                                                                                                     gATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAAt                                                                                      <  2:872283/75‑1 (MQ=255)                                                                                                     gATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAAt                                                                                      >  1:872283/1‑75 (MQ=255)                                                                                                     gATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGa                      <  1:1302567/139‑1 (MQ=255)                                                                                                       tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATcatcat                                                                                                                   <  1:260892/44‑1 (MQ=255)                                                                                                       tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATcatcat                                                                                                                   >  2:260892/1‑44 (MQ=255)                                                                                                       tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTa                                            >  2:1600517/1‑115 (MQ=255)                                                                                                       tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTa                                            <  1:1600517/115‑1 (MQ=255)                                                                                                        aaCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATgg                   >  1:1088335/1‑139 (MQ=255)                                                                                                                tCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    >  1:596374/1‑82 (MQ=255)                                                                                                                tCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                    <  2:596374/82‑1 (MQ=255)                                                                                                                tCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTAc           >  2:585132/1‑139 (MQ=255)                                                                                                                 cTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTAcc          <  1:145167/139‑1 (MQ=255)                                                                                                                    gAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGt       >  1:1099443/1‑139 (MQ=255)                                                                                                                     aaaTATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTc      >  2:413943/1‑139 (MQ=255)                                                                                                                     aaaTATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTc      <  2:1110490/139‑1 (MQ=255)                                                                                                                     aaaTATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTc      <  1:413943/139‑1 (MQ=255)                                                                                                                         ataAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCa  >  1:1516963/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                        AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCA  >  NC_000913/3187292‑3187550

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