Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,004	G→A	*289* (TAG→TAA)	insF1 →	IS3 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,004	0	G	A	100.0%	120.1 / NA	39	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/17);  total (22/17)

TTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGA  >  NC_000913/567914‑568138                                                                                           |                                                                                                                                       ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        >  1:1155814/1‑139 (MQ=255) ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        <  1:362544/139‑1 (MQ=255)    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAAc                                                                                                        <  2:1414151/120‑1 (MQ=11)    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAAc                                                                                                        >  1:1414151/1‑120 (MQ=11)    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                                     >  1:696986/1‑139 (MQ=255)      tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGc                                                                                   >  1:258824/1‑139 (MQ=255)      tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGc                                                                                   >  1:52656/1‑139 (MQ=255)           tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        >  2:529802/1‑129 (MQ=255)           tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        <  1:13928/129‑1 (MQ=255)           tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        >  2:13928/1‑129 (MQ=255)           tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                        <  1:529802/129‑1 (MQ=255)                ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGgcg                                                                         >  1:812732/1‑139 (MQ=255)                      tCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCa                                                                   >  1:1288674/1‑139 (MQ=255)                                cggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTgg                                                         >  2:787632/1‑139 (MQ=255)                                  gcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTc                                                       >  1:247296/1‑139 (MQ=255)                                     caGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGtt                                                    <  2:770072/139‑1 (MQ=255)                                     caGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGtt                                                    <  2:1155814/139‑1 (MQ=255)                                     caGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGtt                                                    <  1:1415491/139‑1 (MQ=255)                                      aGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTg                                                   <  2:247296/139‑1 (MQ=255)                                                gcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGt                                         >  2:27622/1‑139 (MQ=255)                                                  ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATc                                                             >  1:571217/1‑117 (MQ=255)                                                  ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATc                                                             <  2:571217/117‑1 (MQ=255)                                                  ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCt                                       >  1:1274690/1‑139 (MQ=255)                                                   gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTg                                      <  2:258824/139‑1 (MQ=255)                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                                 <  1:101283/139‑1 (MQ=255)                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                                 <  2:1053696/139‑1 (MQ=255)                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                                 <  1:1400310/139‑1 (MQ=255)                                                           ccGGGACAATTTGAAAACAAGATCCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACgg                              <  2:141896/139‑1 (MQ=255)                                                                  aaTTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAAc                                >  2:686603/1‑130 (MQ=255)                                                                  aaTTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAAc                                <  1:686603/130‑1 (MQ=255)                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc                     >  1:498479/1‑139 (MQ=255)                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc                     >  1:858790/1‑139 (MQ=255)                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc                     >  1:773373/1‑139 (MQ=255)                                                                           aCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGc              <  1:27622/139‑1 (MQ=255)                                                                             aaGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTa            >  2:438776/1‑139 (MQ=255)                                                                              aGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTAc           >  2:1437700/1‑139 (MQ=255)                                                                                     cGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGagaa    >  1:412149/1‑139 (MQ=255)                                                                                       cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGagaaga  <  2:696986/139‑1 (MQ=255)                                                                                       cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGagaaga  >  2:700861/1‑139 (MQ=255)                                                                                           |                                                                                                                                       TTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGA  >  NC_000913/567914‑568138

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