Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,837	G→A	D158N (GAT→AAT)	lpxK →	tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,837	0	G	A	100.0%	108.1 / NA	38	D158N (GAT→AAT)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/23);  total (15/23)

CGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGG  >  NC_000913/968599‑969081                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ggTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTg                                                                                                                                                                                                                                          <  1:200692/250‑1 (MQ=255)                   gggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGc                                                                                                                                                                                                                        <  1:238093/251‑1 (MQ=255)                        atGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGGTGAATTATTAAAGCACCAATGCGCCCGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGGAAAAGCCATTCTGGCGCCACACCCTGCTGTGCCAATCCTCTTAACCGACCACGGTTTACCGACTTCCCGTCCGGCCCGTAATGTGGACATTGTTCTTTTTGATTGTGTGCGTccg                                                                                                                                                                                                                   >  2:283105/1‑249 (MQ=255)                                   gCTGAATCTTATCCGCTGGTAGTGTCGGCAGATACCACAACAGCTCAGGCTGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGGGGAGATCATCGTAACCGACGACGGTGTTCAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATgg                                                                                                                                                                                                        <  1:400876/251‑1 (MQ=255)                                      gAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTTGCGCGGAATGTGGGAATTGTCGTTATTGATGGGGTGCGTCGCTTTTGCAATGGCtg                                                                                                                                                                                                     >  1:86396/1‑251 (MQ=255)                                       aaTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGGTTTATCAACGCACCGATGCGCCTGTTGCGGTTTTTCCCGTTC‑GTTCTGATGCGGTAAAAACCATTCTGGCGCAAAACCCTGATGTGCAGATCATCGTAAACGACGACGGTTTACAACCTTACCGGCTGGCGCGTAATGTGGGAATTGTCGTTATTTATGGGGTGGGTCGCTTTTGGAATGGCtgg                                                                                                                                                                                                    >  1:406227/1‑251 (MQ=255)                                                 gCTGTTATTGTCGGCAGATACCACAACAGCACAGGCTGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcgg                                                                                                                                                                                          <  2:170241/251‑1 (MQ=255)                                                            cGGCAGATACCACAACAGCACAGGCGGGTGATGAACCCGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGCCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTTTCGTTTTTTGTGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCgggggccatg                                                                                                                                                                               >  1:149284/1‑245 (MQ=255)                                                               cAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATAAACGCCCTGATGCGCCTGTTGCGGGTTCTCCCGTTCTGTT‑TGATGCGGTAAAAGCCATTCTTGCGCAACACCCTGATGTGCAGAACATCGTAACCGGCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTATTGTTGGTGTTCGTCGCTTTTGCCATGGCTGGTGGGTGCCGGCgggggcccatcgg                                                                                                                                                                            >  1:16816/1‑242 (MQ=255)                                                                 gATACCACAACAGCACAGGCGGGGGATGAAACTGTGATGATTTATCAAAACACCGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGCTGCGGGAAAAACCATTCTGGCGCCACACCCCGCTGTGCAAGTCATAGTAACCGAAGCCGGTTTACAGCATTACAGACTGGCGCGGAATGTTGAAGTTGTTGGTATTAAAGGTGTGGGGCGCTTTGGGAATGGCTAGTGGGTGCCGgggggggcacagtgtgg                                                                                                                                                                          >  2:342887/1‑240 (MQ=255)                                                                  aTACCACAACAGCACAGGCGGGGGAAGAAACCGTGTTGGTTTATTAACGCCCCGGTGGGCCTGCTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGGAAAAGGCATTCTGGCGCAACAACCCGATGGGCGCAGCATCGGAACCGACGTCGGTTTAAAGCATTACCGACTGGCGCGGAATGGGGAAAATTTCGGTATTGATGGGGGGGGGCGCTTTGGGGATGGCAGGTGGTTGGCGGCgggggccatcacgcgg                                                                                                                                                                         >  1:258324/1‑239 (MQ=255)                                                                    aCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagc                                                                                                                                                                        <  2:86396/250‑1 (MQ=255)                                                                                ccggcggggtatgaACCTGTTTTGATTTATCAACGCACTGCTGCGCCTGTTGCTGTTTCTCCCGTTC‑GTTCTGAGGCGGTAAAATCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGATCGGGTCGCtt                                                                                                                                                           <  2:261689/241‑1 (MQ=255)                                                                                caGGCGGGTGATGAACCTGGGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCGGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCtt                                                                                                                                                           <  2:406227/251‑1 (MQ=255)                                                                                 aGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGAGTGAGCGCGCGGGGCGCttt                                                                                                                                                          >  1:293983/1‑250 (MQ=255)                                                                                    cGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCACTGCgggggggggcggggcggttta                                                                                                                                                         >  2:229696/1‑229 (MQ=255)                                                                                      ggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGGGGGAGCGGGGCGCTTAAAGtg                                                                                                                                                     >  1:367525/1‑250 (MQ=255)                                                                                              aCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGGGTGGGCGGGCGGGGCGCTTAAAGTCGGTTGagg                                                                                                                                             >  1:116047/1‑249 (MQ=255)                                                                                               ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGc                                                                                                                                            <  1:113704/251‑1 (MQ=255)                                                                                                    gtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTa                                                                                                                                        <  1:325989/250‑1 (MQ=255)                                                                                                              cAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATTATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcgg                                                                                                                             <  1:119232/251‑1 (MQ=255)                                                                                                                    aCTGATGCGCCTGTTGGGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGgtgt                                                                                                                       <  2:181609/251‑1 (MQ=255)                                                                                                                    aCTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGGTGCGGTAATCGTCAAcggcggggg                                                                                                                       >  1:412663/1‑248 (MQ=255)                                                                                                                                               gttc‑gttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAAAGCGTGAGCGGGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCCCGCAGCGGGGAAAAcccccttcc                                                                                            >  1:224394/1‑246 (MQ=255)                                                                                                                                                     ttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCCCGCAGCGGGGAAATTCCCCTGCATcttc                                                                                       >  1:162283/1‑249 (MQ=255)                                                                                                                                                                                gcAACTCCCTTATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTa                                                            <  2:224394/251‑1 (MQ=255)                                                                                                                                                                                    cacccgtaggtgCTGATCATCGTAACCGACGACTGTTGACAGTATTACCGTCTGGCTCTTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCACTGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATTCTGTAATCGTCAACGGCTGTGTCCCTCGCAGCGTTGAAATCCCCATGCAGCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCgg                                                        <  1:56484/242‑1 (MQ=255)                                                                                                                                                                                                  gATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGTTGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGGGCGGGGCGCTTAAAATCGGTTGGTGCGGGAATCGTCAACGGCGGGGTCCCCCGCCGCGGGGGAATTCCCCTGCATCTGCTTCCGGGGCCGGGGGGGGATTTTCCTACCGGTACGCGGtgtgggg                                          >  1:47456/1‑248 (MQ=255)                                                                                                                                                                                                   atcatcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGt                                         <  1:229696/251‑1 (MQ=255)                                                                                                                                                                                                     catcGTAGCCGACGACGGTGTACAGCATTACCGTCTGGCGCGGAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGTTGAATTTACGTACCGGTACGCGTTGTGACGtt                                        <  2:103489/250‑1 (MQ=255)                                                                                                                                                                                                      atcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGc                                      <  1:283105/251‑1 (MQ=255)                                                                                                                                                                                                                 ggacgGTTTACAGCATTACCGTCTGGCGCGTAATGGGGAAATTGTCGTTATTGAGGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGGCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAAc                           <  2:149284/250‑1 (MQ=255)                                                                                                                                                                                                                   acgGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCACTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACAt                         <  1:353732/251‑1 (MQ=255)                                                                                                                                                                                                                    cgGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATg                        <  1:89432/251‑1 (MQ=255)                                                                                                                                                                                                                          aCAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGAGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTgg                  <  1:338995/251‑1 (MQ=255)                                                                                                                                                                                                                               aTTACCGTATGGCGCGTAAGGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATg             <  2:213903/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTg   <  1:16290/251‑1 (MQ=255)                                                                                                                                                                                                                                          gcgcGTAATGTGTAATTTGTTGTTAGTGATGTTGTGCGTCGCGTTGGCAATGGCTTGTGGTTGCCGGCGGGGCCATTGCGTGAGCGAGGTGGGCGCTTAAAGTCGGTTGATGCGGTAAACGTCAACGGCGGGGTCCCTCGCAGGGTTGGAATCCCCATGCATGTGCTGCCGGGGCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACAGGTAGTCGCGATGGTGGGGATTgg  <  2:112356/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     CGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTC‑GTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGG  >  NC_000913/968599‑969081

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