Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	738,241	T→C	F50F (TTT→TTC)	ybfD →	H repeat‑associated putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	738,241	0	T	C	100.0%	100.4 / NA	38	F50F (TTT→TTC)	ybfD	H repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/17);  total (21/17)

TTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTT  >  NC_000913/738126‑738376                                                                                                                    |                                                                                                                                        ttatCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTATTGATGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACAGATCCCGAtttgtt                                                                                                                  >  1:403236/1‑139 (MQ=11)     ccccTATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAg                                                                                                              <  1:546468/139‑1 (MQ=21)      cccGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAGGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGc                                                                                                             <  2:411730/139‑1 (MQ=255)         gATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGa                                                                                                                        <  2:439961/125‑1 (MQ=25)         gATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGa                                                                                                                        >  1:439961/1‑125 (MQ=25)         gATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAat                                                                                                          >  2:450659/1‑139 (MQ=21)                   aaGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGAttt                                                                                                >  1:126974/1‑139 (MQ=25)                     gCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAg                                                                                                              >  1:262548/1‑123 (MQ=25)                     gCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAg                                                                                                              <  2:262548/123‑1 (MQ=25)                               tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  1:508723/1‑139 (MQ=17)                               tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  2:32281/1‑139 (MQ=17)                                         aTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAaca                                                                                                                                <  1:366116/85‑1 (MQ=14)                                         aTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAaca                                                                                                                                >  2:366116/1‑85 (MQ=14)                                              cGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGCGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATAc                                                                     <  2:403236/139‑1 (MQ=17)                                                     tCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGcc                                                              >  2:82642/1‑139 (MQ=17)                                                      cTAGTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             <  2:425974/139‑1 (MQ=18)                                                            ttGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTg                                                       >  2:4341/1‑139 (MQ=17)                                                              gACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTa                                                     >  2:560125/1‑139 (MQ=17)                                                              gACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTa                                                     >  2:453782/1‑139 (MQ=18)                                                              gACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTa                                                     <  2:323651/139‑1 (MQ=17)                                                              gACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTa                                                     >  1:606672/1‑139 (MQ=17)                                                               aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             <  1:200861/130‑1 (MQ=17)                                                               aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             >  2:200861/1‑130 (MQ=14)                                                                cTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATc                                                   <  1:4341/139‑1 (MQ=17)                                                                 tATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATcc                                                  <  1:453782/139‑1 (MQ=14)                                                                      gtgCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             <  1:399742/123‑1 (MQ=18)                                                                      gtgCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             >  2:399742/1‑123 (MQ=17)                                                                         ccGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTg                                          >  1:588252/1‑139 (MQ=21)                                                                         ccGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTg                                          >  2:377397/1‑139 (MQ=18)                                                                               tttCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTAt                                             >  1:161926/1‑130 (MQ=18)                                                                               tttCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTAt                                             <  2:161926/130‑1 (MQ=18)                                                                                        cAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCa                           >  2:129220/1‑139 (MQ=38)                                                                                         aGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCAt                          <  1:32281/139‑1 (MQ=38)                                                                                             ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATgaga                      <  1:450659/139‑1 (MQ=255)                                                                                              gCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGagag                     <  2:588252/139‑1 (MQ=255)                                                                                                tGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCt                   >  2:279153/1‑139 (MQ=255)                                                                                                tGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCt                   <  1:279153/139‑1 (MQ=255)                                                                                                             agGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGAt      >  2:178788/1‑139 (MQ=255)                                                                                                                 tttcGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCtt  >  2:575875/1‑139 (MQ=255)                                                                                                                    |                                                                                                                                        TTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTT  >  NC_000913/738126‑738376

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