Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	553,161	G→A	intergenic (‑61/+57)	purE ← / ← lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	553,161	0	G	A	100.0%	178.1 / NA	55	intergenic (‑61/+57)	purE/lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/21);  total (34/21)

TGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTG  >  NC_000913/552870‑553446                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                               tGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAagaga                                                                                                                                                                                                                                                                                         >  2:541271/1‑299 (MQ=255)                             gCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaaa                                                                                                                                                                                                                                                            <  1:596954/300‑1 (MQ=255)                             gCAATAATCACATGATAACCGTTCTCTTCGGCGCTTTCGTCGAAGCGGAACAGGTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGTTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGGGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaaa                                                                                                                                                                                                                                                            <  1:596962/300‑1 (MQ=255)                                ataatCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCa                                                                                                                                                                                                                                                         >  2:525321/1‑300 (MQ=255)                                  aataaCCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TATTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATccgc                                                                                                                                                                                                                                                       >  2:570950/5‑297 (MQ=255)                                       ccTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGt                                                                                                                                                                                                                                                  >  1:568342/1‑300 (MQ=255)                                                 gTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGtttt                                                                                                                                                                                                                                        >  2:616819/1‑300 (MQ=255)                                                                            gAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAAc                                                                                                                                                                                                             >  2:646939/1‑300 (MQ=255)                                                                             aaCAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACg                                                                                                                                                                                                            <  1:646939/300‑1 (MQ=255)                                                                                 gTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCAt                                                                                                                                                                                                        >  2:427409/1‑300 (MQ=255)                                                                                   ttATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATcc                                                                                                                                                                                                      <  1:225680/300‑1 (MQ=255)                                                                                         ggggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTCAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGt                                                                                                                                                                                                 >  1:149256/1‑299 (MQ=255)                                                                                         ggggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGt                                                                                                                                                                                                 >  1:149010/1‑299 (MQ=255)                                                                                              gCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGCGACttt                                                                                                                                                                                           >  2:290457/1‑300 (MQ=255)                                                                                                ggTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCa                                                                                                                                                                                         >  1:591773/1‑300 (MQ=255)                                                                                                         aaaCCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTNGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAAc                                                                                                                                                                                <  1:126783/300‑1 (MQ=255)                                                                                                               cTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGCACCTTcc                                                                                                                                                                           >  1:507452/1‑299 (MQ=255)                                                                                                               cTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGAGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGACGTGACTTTCACCATTGAACCTTcc                                                                                                                                                                           >  2:356945/1‑299 (MQ=255)                                                                                                                         gtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATCGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTAAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGcc                                                                                                                                                                    >  2:350895/1‑296 (MQ=255)                                                                                                                         gtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGc                                                                                                                                                                >  1:331726/1‑300 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGTTAGCCCAGTCGCTTTTGTACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  2:624245/300‑1 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGTTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  1:313410/300‑1 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATTGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGTTTGCGCGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACGTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  1:570950/300‑1 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       >  1:27875/1‑300 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       >  2:6436/1‑300 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTGTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAAATCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  1:563413/300‑1 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCACTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  2:536118/300‑1 (MQ=255)                                                                                                                                  aTTCAGGATTTCGAAGATTGCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  1:414527/300‑1 (MQ=255)                                                                                                                                    tcgggATTTCGAAGATTTCGGCGGCGAACTGCGTGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGTAAGACATTCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTGGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGCCATCCGCCTCGACTTTCACCATTGAACCTTCCGTATGCCATGCACCCAGTAc                                                                                                                                                       <  2:180792/295‑1 (MQ=255)                                                                                                                                       ggATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCAACATTGAACCTTCCGTATGCCAGGCACCCAGTACCAcgc                                                                                                                                                  >  2:273896/1‑300 (MQ=255)                                                                                                                                                        gcggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATANCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAAAGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGa                                                                                                                                 >  2:466903/1‑300 (MQ=255)                                                                                                                                                             gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       >  1:94577/1‑273 (MQ=255)                                                                                                                                                             gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc                                                                                                                                                       <  2:94577/273‑1 (MQ=255)                                                                                                                                                             gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGc                                                                                                                            >  2:19266/1‑300 (MQ=255)                                                                                                                                                             gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAAGCGAAAAGCAGGTTGCTGATTGGc                                                                                                                            >  2:650443/1‑300 (MQ=255)                                                                                                                                                                    aTGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGNGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGACTCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGt                                                                                                                     >  2:453966/1‑300 (MQ=255)                                                                                                                                                                     tGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGtt                                                                                                                    >  1:23980/1‑300 (MQ=255)                                                                                                                                                                     tGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGtt                                                                                                                    >  2:514236/1‑300 (MQ=255)                                                                                                                                                                     tGGTAGCCCAATCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGtt                                                                                                                    >  1:65112/1‑300 (MQ=255)                                                                                                                                                                                                tCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTACGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGACGGGCGA‑TGGGTATGcc                                                                                        >  2:341039/1‑300 (MQ=255)                                                                                                                                                                                                tCACGAGGGCGACACGCGCCGGATTATTGCTGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGcc                                                                                        <  2:268197/300‑1 (MQ=255)                                                                                                                                                                                                        gCGACACGCGCCGGATTATTGCGGTAAGACGTGCG‑TCTTAAAACTCCTTGGGTGCACAACTCTCGGCTTTATATGGCACAGAGAATAGCACGAAAAGAGAGCAATGAAAACAGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGcc                                                                                        <  2:268128/292‑1 (MQ=255)                                                                                                                                                                                                          acccgcgtcggattcnctGCGGGAAGACATGCGCNCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGANTGGGTATGCCCGTGGa                                                                                  <  1:537274/281‑1 (MQ=255)                                                                                                                                                                                                                                              ttAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGacc                                            <  1:525321/299‑1 (MQ=255)                                                                                                                                                                                                                                              ttAAAACTCCTGTGGTACACAACTCTCGGCTTTAGAGGGCACAGAGGATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGGGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAAGCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGTGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGacc                                            <  2:85707/299‑1 (MQ=255)                                                                                                                                                                                                                                                         gtggtgCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGtttt                                  >  2:379664/1‑298 (MQ=255)                                                                                                                                                                                                                                                          tggtgCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTGGGt                               <  1:40934/300‑1 (MQ=255)                                                                                                                                                                                                                                                            gtgCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTa                             >  2:583854/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          gCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGc               >  1:23492/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          gCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCACGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTACCGTCCATGAt                  >  2:89022/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                  gggCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGTGTTTTGGTTAACGTCCATGATCGCCAGCGAtt       >  2:96595/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                    gCACAGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTa     >  2:548055/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                      acaGAGAATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGa                   <  1:518126/284‑1 (MQ=255)                                                                                                                                                                                                                                                                                          agaATAGCACGANAAGAGAGCAAGGAAAACGNTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATNATCGCCAGCGATTTAc    <  1:541271/295‑1 (MQ=255)                                                                                                                                                                                                                                                                                           gaATAGCACGAAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAAGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTActg  >  2:357591/1‑296 (MQ=255)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                               TGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCG‑TCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGA‑TGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTG  >  NC_000913/552870‑553446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑