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%	86.3 / NA	28	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/10);  total (18/10)

GGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTA  >  NC_000913/567899‑568128                                                                                                          |                                                                                                                             ggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAAcc                                                                                             >  1:71991/1‑139 (MQ=14)                ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                              <  1:219232/139‑1 (MQ=255)                   tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                           >  2:591145/1‑139 (MQ=255)                   tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                           >  2:530877/1‑139 (MQ=255)                   tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                           >  2:4928/1‑139 (MQ=255)                     tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGc                                                                         >  1:657748/1‑139 (MQ=255)                          tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTaaaa                                                                    >  1:546949/1‑139 (MQ=255)                                      cGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAg                                                        >  2:70193/1‑139 (MQ=255)                                        gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGa                                                      >  2:154944/1‑139 (MQ=255)                                         tGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGaa                                                     >  1:742371/1‑139 (MQ=255)                                           gcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGACCATTTGGAAAAAAAAAACCTCGCTTAGGCTTGGTTCCTAATTCCGGGGGAAGGTCCACCCAGCTTAAACAAGCCTGCGGCAAAAATGCCGCGGGAAGGAATc                                                   >  1:307765/1‑139 (MQ=9)                                           gcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATc                                                   >  1:18243/1‑139 (MQ=255)                                                   acaGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGt                                           >  2:332458/1‑139 (MQ=255)                                                          gtgtGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGt                                    <  2:71991/139‑1 (MQ=255)                                                                  gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCgg                                                                          <  1:606766/93‑1 (MQ=255)                                                                  gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCgg                                                                          >  2:606766/1‑93 (MQ=255)                                                                       aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                       <  1:4928/139‑1 (MQ=255)                                                                       aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                       <  1:480094/139‑1 (MQ=255)                                                                       aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                       <  2:546949/139‑1 (MQ=255)                                                                                   tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc           >  2:421719/1‑139 (MQ=255)                                                                                   tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc           >  1:205376/1‑139 (MQ=255)                                                                                   tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc           >  1:100718/1‑139 (MQ=255)                                                                                        aaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTc      >  1:701154/1‑139 (MQ=255)                                                                                           cAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGt                               >  2:115376/1‑111 (MQ=255)                                                                                           cAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGt                               <  1:115376/111‑1 (MQ=255)                                                                                            aaGAACCTCGCTTAAGCCTGTGTCCATATTTCGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTa  <  2:18243/139‑1 (MQ=255)                                                                                            aaGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTa  <  2:307765/139‑1 (MQ=255)                                                                                             aGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGc                                                                            >  1:519610/1‑64 (MQ=255)                                                                                             aGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGc                                                                            <  2:519610/64‑1 (MQ=255)                                                                                                          |                                                                                                                             GGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTA  >  NC_000913/567899‑568128

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