Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	738,853	G→A	*254* (TAG→TAA)	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,853	0	G	A	100.0%	135.7 / NA	42	*254* (TAG→TAA)	ybfD	H repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/22);  total (20/22)

GAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCC  >  NC_000913/738727‑738978                                                                                                                               |                                                                                                                              gAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTCAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                   <  2:419602/139‑1 (MQ=255)          aaaTAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGaa                                                                                                          <  1:58311/139‑1 (MQ=255)          aaaTAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGaa                                                                                                          <  1:4624/139‑1 (MQ=255)                tCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                    >  2:721146/1‑139 (MQ=255)                            cAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTcc                                                                                        >  2:10954/1‑139 (MQ=255)                               ttACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTg                                                                                     >  2:535908/1‑139 (MQ=255)                                 aCGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   <  2:161307/139‑1 (MQ=255)                                    cAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTa                                                                                                             <  1:406405/110‑1 (MQ=255)                                    cAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTa                                                                                                             >  2:406405/1‑110 (MQ=255)                                    cAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCACATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCgggg                                                                                >  1:381402/1‑139 (MQ=255)                                              aaaGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCttt                                                                      <  2:525119/139‑1 (MQ=255)                                                  agTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                   >  2:544628/1‑106 (MQ=255)                                                  agTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                   <  1:544628/106‑1 (MQ=255)                                                         gCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           <  1:412368/139‑1 (MQ=255)                                                         gCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           <  1:721146/139‑1 (MQ=255)                                                          cAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                          <  2:177738/139‑1 (MQ=255)                                                                   aaCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                 <  1:55703/139‑1 (MQ=255)                                                                    aCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaac                                                >  2:185735/1‑139 (MQ=255)                                                                        gTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGa                                            <  1:286552/139‑1 (MQ=255)                                                                        gTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGa                                            <  1:61606/139‑1 (MQ=255)                                                                         tCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGAt                                           >  2:374662/1‑139 (MQ=255)                                                                                   gTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTcaca                                 <  2:493399/139‑1 (MQ=255)                                                                                     ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAg                                                                                           >  2:22136/1‑79 (MQ=255)                                                                                     ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAg                                                                                           <  1:22136/79‑1 (MQ=255)                                                                                      tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAat                              >  1:147333/1‑139 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                  >  1:56605/1‑111 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                  <  2:56605/111‑1 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtc                      >  2:235299/1‑139 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtc                      >  2:280178/1‑139 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtc                      >  2:597848/1‑139 (MQ=255)                                                                                              cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtc                      >  2:567910/1‑139 (MQ=255)                                                                                                     gaGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGt               >  2:378698/1‑139 (MQ=255)                                                                                                        cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           <  2:320961/92‑1 (MQ=255)                                                                                                        cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           >  1:320961/1‑92 (MQ=255)                                                                                                          tATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTaca          <  1:517637/139‑1 (MQ=255)                                                                                                           aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                  <  1:312769/98‑1 (MQ=255)                                                                                                           aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                  >  2:312769/1‑98 (MQ=255)                                                                                                            ttGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACt        <  1:374662/139‑1 (MQ=255)                                                                                                              gATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTa      <  1:535908/139‑1 (MQ=255)                                                                                                              gATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTa      <  1:10954/139‑1 (MQ=255)                                                                                                               aTTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAACCTACCTGACGTCAGTCTTTGCGGGGAGAAGGCTTTCGTAAATTTGTCCTGTTACACAAGGATTAATTTCCCAAATATCATTTTTAAACGTTTACACTTAc     >  1:289010/1‑139 (MQ=255)                                                                                                                  tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTcc  >  1:56934/1‑139 (MQ=255)                                                                                                                               |                                                                                                                              GAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCC  >  NC_000913/738727‑738978

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