Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,426,929	C→T	C165Y (TGC→TAC)	yjfZ ←	DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,426,929	0	C	T	100.0%	150.6 / NA	45	C165Y (TGC→TAC)	yjfZ	DUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/26);  total (19/26)

GCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGA  >  NC_000913/4426686‑4427170                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  gCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAAGCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAg                                                                                                                                                                                                                                             <  2:269362/250‑1 (MQ=255)    ggCGTTAGATACATACGAGCTCAACTGATGTCTCGCTCCTGCGGCACTCTCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACTAATTTTCAAATCATGTGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGTCATTAGAATCGCCTGTAACAGTGCGAAGTTTGTCGATATAAAGATCGATTGTTGGAAGTTCAGATGACCAGTAATTAGGCAt                                                                                                                                                                                                                                         <  1:583807/251‑1 (MQ=255)    ggCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGc                                                                                                                                                                                                                                           <  2:180510/249‑1 (MQ=255)      cGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCAttt                                                                                                                                                                                                                                       >  2:262751/1‑251 (MQ=255)                  acGAGCGCAACTGATGTCACGCTCCTGCGGCACTCGCTCGAAACCCAAACGCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATGGTTGGTAGTTCAGATGACCAGTAATTAGGCAGCTGGGCAGGAAggg                                                                                                                                                                                                                           <  2:262026/251‑1 (MQ=255)                   cGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGt                                                                                                                                                                                                                          >  2:140800/1‑251 (MQ=255)                                   ctcGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATc                                                                                                                                                                                                          >  1:297053/1‑251 (MQ=255)                                      gCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACATCGTATCCTATCGtt                                                                                                                                                                                                       <  2:253723/251‑1 (MQ=255)                                      gCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGtt                                                                                                                                                                                                       >  1:295805/1‑251 (MQ=255)                                               gCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGc                                                                                                                                                                                                                                           >  2:398210/1‑206 (MQ=255)                                               gCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGc                                                                                                                                                                                                                                           <  1:398210/206‑1 (MQ=255)                                                     gctcgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCt                                                                                                                                                                                        <  1:452819/251‑1 (MQ=255)                                                            acccaaacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTAc                                                                                                                                                                                 <  2:299424/251‑1 (MQ=255)                                                                  acccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAAc                                                                                                                                                                           >  1:478752/1‑251 (MQ=255)                                                                      gAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAgg                                                                                                                                                                       >  1:59724/1‑251 (MQ=255)                                                                                       aCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCTGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgt                                                                                                                                                      <  1:449628/251‑1 (MQ=255)                                                                                            aTAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGAGTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCttgtatt                                                                                                                                                 <  2:233995/251‑1 (MQ=255)                                                                                                    cGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGACACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTACTTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGACGGTATAGTCTGGAGTACCATCACCAGGTCCCAAAACCACTTCGTGTATTGTGTAttt                                                                                                                                         >  1:373143/1‑251 (MQ=255)                                                                                                               tttCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGt                                                                                                                               <  1:27291/250‑1 (MQ=255)                                                                                                                   aaaTCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGa                                                                                                                           >  2:186804/1‑250 (MQ=255)                                                                                                                          ggggTAAAGCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTTGTCGATATAAAGAGCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTc                                                                                                                    <  2:466804/250‑1 (MQ=255)                                                                                                                              tAAATCAGAATACGCGTTTAGCACTGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACAt                                                                                                               >  2:206581/1‑251 (MQ=255)                                                                                                                                 aTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGAATTGTGTATTTACTGTTGGGTTCGATGGATTCGACATCtt                                                                                                            >  1:330064/1‑251 (MQ=255)                                                                                                                                     aaaTACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaa                                                                                                        <  2:179289/250‑1 (MQ=255)                                                                                                                                      aaTACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCTTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaat                                                                                                       <  1:566500/251‑1 (MQ=255)                                                                                                                                       aTACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGTAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGTCAGCTTGCaata                                                                                                      <  1:488246/251‑1 (MQ=255)                                                                                                                                            cgTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTCATTAGGCATTTGGGCAGGAAGGGCTACAACGTATCATATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAatat                                                                                                 >  1:80355/1‑251 (MQ=255)                                                                                                                                                 aGCAATGCAACAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTTCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTGCGAAGGATTCGACATCTTGCAATAAATATTTCAg                                                                                            <  2:477657/251‑1 (MQ=255)                                                                                                                                                     aTGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGCTCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCa                                                                                                         >  1:360453/1‑234 (MQ=255)                                                                                                                                                     aTGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGGATTTACTGTGGGGTTCGAAGGATGCGACAGCTTGCa                                                                                                         <  2:360453/234‑1 (MQ=255)                                                                                                                                                               gatgaATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAAcc                                                                              >  1:542253/1‑251 (MQ=255)                                                                                                                                                                      ccGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCTATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCGCTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTTCAATAAATATTTCAGTTCTGCGGTTAACCTCGCAcc                                                                       <  2:404576/251‑1 (MQ=255)                                                                                                                                                                       cGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGGCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATATCAGTTTGGCGGTTAACCTCGCAccc                                                                      <  2:59724/251‑1 (MQ=255)                                                                                                                                                                             cATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTa                                                                 >  2:411509/1‑250 (MQ=255)                                                                                                                                                                                      cGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTg                                                       <  2:584438/251‑1 (MQ=255)                                                                                                                                                                                        ccTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTa                                                     <  1:186804/251‑1 (MQ=255)                                                                                                                                                                                             aaCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCt                                                >  2:113143/1‑251 (MQ=255)                                                                                                                                                                                              aCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTg                                               <  1:411509/251‑1 (MQ=255)                                                                                                                                                                                               cAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTgg                                              >  2:36058/1‑251 (MQ=255)                                                                                                                                                                                                            gTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAAtattat                                 <  1:59865/251‑1 (MQ=255)                                                                                                                                                                                                               gATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCt                               >  1:539134/1‑250 (MQ=255)                                                                                                                                                                                                                  ataAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATg                           <  2:295805/251‑1 (MQ=255)                                                                                                                                                                                                                            aTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTa                                                                 <  1:182726/203‑1 (MQ=255)                                                                                                                                                                                                                                     aGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGTGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTtccatc         <  1:61884/250‑1 (MQ=255)                                                                                                                                                                                                                                            aTGACCAGTAATTAGGCATTTGTGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGa  >  2:441265/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  GCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGA  >  NC_000913/4426686‑4427170

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