Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,496,885	G→A	intergenic (+396/‑1387)	leuX → / → insC6	tRNA‑Leu/IS2 insertion element repressor InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,496,885	0	G	A	100.0%	64.4 / NA	22	intergenic (+396/‑1387)	leuX/insC6	tRNA‑Leu/IS2 insertion element repressor InsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/11);  total (11/11)

GGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTA  >  NC_000913/4496667‑4497094                                                                                                                                                                                                                           |                                                                                                                                                                                                                  ggaATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCg                                                                                                                                                                                   <  1:62054/251‑1 (MQ=255)      aCGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGgag                                                                                                                                                                               <  1:221585/250‑1 (MQ=255)       cGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCAATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGAc                                                                                                                                                                             >  2:201710/1‑251 (MQ=255)                        atCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATaa                                                                                                                                                            <  2:91281/251‑1 (MQ=255)                            aaGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGtt                                                                                                                                                        <  2:60063/251‑1 (MQ=255)                                  gagCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGAACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAAc                                                                                                                                                    >  2:259803/1‑249 (MQ=255)                                     gagccaagccaACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATATGCTTGTCCATCCAAAGGGTTCTAAGTACTGGCGTTTGCAGTACCGTGAGGAGGGAAAGCAAAAAATGCTGGCTCTTGGGGTTTATCCTGACATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGATTGCGAATGGCGGCGATCCGGGAGACAAAAAGAAAAATGATAAGGtttaccttgag                                                                                                                                               <  2:137648/250‑10 (MQ=255)                                              caaCGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACg                                                                                                                                       >  1:221749/1‑250 (MQ=255)                                                      aaGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAAAGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGGCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTaa                                                                                                                              <  2:22953/251‑1 (MQ=255)                                                              ttAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCACATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTCTGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGCGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTcg                                                                                                                      >  2:115800/1‑251 (MQ=255)                                                                          tggtggCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGtggc                                                                                                          >  1:71151/1‑251 (MQ=255)                                                                                              ttGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAACAAATGCTGGCACTTCGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATAGCGTCGATCCGGGAGACAAACAGAAAAATGATAAGGTTGAACAGCGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGCGTGGAATGGCACCAATAAAAAgagg                                                                                      >  1:30120/1‑248 (MQ=255)                                                                                                                  aaGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCAt                                                                    <  1:197406/249‑1 (MQ=255)                                                                                                                                     aCCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGa                                                >  2:24542/1‑250 (MQ=255)                                                                                                                                      ccGTTATGAGGTAAAGCAAAAACTGCTGGCACTTGGGGTTTCTCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAATCTGCTTGCGAATGGCGTCGCTCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCTCGATTGAGTGGCATGGCACCACTAACAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGa                                                >  2:72491/1‑249 (MQ=255)                                                                                                                                       cGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCCGAGCACGTCGTGACCAGACACGTAAGCTGCTTGCGAATCGCGTCGATCCGGGAGACAAACAGAAACATGATAAGGTTGAACAGAGTCAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTCAAAAGTCTTGa                                                >  1:125699/1‑248 (MQ=255)                                                                                                                                        gTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCTTGTGCTAAAAAGTCTTGAAGAt                                            <  1:102004/251‑1 (MQ=255)                                                                                                                                             gAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGAGGATGCCAGAGTACGGCGTGACGAGACGCGTAAGATGCTTGCGAAGGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGGATGGTACCAATAAAAAGTGGTCTGAAGATCTCGCCCATCGTGTGTTAAAAAGTCTGGAATATAATc                                        <  2:182470/250‑1 (MQ=255)                                                                                                                                                         aaaaTGCTGGCACTTGGGGTGTATCCGGAAATCACACTAGCGGATGCAAGAGTACGTCGTGACGAAACGCGTAAGCTGCTTGCGAAGGGCGTCGATACGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACTAACCTGTAAAGAAGTCGAGATTGAGTGGCATGGCACCAATAAAAAGTGGTCt                                                                                   <  2:220659/195‑1 (MQ=255)                                                                                                                                                                 ggCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAAc                   >  2:228952/1‑251 (MQ=255)                                                                                                                                                                            tttATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCt        >  1:171948/1‑251 (MQ=255)                                                                                                                                                                                  cggAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGACGCGTAAGCTGCTTGCGAATGGCGGCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTa  <  2:171948/249‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                  GGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAAAAGTGGTCTGAAGATCACGCCCATCGTGTGCTAAAAAGTCTTGAAGATAATCTTTTTGCAGCGCTTGGTGAACGTAATATCGCTGAGTTA  >  NC_000913/4496667‑4497094

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