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

GGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCA  >  NC_000913/567898‑568131                                                                                                           |                                                                                                                                gggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAAc                                                                                                 <  1:780023/139‑1 (MQ=11) gggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAAc                                                                                                 <  1:279054/139‑1 (MQ=11)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                 >  1:638383/1‑136 (MQ=255)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                 <  2:638383/136‑1 (MQ=255)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                              >  2:124659/1‑139 (MQ=255)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCg                                                                              >  2:809368/1‑139 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGc                                                                            >  1:456301/1‑139 (MQ=255)                         aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCAt                                                                                           >  1:947253/1‑121 (MQ=21)                         aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCAt                                                                                           <  2:947253/121‑1 (MQ=21)                         aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTaa                                                                         >  2:201209/1‑139 (MQ=255)                          aTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTaaa                                                                        >  2:488254/1‑139 (MQ=255)                                ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGgcg                                                                  <  1:488254/139‑1 (MQ=255)                                       cGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAg                                                           <  1:743240/139‑1 (MQ=255)                                                cggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTgg                                                  >  1:828239/1‑139 (MQ=255)                                                    acaGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGt                                              >  2:132593/1‑139 (MQ=255)                                                            tgtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                 <  2:39407/96‑1 (MQ=255)                                                            tgtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCt                                                                                 >  1:39407/1‑96 (MQ=255)                                                            tgtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTa                                      <  1:201209/139‑1 (MQ=255)                                                             gtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTaa                                     >  1:198230/1‑139 (MQ=255)                                                                  ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCt                                >  1:734627/1‑139 (MQ=255)                                                                     cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGc                                                                                  <  2:1061473/86‑1 (MQ=255)                                                                     cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGc                                                                                  >  1:1061473/1‑86 (MQ=255)                                                                      tCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTa                                      <  1:931612/129‑1 (MQ=255)                                                                      tCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTa                                      >  2:931612/1‑129 (MQ=255)                                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          <  1:732913/139‑1 (MQ=255)                                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          <  1:124659/139‑1 (MQ=255)                                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          <  2:456301/139‑1 (MQ=255)                                                                        aGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          <  2:734627/139‑1 (MQ=255)                                                                          tCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACg                        >  2:351755/1‑139 (MQ=255)                                                                              gAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGc                                                                                             <  2:152446/66‑1 (MQ=255)                                                                              gAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGc                                                                                             >  1:152446/1‑66 (MQ=255)                                                                              gAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGgtgt                    >  1:10101/1‑139 (MQ=255)                                                                              gAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGgtgt                    >  1:987185/1‑139 (MQ=255)                                                                              gAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGgtgt                    >  1:547728/1‑139 (MQ=255)                                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc              <  2:987185/139‑1 (MQ=255)                                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc              <  1:132593/139‑1 (MQ=255)                                                                                    tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTAcc              <  2:547728/139‑1 (MQ=255)                                                                                                aaCCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCa  <  1:701400/139‑1 (MQ=255)                                                                                                     cGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGt                                  <  2:48890/102‑1 (MQ=255)                                                                                                     cGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGt                                  >  1:48890/1‑102 (MQ=255)                                                                                                      gCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          >  2:810779/1‑109 (MQ=255)                                                                                                      gCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGaa                          <  1:810779/109‑1 (MQ=255)                                                                                                       cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGtt                                             >  1:475606/1‑89 (MQ=255)                                                                                                       cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGtt                                             <  2:475606/89‑1 (MQ=255)                                                                                                           |                                                                                                                                GGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCA  >  NC_000913/567898‑568131

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: ‑