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%	132.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)

GCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCC  >  NC_000913/738724‑738978                                                                                                                                  |                                                                                                                              gCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAg                                                                                                                      <  1:1020517/139‑1 (MQ=255)    gAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                   <  2:208324/139‑1 (MQ=255)                   tCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                    >  1:942018/1‑139 (MQ=255)                              aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCt                                                                                                                    >  2:266322/1‑112 (MQ=255)                              aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCt                                                                                                                    <  1:266322/112‑1 (MQ=255)                              aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTc                                                                                         >  2:298946/1‑139 (MQ=255)                              aCAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTc                                                                                         >  2:82732/1‑139 (MQ=255)                                    aCGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   <  1:82732/139‑1 (MQ=255)                                              tGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGc                                                                         <  2:942018/139‑1 (MQ=255)                                                               gagaAGACACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCt                                                                                                                    <  2:742302/79‑1 (MQ=255)                                                               gagaAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCt                                                                                                                    >  1:742302/1‑79 (MQ=255)                                                                     aaaCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTaca                                                  <  1:298946/139‑1 (MQ=255)                                                                                cATATTGTTTGCGATGTCCCTGGTGAGCTTATTGATTTCACGTTTGAATTAAAAGGGCAGCTATGGTTAGAAAc                                                                                                        <  1:792390/74‑1 (MQ=18)                                                                                cATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAAc                                                                                                        >  2:792390/1‑74 (MQ=255)                                                                                      gTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           <  2:840090/113‑1 (MQ=255)                                                                                      gTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           >  1:840090/1‑113 (MQ=255)                                                                                       tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaa                                >  1:1055987/1‑139 (MQ=255)                                                                                        ttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                               >  1:1066014/1‑139 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                         <  1:82789/48‑1 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                         >  2:82789/1‑48 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGc                                                                                                                        >  1:265289/1‑49 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGc                                                                                                                        <  2:265289/49‑1 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                          <  1:392279/111‑1 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                          >  2:392279/1‑111 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTc                                    <  1:577152/133‑1 (MQ=255)                                                                                         tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTc                                    >  2:577152/1‑133 (MQ=255)                                                                                          gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   <  1:345479/85‑1 (MQ=255)                                                                                          gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   >  2:345479/1‑85 (MQ=255)                                                                                             aTGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCa                          >  1:471304/1‑139 (MQ=255)                                                                                                    tgatgaGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGCGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCa                   <  2:1120819/139‑1 (MQ=255)                                                                                                          gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   <  2:537033/69‑1 (MQ=255)                                                                                                          gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                   >  1:537033/1‑69 (MQ=255)                                                                                                          gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCt             <  2:471304/139‑1 (MQ=255)                                                                                                              aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                    >  1:950298/1‑48 (MQ=255)                                                                                                              aTTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTAcc                                                                                                    <  2:950298/48‑1 (MQ=255)                                                                                                               ttGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACt        <  2:1055987/139‑1 (MQ=255)                                                                                                                  aTTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTAc     <  2:499330/139‑1 (MQ=255)                                                                                                                     tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTcc  >  2:628002/1‑139 (MQ=255)                                                                                                                        cGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTaca          >  2:1235273/1‑128 (MQ=255)                                                                                                                        cGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTaca          <  1:1235273/128‑1 (MQ=255)                                                                                                                            tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           <  2:320462/75‑1 (MQ=255)                                                                                                                            tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                           >  1:320462/1‑75 (MQ=255)                                                                                                                                  |                                                                                                                              GCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCC  >  NC_000913/738724‑738978

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