Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,791,057	C→T	G63D (GGT→GAT)	selO ←	protein adenylyltransferase SelO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,791,057	0	C	T	100.0%	124.5 / NA	38	G63D (GGT→GAT)	selO	protein adenylyltransferase SelO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/22);  total (16/22)

GTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGCCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATACCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCA  >  NC_000913/1790811‑1791293                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             gtcgtcGGAATGCCCAGATAATGCATAGCCTCACTGGCGAGACTTTCTCGTATCGTCGAAAGAACAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTAAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGCCCACTGTAAACCTGCGCCAATGGTGACATATCAgg                                                                                                                                                                                                                                          >  2:277809/1‑251 (MQ=255)         aaTGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAggg                                                                                                                                                                                                                                  >  1:231504/1‑251 (MQ=255)             cccAGATAATGCATCGCCTCACTGGCGAGCCTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTc                                                                                                                                                                                                                              <  2:346994/251‑1 (MQ=255)              cctgATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCAGTCGCGAATAAGGCGTCAGGCCAGCACCTGTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTcg                                                                                                                                                                                                                             <  2:231504/248‑1 (MQ=255)               cAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTcgc                                                                                                                                                                                                                            >  2:476239/1‑251 (MQ=255)                                cACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATACCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTgcg                                                                                                                                                                                                           >  1:10366/1‑251 (MQ=255)                                       gagaCTTTCTCGTATCGTCGAACGTACAACCGCCCGTCCAGCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCCTTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACTTGCGCCAGTGGTGACATATCAGGCAGTAGAGCTTCGCCGCCCCCGACCCCTGCGCCAtttt                                                                                                                                                                                                    >  2:435312/1‑251 (MQ=255)                                        agaatttCTCGTATCGTGGAACGTAACACCGCCCGTCCATTCCCCACTCGCGAATAAGGCGTCAGGCCCGCACCATTCATATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCCCCCAGTTGGCCCGCCCAGACGCCGAACTGATGACTACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCAttttt                                                                                                                                                                                                   <  2:426934/247‑1 (MQ=255)                                                    aTCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGcgacg                                                                                                                                                                                        <  1:449255/250‑1 (MQ=255)                                                     tcgtcgAACGTAAAACCGCCCGTCCATCACCAATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGAt                                                                                                                                                                                      >  1:430882/1‑251 (MQ=255)                                                                aaaaCCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTc                                                                                                                                                                             >  1:333398/1‑249 (MQ=255)                                                                      gCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGtt                                                                                                                                                                     <  1:94435/251‑1 (MQ=255)                                                                       cccGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTa                                                                                                                                                                    <  1:476239/251‑1 (MQ=255)                                                                          gTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGcc                                                                                                                                                                 <  1:276523/251‑1 (MQ=255)                                                                            ccATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCa                                                                                                                                                                <  1:505774/250‑1 (MQ=255)                                                                            ccATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTACTGTCATCAGCAAGCAGTTGTTCGCCTAGTAAAATGCCGCGCCCATCTCACAGTTGGCCCGCCCATACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCa                                                                                                                                                                <  1:347422/250‑1 (MQ=255)                                                                                cACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTAGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGACAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCCGACACCTGCGCCATTCTTAAACAGCGACGAAGGAATACTCAGCGTGTTAGCCAGCTCg                                                                                                                                                           >  2:148784/1‑251 (MQ=255)                                                                                    cATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCAGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGAGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGtat                                                                                                                                                       >  2:225941/1‑251 (MQ=255)                                                                                          cgAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCATGCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGGGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGGATTATGcc                                                                                                                                                 <  2:289479/251‑1 (MQ=255)                                                                                             aTAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGGCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCaaa                                                                                                                                              <  2:19365/251‑1 (MQ=255)                                                                                              tAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGAAAGCAGTTGTTCGCCGAGTCAACTGCCGCGCCCATCACCCAGTTAGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCCTTTTTAAACAGCGACGATGAAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAAt                                                                                                                                             >  1:295649/1‑251 (MQ=255)                                                                                                  gCGTCAGGCCAGCGCCTTTCAGATGCCAGTCCATTTTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCTCCCAGACGCCTAACTGATGACCACTGTAAACCTGCGCCAGTGGTGATATATCAGGCAGTATGGCTTCGCCGCCCCAGACTCCTGCGCCATTTCTAAACAGCGACTATGGAATACTCAGCGTGTTAGCCAGTGCGGTATTATGCCAAATCAGc                                                                                                                                         <  2:295382/251‑1 (MQ=255)                                                                                                             gCACCTTTCAGATGCCAGTCCATTGTAGTGTCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCcttatt                                                                                                                              >  2:493776/1‑251 (MQ=255)                                                                                                                            cAGTCCATTGTAGTATCATCAGCAAGCAGTTGTTCGCCAGGCAAACTGCAACGACCATCACACAGTTGGCCCGCCCAGACGCCGTACTGATGACCACTGTAAACCTACTCCAGTGGGGACATATCAGACAGTAGGCCTTCGCCGCCCCAGACACCTGCCCCATTTTTAGACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGGATTATGCCAACTCAGCCGGGCATTATTTAACGGCGTAggggc                                                                                                               >  2:517066/1‑250 (MQ=255)                                                                                                                                          gtcATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAgg                                                                                                 <  2:400587/249‑1 (MQ=255)                                                                                                                                                       aGTTGTTCGACGAGTAAAATGCCGCGCCCATCGCCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGTGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATTGAATACTCAGCGTTTTAGCCAGTTCGGTATTATGCCAAATCATGCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTGTTTGGCAATTc                                                                                    <  2:390546/251‑1 (MQ=255)                                                                                                                                                             tCGCCGAGTAAAATGCCGCGCCCATCACCCAGTGGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATcgc                                                                               <  2:178818/250‑1 (MQ=255)                                                                                                                                                                      aaaaTGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCggg                                                                      <  2:333398/250‑1 (MQ=255)                                                                                                                                                                         aggCCGCGCCCATCCACCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAGCCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAAc                                                                  <  2:88627/249‑1 (MQ=255)                                                                                                                                                                              gcgcCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACaaa                                                               >  1:88419/1‑249 (MQ=255)                                                                                                                                                                               cgcCCATCCCCCAGTTGGCCCGCCCAGCCGCCGAACTGATGACCCCTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGa                                                            <  2:528830/251‑1 (MQ=255)                                                                                                                                                                                        cccAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACTATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAGATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTTCTGTAGAGGTTTCTGGCAATTCATCGCGCCAGCTGTTAACAAAAGACAGTGTCAt                                                   <  2:287550/251‑1 (MQ=255)                                                                                                                                                                                                        aGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGCCATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAACAGACAGGGTCATAGATCCTCCTGTTTg                                    >  2:52525/1‑250 (MQ=255)                                                                                                                                                                                                         gACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGAt                                  <  1:202888/251‑1 (MQ=255)                                                                                                                                                                                                                       gaCCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGtt                     <  2:88419/250‑1 (MQ=255)                                                                                                                                                                                                                                  aCCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGCCGTTTACTCTCGTCtt          >  2:16009/1‑250 (MQ=255)                                                                                                                                                                                                                                         ccAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCa  <  1:171922/251‑1 (MQ=255)                                                                                                                                                                                                                                          cAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACAcc   >  1:204234/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             GTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGCCATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATACCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCA  >  NC_000913/1790811‑1791293

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