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%	174.6 / NA	62	F50F (TTT→TTC)	ybfD	H repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/28);  total (34/28)

TATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTT  >  NC_000913/738124‑738376                                                                                                                      |                                                                                                                                        tattatCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGAtttt                                                                                                                    >  2:435745/1‑139 (MQ=21) tattatCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGAtttt                                                                                                                    >  1:231287/1‑139 (MQ=25)   ttatCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGAtttttt                                                                                                                  >  1:666285/1‑139 (MQ=25)    tatCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTg                                                                                                                 <  2:760376/139‑1 (MQ=21)      tCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGaa                                                                                                               >  1:223627/1‑139 (MQ=21)        cccGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGc                                                                                                             >  2:570237/1‑139 (MQ=25)           gATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAat                                                                                                          >  1:1062670/1‑139 (MQ=37)           gATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAat                                                                                                          >  2:286402/1‑139 (MQ=21)            aTTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAata                                                                                                         <  2:666285/139‑1 (MQ=21)                  gACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGa                                                                                                   <  2:1062670/139‑1 (MQ=25)                       gCCTGGAAAGTGGAACATAAATTATCGGATATTCTAGTGTTGACTATGTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGCGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTg                                                                                                    >  1:786533/1‑133 (MQ=255)                       gCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTg                                                                                                    <  2:786533/133‑1 (MQ=21)                         cTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGaa                                                                                            <  1:570237/139‑1 (MQ=21)                                 tGGAACATAAATTATCGGATATTCTCCTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  1:152645/1‑139 (MQ=17)                                 tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAg                                                                                                              >  1:357009/1‑113 (MQ=17)                                 tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAg                                                                                                              <  2:357009/113‑1 (MQ=17)                                 tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  1:580220/1‑139 (MQ=17)                                 tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  2:123016/1‑139 (MQ=17)                                 tGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                    >  1:227069/1‑139 (MQ=18)                                  ggAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAtt                                                                                   >  1:12433/1‑139 (MQ=18)                                      cATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTg                                                                               <  2:1190007/139‑1 (MQ=17)                                        tAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGtt                                                                             <  1:1008507/139‑1 (MQ=18)                                         aaaTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTc                                                                            >  1:934391/1‑139 (MQ=17)                                               tCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATa                                                                      >  2:1228355/1‑139 (MQ=14)                                                cGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATAc                                                                     >  1:575114/1‑139 (MQ=17)                                                  gATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCa                                                                   >  1:1080944/1‑139 (MQ=17)                                                  gATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCa                                                                   >  1:466170/1‑139 (MQ=14)                                                        cTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             <  2:12433/139‑1 (MQ=18)                                                        cTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                                             <  2:263322/139‑1 (MQ=17)                                                          aCTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCaga                                                           <  1:435745/139‑1 (MQ=17)                                                             gTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGtt                                                        <  2:227069/139‑1 (MQ=17)                                                              ttGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTg                                                       <  2:1080944/139‑1 (MQ=18)                                                               tGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGt                                                      <  2:309456/139‑1 (MQ=18)                                                                 aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGc                                                               <  2:317407/128‑1 (MQ=18)                                                                 aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGc                                                               >  1:317407/1‑128 (MQ=18)                                                                 aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTAt                                                    >  1:783478/1‑139 (MQ=18)                                                                 aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTAt                                                    >  2:1138129/1‑139 (MQ=17)                                                                 aCTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTAt                                                    >  1:1068672/1‑139 (MQ=17)                                                                        gtgCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTAt                                             <  2:523312/139‑1 (MQ=18)                                                                          gCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCt                                           >  2:1081345/1‑139 (MQ=18)                                                                             gTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGAGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTc                                        <  2:152645/139‑1 (MQ=21)                                                                                        tGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAAttt                             >  2:151026/1‑139 (MQ=37)                                                                                        tGCAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAAttt                             >  1:650254/1‑139 (MQ=255)                                                                                          cAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCa                           >  2:963837/1‑139 (MQ=255)                                                                                          cAGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCa                           >  1:476148/1‑139 (MQ=37)                                                                                           aGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCAt                          <  2:189446/139‑1 (MQ=255)                                                                                           aGAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCAt                          <  1:1138129/139‑1 (MQ=38)                                                                                            gAAGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCTTGCGAAATTTCATg                         <  2:580220/139‑1 (MQ=255)                                                                                             aaGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATga                        <  2:213062/139‑1 (MQ=255)                                                                                             aaGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATga                        <  1:123016/139‑1 (MQ=255)                                                                                             aaGGCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATga                        <  1:963837/139‑1 (MQ=255)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTg                                                       <  2:518808/106‑1 (MQ=14)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTg                                                       >  1:518808/1‑106 (MQ=17)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTAt                                                    <  1:879394/109‑1 (MQ=17)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATgaga                      <  2:575114/139‑1 (MQ=255)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATgaga                      <  2:934391/139‑1 (MQ=255)                                                                                               ggCTGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTATTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTAt                                                    >  2:879394/1‑109 (MQ=9)                                                                                                  tGGGAAGATATAGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCt                   >  2:458932/1‑139 (MQ=255)                                                                                                           ataGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACt          >  1:777476/1‑139 (MQ=255)                                                                                                           ataGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTCGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATga                        <  2:442673/125‑1 (MQ=255)                                                                                                           ataGAGGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTCGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATga                        >  1:442673/1‑125 (MQ=255)                                                                                                              gagGATTTCGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGa       >  2:629287/1‑139 (MQ=255)                                                                                                                   tttcGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACt          <  1:586188/131‑1 (MQ=255)                                                                                                                   tttcGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCtt  >  2:601931/1‑139 (MQ=255)                                                                                                                      |                                                                                                                                        TATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTATTAACTGGATGCTT  >  NC_000913/738124‑738376

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