Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,379,942	A→G	R192R (CGT→CGC)	frdB ←	fumarate reductase (anaerobic), Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,379,942	0	A	G	100.0%	201.5 / NA	56	R192R (CGT→CGC)	frdB	fumarate reductase (anaerobic), Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/30);  total (26/30)

CAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATACGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGT  >  NC_000913/4379795‑4380087                                                                                                                                                    |                                                                                                                                                  cAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATgcgc                                                                                                                                                >  2:726678/1‑151 (MQ=255)     gTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCct                                                                                                                                            <  1:547268/151‑1 (MQ=255)     gTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCct                                                                                                                                            >  1:679574/1‑151 (MQ=255)        gcgATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCcttct                                                                                                                                         >  1:267732/1‑151 (MQ=255)        gcgATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCcttct                                                                                                                                         <  1:575169/151‑1 (MQ=255)         cgATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCcttctt                                                                                                                                        >  1:592069/1‑151 (MQ=255)               gAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTg                                                                                                                                  <  2:3934/151‑1 (MQ=255)                aaaGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTgg                                                                                                                                 >  1:255377/1‑151 (MQ=255)                  aGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTc                                                                                                                               <  1:348868/151‑1 (MQ=255)                      tttCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCgg                                                                                                                           >  1:47487/1‑151 (MQ=255)                       ttCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTcgcg                                                                                                                            <  1:659628/149‑1 (MQ=255)                            aCTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCt                                                                                                                     <  2:297258/151‑1 (MQ=255)                            aCTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCt                                                                                                                     <  1:647562/151‑1 (MQ=255)                            aCTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCt                                                                                                                     <  1:415655/151‑1 (MQ=255)                              tttCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTc                                                                                                                   >  1:266993/1‑151 (MQ=255)                                tCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCg                                                                                                                  <  1:449454/150‑1 (MQ=255)                                  tACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGtt                                                                                                                >  2:709325/1‑150 (MQ=255)                                     tttGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTtata                                                                                                             >  1:686319/1‑150 (MQ=255)                                      ttGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTtata                                                                                                             <  1:323646/149‑1 (MQ=255)                                       tGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACg                                                                                                          >  1:134976/1‑151 (MQ=255)                                        gCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACg                                                                                                          >  1:458240/1‑150 (MQ=255)                                         cccTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGAt                                                                                                        <  2:851562/151‑1 (MQ=255)                                             gctgAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATgc                                                                                                      <  1:239302/149‑1 (MQ=255)                                             gctgAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATgc                                                                                                      <  1:326061/149‑1 (MQ=255)                                               tgAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGcc                                                                                                   >  2:628259/1‑150 (MQ=255)                                                gAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATgcg                                                                                                     >  2:580431/1‑147 (MQ=255)                                                  aTGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCg                                                                                               <  2:671124/151‑1 (MQ=255)                                                  aTGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCg                                                                                               <  1:723171/151‑1 (MQ=255)                                                    ggCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGt                                                                                              >  1:847782/1‑150 (MQ=255)                                                    ggCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTa                                                                                             <  1:52610/151‑1 (MQ=255)                                                      ccGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAAt                                                                                           <  2:592107/151‑1 (MQ=255)                                                          aGCCGGATCGACGTGTTTCGGGCAGACTTCGTAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCATACCGTGGTCGCGGCTAACTTCGTTATAACGATGCGCCAGCGTAATGGCa                                                                                       <  1:477210/151‑1 (MQ=255)                                                                          ttCGGGCAGACTTCGGAGCAGTAGCCCACGAGAGTACAGCTCCATACGCCGTTCTGTCTGTCCAACTGCGCCATGCGCTCCTACTTACCGTGGTCGGGGCTATCTTAGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGaa                                                                         <  1:550097/149‑1 (MQ=255)                                                                                cAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGt                                                                 >  1:395972/1‑151 (MQ=255)                                                                                 aGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGtt                                                                >  1:809665/1‑151 (MQ=255)                                                                                     ttCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAgg                                                            >  1:670358/1‑151 (MQ=255)                                                                                          aGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCaaa                                                       <  1:80396/151‑1 (MQ=255)                                                                                              gTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGc                                                   >  2:708827/1‑151 (MQ=255)                                                                                                aGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCgg                                                 >  2:239622/1‑151 (MQ=255)                                                                                                aGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCgg                                                 <  2:267124/151‑1 (MQ=255)                                                                                                    cACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCat                                             <  1:124893/151‑2 (MQ=255)                                                                                                       gAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCAcgcg                                          >  1:351985/1‑151 (MQ=255)                                                                                                        aaaGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACgcgg                                         >  2:538309/1‑151 (MQ=255)                                                                                                         aaGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACgcggc                                        <  1:667207/151‑1 (MQ=255)                                                                                                             aCAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAg                                    <  1:370170/151‑1 (MQ=255)                                                                                                                 cTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGCTGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGcac                                 >  1:525384/1‑150 (MQ=255)                                                                                                                    cATACGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAAc                             <  1:407858/151‑1 (MQ=255)                                                                                                                       aCGCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACca                           <  1:842948/150‑1 (MQ=255)                                                                                                                         gCCGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAg                        <  1:538272/151‑1 (MQ=255)                                                                                                                          ccGTTCTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGt                       >  2:302515/1‑151 (MQ=255)                                                                                                                               cTGGCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATg                  >  1:248844/1‑151 (MQ=255)                                                                                                                                 ggCTGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGAt                   <  2:448592/148‑1 (MQ=255)                                                                                                                                    tGTTCAACTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAAcc             >  2:302200/1‑151 (MQ=255)                                                                                                                                         aaCTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCgg           >  1:298981/1‑148 (MQ=255)                                                                                                                                           cTGCGCCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAAc      >  2:15364/1‑151 (MQ=255)                                                                                                                                             gcgcCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACt     <  2:164436/150‑1 (MQ=255)                                                                                                                                               gcCATGCGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGt  <  2:318724/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGCAGTAGCCCACGAAAGTACAGCTCCATACGCCGTTCTGGCTGTTCAACTGCGCCATACGCTCCTTCTTACCGTGGTCGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGT  >  NC_000913/4379795‑4380087

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