Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,590	C→T	V25I (GTC→ATC)	dcuS ←	sensor histidine kinase DcuS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,590	0	C	T	100.0%	141.0 / NA	42	V25I (GTC→ATC)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/21);  total (21/21)

GAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATA  >  NC_000913/4351344‑4351822                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                       <  1:107510/251‑1 (MQ=255) gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                       <  2:421760/251‑1 (MQ=255)       tttGCATATCGGTAACGACAATAAACAGCAGATA‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATAAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGAt                                                                                                                                                                                                                                  >  1:151670/1‑250 (MQ=255)          gcatgacGGTAACGACAATAAACAGCAG‑TCTGTTGCGTTTGCGTCCGGTGTCCGGGTGGGCCTGGATGCTCCTCTCCTGCGGTTTTTGCGGCAAGCCCGGCCGGATGTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTACCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGGACCCCCAATAGCACCGACAACAGTACCGCCCTGATCATTAATATCAc                                                                                                                                                                                                                               <  2:3964/244‑1 (MQ=255)          gCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCTCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCAc                                                                                                                                                                                                                               <  1:60700/250‑1 (MQ=255)             aatCGGTAACGACAATAAACAGCAGTTC‑GTTGCGTTGGCTTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGTACCGAGAACAGTACCGCACTGATCAGTAAGATCACtgtg                                                                                                                                                                                                                           <  2:381563/250‑1 (MQ=255)                        aCAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAttt                                                                                                                                                                                                                <  1:466371/251‑1 (MQ=255)                           aTAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCAt                                                                                                                                                                                                             >  1:213750/1‑251 (MQ=255)                              aaCAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATGTCATCgg                                                                                                                                                                                                          <  1:332705/251‑1 (MQ=255)                                        c‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGt                                                                                                                                                                                                >  1:109957/1‑251 (MQ=255)                                               gtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATg                                                                                                                                                                                          >  1:334517/1‑251 (MQ=255)                                                                 cgATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGt                                                                                                                                                                        >  2:94987/1‑251 (MQ=255)                                                                  gATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtc                                                                                                                                                                       >  1:70117/1‑251 (MQ=255)                                                                   aTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGCGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGttt                                                                                                                                                                      >  2:119843/1‑249 (MQ=255)                                                                               ccACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                                                                                            <  2:221283/201‑1 (MQ=255)                                                                               ccACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                                                                                            >  1:221283/1‑201 (MQ=255)                                                                                cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGCGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCtt                                                                                                                                                         >  2:179058/1‑251 (MQ=255)                                                                                cACTCTCCTGCGGTTTTTTATGCAAGCCGTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCCGTGCCTTGTTGGCTAGCACATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCACTAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGCATGTCTCATCAGCTTCCtt                                                                                                                                                         >  2:120164/1‑251 (MQ=255)                                                                                        tGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGAc                                                                                                                                                  >  2:112129/1‑250 (MQ=255)                                                                                                  tCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATTCGGTAGTGCAATTAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTa                                                                                                                                        <  2:70117/250‑1 (MQ=255)                                                                                                         gCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAt                                                                                                                                <  1:87614/251‑1 (MQ=255)                                                                                                          cccTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACACATTTCTTAAGCATTAtc                                                                                                                               >  2:118325/1‑251 (MQ=255)                                                                                                              gACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtctc                                                                                                                             <  2:1039/249‑1 (MQ=255)                                                                                                                cGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTGGTTGGCTATCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAAACAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgat                                                                                                                         <  1:255892/251‑1 (MQ=255)                                                                                                                        cGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCggg                                                                                                                 <  2:213750/251‑1 (MQ=255)                                                                                                                            gAGTCGGCGAGGGAACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACCGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGCACTCAATTTCATCGGACGTTTGCGGAACATGCGGTAGGGCAGTGAATGTCTCATCAGCTTCCTTGTGTGCCAAATTTCTTAATCATTATCTCTGATGAGGCggtgatt                                                                                                             >  2:255220/1‑246 (MQ=255)                                                                                                                                      gggTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGGAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAg                                                                                                   <  1:291258/251‑1 (MQ=255)                                                                                                                                          aCGCGCCGCTGCCAGTGCCTTGTGGGCTAGCCCATCTCGGGTCATATCACTGATTTGCTAGAAGTAAAGCAGATGCACCACCAATGGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACTGTTGCGTAACATGCTGTAGGGCAATGAATGTCTCATCAGCTTCATTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTaa                                                                                                <  1:35847/250‑1 (MQ=255)                                                                                                                                             cgcCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAAt                                                                                            <  1:162220/251‑1 (MQ=255)                                                                                                                                                    gCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGAGTAAGAATGATTggg                                                                                     >  1:345581/1‑250 (MQ=255)                                                                                                                                                       aGTGCCTTGTTGGCTAGCCCATCTCGTTTCATATCTCTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCTCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCTGTAGGTCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCtat                                                                                  <  2:357608/251‑1 (MQ=255)                                                                                                                                                               gttgGATAGACCATCTCGTGTCATATCACGGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACGCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGa                                                                         <  1:347136/251‑1 (MQ=255)                                                                                                                                                                                 tgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGAACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTCATTCACAGAGAGTAAGAATAATTGGCTATTTAGAGGGTA‑AGACTCTGGCGACggg                                                       >  2:18944/2‑250 (MQ=255)                                                                                                                                                                                       atCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGCATGAGTTGCTCTATTG‑GGGAAAGGACTCTGGCAAAGGGAACTGc                                                 >  2:136380/1‑251 (MQ=255)                                                                                                                                                                                       atCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAAAATAATTGGCTTTATAG‑GGGATGAGACTCTGGCAAAGGGAACTGc                                                 >  2:48893/1‑251 (MQ=255)                                                                                                                                                                                         cACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGAAGAGACTCTGGCAACGGAAACTGCCa                                               <  1:387343/251‑1 (MQ=255)                                                                                                                                                                                                           aaaTCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTGATTCCAAGGGGGTAAGAATGATTGGCTATATTG‑GGGGAGAGACTCTGGCAACGGACACTGGCAGTGCTGTATGAAGATTcc                             >  2:152614/1‑251 (MQ=255)                                                                                                                                                                                                                     cactaccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTTTGGTACTCAATTTCATCGTACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCGATATAG‑GGGAAGAGACTCTGGCAGCGGAAACTGCCAGTGCTGTATGAACATTCCGGGGCTATGc                   <  2:404382/247‑1 (MQ=255)                                                                                                                                                                                                                          ccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTtat               >  1:218945/1‑250 (MQ=255)                                                                                                                                                                                                                              tAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGAACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGAAAGAATGGATGGCTTTATAG‑GGGAAGAGACTCTAGCAACGGAAACTTCCAGTGCTGTGTGAAGATTCCCGGGCTATGGTTATTGCg           >  1:245586/1‑250 (MQ=255)                                                                                                                                                                                                                                 cACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGGAAGAATGATTGGCTATATAG‑GGGGAGGGACTCTGGGAACGGAAAATGACCGTGCTGTATGAAGGTTCCGGGGCTATGCttttaacggta        >  1:395821/1‑241 (MQ=255)                                                                                                                                                                                                                                      agaACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATa  <  1:152614/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          GAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAG‑GGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATA  >  NC_000913/4351344‑4351822

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