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%	238.5 / NA	71	D158N (GAT→AAT)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/33);  total (38/33)

CGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATT  >  NC_000913/968597‑969079                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    cgcgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGACTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTTTGGCGCAACACCCTTGTGTGCAGATCATCGTAACCGACGGCGGGTTACAGCCTTACCGTCTTGGGCGTAATGTgggact                                                                                                                                                                                                                                           >  2:518668/1‑247 (MQ=255) cgcgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGGTTATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGGTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTCACCGACGACGGTTTACAGAATTACCGTCTGGCGGGTAATGTGGAAAt                                                                                                                                                                                                                                           >  1:94093/1‑251 (MQ=255)             gTATCGCGGGGATATGGAGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTATGTTGATTTATCCACGCACTTATGCGCCTGTTGCGGGTTCTCCCGTTCGTTTTGGTGCGGTAAAAGCCCTTATGGCGCAACACCCTAATGTGCCGATCATCGGAACCGACAACGGGTTACAGCATTACCTTCTGGCGCGTAATGTTGGAATTTTCGTTTTTgg                                                                                                                                                                                                                               >  1:31818/1‑250 (MQ=255)                   cgGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTTATTTATCAACGCACTGCTGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGGCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGgtgt                                                                                                                                                                                                                         >  1:452911/1‑251 (MQ=255)                   cgGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGTTTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGgtgg                                                                                                                                                                                                                         >  1:100705/1‑250 (MQ=255)                         tatGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTTTGATGCGGTAAAAGCCATTTTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGCCGCCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCg                                                                                                                                                                                                                   >  1:112335/1‑251 (MQ=255)                             gtggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCttt                                                                                                                                                                                                               <  1:259203/251‑1 (MQ=255)                                  aaGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTAACGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGGGTGGGTCGCTTTGGCa                                                                                                                                                                                                           >  1:626225/1‑250 (MQ=255)                                        gAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTTTCAAAGCACTGATGCGCCTGTTGAGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAACGACGACGGGTTACAGCATTAACGTTTGGCGCGAAATGTGGACATTGTCGTTATTGATGGGGTGCGTCGCTTTGGCAATGGCtg                                                                                                                                                                                                    >  1:517603/1‑251 (MQ=255)                                         aaTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATAAACGCCCTGATGCGCCTGTTGCGGTTTCTCCCCTTCGTTCTGGTGCGGTAAAAACCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGGAACCGAAGACGGGTTTCAGCATTACCGTCTGGCGCGCAATGTGGAAATTGTCGGTATTGATGGTGTGCGGCGCTTTGTCAATTGCtgg                                                                                                                                                                                                   >  2:2940/1‑251 (MQ=255)                                            cTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGCTGAACCTGTGCTGATTTATCAACGCAATGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTCATGTGCAGATCATCGTAACCGACAACGGCTTACAGCACTACCGTCTTGCGCGTAATGTGGAAATTGTCGTTATTGGTGGGTTGCGTCGCTTTGGGAAttgccagggg                                                                                                                                                                                                >  1:404326/1‑242 (MQ=255)                                                      gTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCcggc‑gg                                                                                                                                                                                     >  1:339386/1‑251 (MQ=255)                                                             tCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTTTCAACGCACTGATGCGCCTGTTGCGGTTTATTCCGTTCGTTCTGATGCGGTTAAAGCCCTTATGGGGCAACACCCTGCTGTGCAGATCATAGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTTTTGGTGGGGTGCGTCGGTTTGGCAATGGGTGGGGGTTGCCGGC‑gggggcca                                                                                                                                                                               >  1:726563/1‑246 (MQ=255)                                                                       ccACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTgagcga                                                                                                                                                                    <  1:518668/251‑1 (MQ=255)                                                                        cacaACAGCACAGGCGGGTGATGAACCCGTGTTTATTTATCCACCCACTCATGCGCCTGTTGCGGTTTCTCCCCTTCGTTCTGATGCGGTAAAAAGCATTCTGGCGCAACACCCTGGTGTGTAGGTCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTTTTGATGGGGTGCGTCGCTTTGGCAATGGCTGGTGGTTGTCCGC‑gggggcacagcgggggcg                                                                                                                                                                     >  1:32386/1‑235 (MQ=255)                                                                          caacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCCCTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCCTTCTGGCGCAACACCCTGGTGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGGAATGTGGAAATTGTCGTTATTGATGGTGTGGGGCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGcccaagcgggggcggggg                                                                                                                                                                 >  1:61094/1‑235 (MQ=255)                                                                           aacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCgg                                                                                                                                                                <  1:721477/251‑1 (MQ=255)                                                                           aacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCgg                                                                                                                                                                <  1:408794/251‑1 (MQ=255)                                                                              aGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGCTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGgc                                                                                                                                                             <  1:712039/251‑1 (MQ=255)                                                                               gCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGAACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGACAATGCGTGAGCGAGCGGGgcg                                                                                                                                                            <  1:543820/251‑1 (MQ=255)                                                                                  caGGCGGGAGCTGAAGCTGTGTTGATTTATCAACGCCCTGGTGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCCACCCCCTGATGTGCAGATCCTCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGTTGGGGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCTGGGGcgaaggggcgtctgtgggggggctt                                                                                                                                                          >  2:481339/1‑227 (MQ=255)                                                                                      cGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGGGCGAGCGGGGCGCTTaa                                                                                                                                                       >  1:641051/1‑249 (MQ=255)                                                                                      cGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAg                                                                                                                                                     <  1:466419/251‑1 (MQ=255)                                                                                       gggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGagcgagcg                                                                                                                                                                 <  1:410168/238‑1 (MQ=255)                                                                                       gggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGGGAGGGAGggggatcggaagagg                                                                                                                                                    >  2:410168/1‑240 (MQ=255)                                                                                              gaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTg                                                                                                                                              <  1:660455/250‑1 (MQ=255)                                                                                                 ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCATCTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGc                                                                                                                                          <  1:263992/251‑1 (MQ=255)                                                                                                   tgtgTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAAAGCGTGGGCGAGCGGGGCGCTTTAAGTCGGTTgcggcgg                                                                                                                                        >  2:489367/1‑245 (MQ=255)                                                                                                     tgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGt                                                                                                                                       >  1:541323/1‑250 (MQ=255)                                                                                                       ttGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGGTGCGGGAAt                                                                                                                                    >  1:739010/1‑251 (MQ=255)                                                                                                              aTCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCaa                                                                                                                              <  1:583144/250‑1 (MQ=255)                                                                                                              aTCAACGCACTGATGCGCCTGTGGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCaa                                                                                                                              <  2:100705/250‑1 (MQ=255)                                                                                                                cAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcg                                                                                                                            <  2:339386/250‑1 (MQ=255)                                                                                                                 aaCGCACTGATGAGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGGGCGGGGCGCTTAAAGTCGGGTGATGCGGGAATCGTCAAcggc                                                                                                                          >  1:540755/1‑251 (MQ=255)                                                                                                                   cGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggcgg                                                                                                                        <  1:191522/251‑1 (MQ=255)                                                                                                                    gaaCTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGGGGGTGCCGGC‑GGGGCCAATGCGTGAGGGGGCGGGGCGCTTTAAGGCGGTTGATGCGGGAAATGTCAAcggaggg                                                                                                                       >  2:572149/3‑247 (MQ=255)                                                                                                                          aTGCGCCTGTTGCGGTTTCTCCCGTGCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGGTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTGAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCt                                                                                                                 <  2:61094/251‑1 (MQ=255)                                                                                                                            gcgcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGGGGTTAAAGGCGGTTGTGGCGGGAATCGTCAACGGCGGGGTcccccc                                                                                                               >  2:682969/1‑248 (MQ=255)                                                                                                                              gcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAAACGTCAACGGCGGTGGCCCCCGcc                                                                                                             >  2:728205/1‑250 (MQ=255)                                                                                                                               ccTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGTGGGTGATGCGGTAATCGTCCACGGCGGGGTCCCTCGc                                                                                                              >  1:561477/1‑249 (MQ=255)                                                                                                                               ccTGTGGCGGTGTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCTTCGCGTTGGCAATGGCTGGTGGTTGCAGGC‑GGGGCCAATGAGTGAGCGAGCGGGGCGCTTAAAGTCGGTGGATGCGGTAATCGTCAACGGCGTTGTCCCGCGCAg                                                                                                            <  2:38787/251‑1 (MQ=255)                                                                                                                                  gtggCGGTTTCCCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCg                                                                                                          <  2:626225/247‑1 (MQ=255)                                                                                                                                      cGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGgtggtccccccagcgggga                                                                                                      >  1:23335/1‑234 (MQ=255)                                                                                                                                       ggTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCCACGGCGGGGGCCCCCGCAGCGGTgga                                                                                                     >  1:700756/1‑248 (MQ=255)                                                                                                                                             tcCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGCCGGTTTACAGCATTACCGTCTGCCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCCATGCGTGAGCGTGCGGGGGGCTTAGAGACGGTGGATGGGGTAATCTGAAACGGCGGGGTcccgccccgggggggaaaacccc                                                                                              >  2:290849/1‑231 (MQ=255)                                                                                                                                               ccGTTCGTTCTGAAGCGTTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGGCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCTTGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAGCGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATg                                                                                            <  1:336001/251‑1 (MQ=255)                                                                                                                                                  ttcgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCCCGCAGCGGGGAAATCCCCATGCAt                                                                                         >  1:104479/1‑251 (MQ=255)                                                                                                                                                     gttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctg                                                                                      <  1:514658/251‑1 (MQ=255)                                                                                                                                                     gttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCCCAGCGGTGAAATCCCCATGCCTctg                                                                                      >  1:16361/1‑251 (MQ=255)                                                                                                                                                          gATGCGGTAAAAGCCATTCTTGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAAGGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGcc                                                                                 <  2:452911/251‑1 (MQ=255)                                                                                                                                                          gATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTGTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctgctg                                                                                   <  1:692318/249‑1 (MQ=255)                                                                                                                                                               ggTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTACCCGCCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGCAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGGGCGAGCGGGGCCCTTACAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGGCCCTCGCAGCGGGGAAACCCCCATGCCTCTGCCGCCggggc                                                                            >  1:492108/1‑249 (MQ=255)                                                                                                                                                                    aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTTCTGCCCGGGCAGGCg                                                                       >  1:516757/1‑251 (MQ=255)                                                                                                                                                                    aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCg                                                                       <  2:390701/251‑1 (MQ=255)                                                                                                                                                                    aaGCCAGTCTGGCGGAACACCCTGAGGTGCAGATCATCGTAACCGGCGACGGTTTACAGCATTACCGTCTGGCGCGGAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCTGCGGTGAAATCCCCATGCACCTGCTGCCGGGTCAGGCg                                                                       <  2:118077/251‑1 (MQ=255)                                                                                                                                                                     aGCCAATCTGGCGCAACAACCTGATGTGCAGATCATCGAAAAAGACGACGGATTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTCATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCAGA‑GGGGCCAGTGCCTGAGCGAGGGGGGCGCTTAAAATTGGTTGCTGGGGCAATCGTCAACGGCGGGATCCCACCCACCGGTGAAAAGCCCATGcgcctgtcgtcaggtcactggg                                                                      >  2:678537/1‑230 (MQ=255)                                                                                                                                                                       ccATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCCGCGGGGAAATCCCCCTGCCTCTTCTGCCCGGTCAGGCgggg                                                                    >  1:724753/1‑249 (MQ=255)                                                                                                                                                                           tCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAAt                                                                 <  2:165917/250‑1 (MQ=255)                                                                                                                                                                               gcgcAACGCCCTGGTGTGCAGATCATCGTAACCGACGTTGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTTTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGTGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATTGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACg                                                            <  2:619174/251‑1 (MQ=255)                                                                                                                                                                               gcgcAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGTGTCCCTCGCAGCGGGGGAATCCCCATGCCCCTGCTGGCGGGGCAAGCGGGGAATTTAcc                                                            >  1:555445/1‑250 (MQ=255)                                                                                                                                                                                  cAACACCCTGATGTTCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTAc                                                         <  1:441071/251‑1 (MQ=255)                                                                                                                                                                                  cAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTAc                                                         <  1:27214/251‑1 (MQ=255)                                                                                                                                                                                      acCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGt                                                     <  1:709758/251‑1 (MQ=255)                                                                                                                                                                                        ccggATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAc                                                   <  1:544292/248‑1 (MQ=255)                                                                                                                                                                                             tgtgCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCACTTCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCCACGGCGGTGTCCCTCGCCGCGGTGGAATCCCCCTGCATCTGCTGCCGGGGCCGGCGGGGGATTTACCTACCGGGACGCGtt                                              >  1:307224/1‑251 (MQ=255)                                                                                                                                                                                                                       gtgtACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTa                    <  1:368548/248‑1 (MQ=255)                                                                                                                                                                                                                          tACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGTGAAATCCCCCTGCATCTGCTGCCGGGTCAGGCGGGGGAATTACGGACCCGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTaagg                 >  1:20373/1‑248 (MQ=255)                                                                                                                                                                                                                            cAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTGATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGc               <  1:279448/251‑1 (MQ=255)                                                                                                                                                                                                                                  tACCGTCTGGCGCGTAATGGGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATgg          <  2:738384/250‑1 (MQ=255)                                                                                                                                                                                                                                         tGGCGCGTAATGTGGACATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGGGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGGGAAATCCCCCTGCCTCTGCCGCCGGGGCCGGCGGGGGATTTTCGTACCGGGACGGGGTGTTGGGTTTGTCAGCTTGAACATGTAggggggaagggggggggg  >  2:609427/1‑234 (MQ=255)                                                                                                                                                                                                                                         gggCGCGTAATGGGGAAATTGTCTTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAAGGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAtt  <  2:641051/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    CGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGC‑GGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATT  >  NC_000913/968597‑969079

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