Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	738,456	G→A	R122H (CGT→CAT)	ybfD →	H repeat‑associated putative transposase YbfD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	738,456	0	G	A	100.0%	60.4 / NA	19	R122H (CGT→CAT)	ybfD	H repeat‑associated putative transposase YbfD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/11);  total (8/11)

TGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTG  >  NC_000913/738241‑738674                                                                                                                                                                                                                        |                                                                                                                                                                                                                           tGGGGAAACACATCCCTATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCAGTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTGCAGATGATAAAGACGGCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTca                                                                                                                                                                                          <  1:246247/250‑1 (MQ=255)             tcaagATTTTTTCAAGCAATATGGTGATTTTGAAAATGGCCTTCCTGTTCACGATCCCATTGCCAGAGTTGTATCATGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAACGACGTCATCGCACTTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTc                                                                                                                                                                             >  2:193349/5‑251 (MQ=17)                  aTTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCGGCAGATGATAAAGACTTCATCGCAATTGATTGAAAAATACACCGGCATTCGTATGACAATAGGCGCCATAAGGGAGCGAGTCATGTCAGTAGTGCGTTCTCAACAATGCACAGTCTGGTc                                                                                                                                                                        <  1:213621/251‑1 (MQ=17)                         gAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGa                                                                                                                                                                  <  1:132705/250‑1 (MQ=255)                                                                cGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCt                                                                                                                           >  2:120990/1‑250 (MQ=255)                                                                   tACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATcc                                                                                                                       >  1:102151/1‑251 (MQ=255)                                                                     ccATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCt                                                                                                                     >  2:236123/1‑251 (MQ=255)                                                                                         tGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTACCTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACAGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGACCTTCTTAACATGCTgg                                                                                                  >  2:367899/1‑250 (MQ=255)                                                                                                               aTGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCaaa                                                                           <  1:38243/251‑1 (MQ=255)                                                                                                                 gagaGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCaaaa                                                                          <  1:335015/250‑1 (MQ=255)                                                                                                                 gagaGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCaaaa                                                                          <  1:372859/250‑1 (MQ=255)                                                                                                                     gCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCGGGATATTAAAGGAAAAATCATCAAAACcgat                                                                     <  2:212932/251‑1 (MQ=255)                                                                                                                         tattaACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACcgatgcg                                                                  >  1:66748/1‑250 (MQ=255)                                                                                                                             aaCTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATgg                                                              <  1:327291/250‑1 (MQ=255)                                                                                                                               cTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGtt                                                           >  2:200722/1‑251 (MQ=255)                                                                                                                                                               aaGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACaaaaa                           <  1:111837/251‑1 (MQ=255)                                                                                                                                                                aGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAAc                          >  2:106949/1‑251 (MQ=255)                                                                                                                                                                                gATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGAttattt          <  2:12083/251‑1 (MQ=255)                                                                                                                                                                                        aaTACACCGGCATTCTTATGACAAGAGTCGCCATAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATTAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTg  <  2:190553/251‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                           TGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTTGACTACCATTCTTCAGATGATAAAGACGTCATCGCAATTGATGGAAAAATACACCGGCATTCTTATGACAAGAGTCGCCGTAAGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCAAAACCGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTG  >  NC_000913/738241‑738674

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