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%	193.4 / NA	69	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 (37/32);  total (37/32)

CTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGC  >  NC_000913/552923‑553397                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ctTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACCCGCGACGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCCCAACTATCGGCTTTAGAGGGAACAGAGAATAGCACGAAAAGAGAGCa                                                                                                                                                                                                                                         >  2:863295/1‑249 (MQ=255)    cGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGaa                                                                                                                                                                                                                                   <  1:824349/251‑1 (MQ=255)     ggCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGAGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCAGCGGCGAACTGCATGGTAGCCAAGTCGCTTTTGGACCCCATCACAATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAACACTCCTGCGGTGCACAACTCTCGGCTTTAGAGGGCCCAGAGAATAGCACGAAAAGAGAGAA‑AGGaa                                                                                                                                                                                                                                   >  1:492307/1‑250 (MQ=255)      gcgcTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGaaaa                                                                                                                                                                                                                                 <  1:715331/251‑1 (MQ=255)       cgcTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGCGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCCCGAAAAGAGAGCA‑AGGaaaa                                                                                                                                                                                                                                 >  2:83166/1‑250 (MQ=255)           ttCGGCGAAGCTGAACAGTTTATCGGGGGTGCAGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACgg                                                                                                                                                                                                                              <  1:284585/249‑1 (MQ=255)            tCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGtt                                                                                                                                                                                                                            <  1:628544/250‑1 (MQ=255)                  aaGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACACCTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAACACGGTTGCGGGG‑c                                                                                                                                                                                                                    >  2:236959/1‑251 (MQ=255)                       gAACATTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑Ctgt                                                                                                                                                                                                                 <  1:26873/249‑1 (MQ=255)                       gAACAGTGTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑Ctgt‑g                                                                                                                                                                                                               <  2:825710/250‑1 (MQ=255)                           aGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCCCAACTCTCGGCTTTAGAGGGAACAGAGAATAACACGAAAAGAGA‑CACAGGAAAACGGTTGCGAGG‑CTGT‑Gaactc                                                                                                                                                                                                          >  1:858120/1‑248 (MQ=255)                             tttATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGCAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATc                                                                                                                                                                                                          >  1:411121/1‑249 (MQ=255)                              ttATCGGGGGTGCGGTGAGCAGAAACCACCTCAACGTGGGGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACGGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCTGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGATGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAg                                                                                                                                                                                                        <  2:182987/250‑1 (MQ=255)                                  cGGGGGTGCGTGGAGCAGTAACCACTTCAACGGGGTGCGGTACATTCAGGATTTCGAAGATGTCGGCGGCGAACGGCATGGTAGCACAGTCGCGTTTGGACCCCAGCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGTGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCa                                                                                                                                                                                                      <  2:437658/248‑1 (MQ=255)                                   gggggAGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCCCAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑t                                                                                                                                                                                                 >  2:692272/1‑251 (MQ=255)                                      ggTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACAGCATGGTAGCCCAGTCGCTTTTGCACCCCATCACGAAG‑GAGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGAACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAACCGGTTGCGAGA‑CTGT‑GAAATCAGCAAAG‑TTGc                                                                                                                                                                                              >  1:109555/1‑251 (MQ=255)                                             tGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGGGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGttt‑t                                                                                                                                                                                      >  2:201399/1‑251 (MQ=255)                                                    aaaCCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAAc                                                                                                                                                                               <  1:744492/251‑1 (MQ=255)                                                     aaCCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACg                                                                                                                                                                              <  1:564511/251‑1 (MQ=255)                                                      aCCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACgg                                                                                                                                                                             <  2:220709/251‑1 (MQ=255)                                                        caattCAACGTGGGGTGGGACCTTCAGGATTTCGAAGATGTCGTCGGCGAGCGTCAAGTGAGCCCATTCGCTCTGGCACCCCATCACGATG‑GCGACACGCGTCTGATTATTTCGGGAAGACATGCGTCTTAAATCTCCTGTGGTGCTCAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGTTTGTTTGG‑TTGT‑GAAATCGGCAAAGGTTG‑GGGGTT‑TTTAAAACGGaa                                                                                                                                                                           <  2:784382/248‑1 (MQ=255)                                                        cACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGaa                                                                                                                                                                           <  1:83166/251‑1 (MQ=255)                                                            tCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGaaaa                                                                                                                                                                         <  1:766661/249‑1 (MQ=255)                                                            tCAACGGGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCACAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGAGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGAA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGaaaa                                                                                                                                                                         <  1:244591/249‑1 (MQ=255)                                                             cAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATg                                                                                                                                                                       <  1:595380/250‑1 (MQ=255)                                                              aaCGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCCCGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGA‑GAAATCAGAAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGa                                                                                                                                                                      >  1:371867/1‑250 (MQ=255)                                                                 gtggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑CGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATc                                                                                                                                                                   >  1:803359/1‑250 (MQ=255)                                                                   ggtgCGGGACATTCAGGATTTCGACGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGc                                                                                                                                                                >  1:292713/1‑251 (MQ=255)                                                                    gtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCt                                                                                                                                                               <  1:735672/251‑1 (MQ=255)                                                                      gCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTc                                                                                                                                                              <  1:801717/250‑1 (MQ=255)                                                                        gggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAAc                                                                                                                                                           >  1:115710/1‑251 (MQ=255)                                                                                  gggtttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGACGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCCGCAAAG‑GTGAGGGGTT‑TTTCAAACGGAAAATGAATCAGCCCaaagtactccgcc                                                                                                                                                 >  2:12575/4‑240 (MQ=255)                                                                                   gATTTCGAAGATTTCGGCGGCGAACTGCATGGACGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGCGGTGAACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGCGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATccgccg                                                                                                                                                >  1:266641/1‑251 (MQ=255)                                                                                     tttcgcagccttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGCTG‑GCGACACGCGCCGGATTATTGCGGGAAGCCATGCGTCTTTAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGATTAGCACGAAAAGCGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCGGCTCAACGTCACCCGCCGt                                                                                                                                               >  2:147295/11‑250 (MQ=255)                                                                                            gATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAAAGGAAAATGAATAAGCTCAACGTCATCCGACGTGACTTTc                                                                                                                                        >  1:801223/1‑250 (MQ=255)                                                                                             aTTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCAc                                                                                                                                      >  1:606856/1‑251 (MQ=255)                                                                                             aTTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGAGG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAT‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCAc                                                                                                                                      <  1:29077/251‑1 (MQ=255)                                                                                                 cggcggCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCAtt                                                                                                                                  >  1:57219/1‑251 (MQ=255)                                                                                                    cggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGa                                                                                                                                >  1:471213/1‑250 (MQ=255)                                                                                                     ggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCACATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAAc                                                                                                                              >  1:308662/1‑251 (MQ=255)                                                                                                       cgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAAAGTCATCCGCCGTGACTTTCACCATTGAACCt                                                                                                                            >  1:71164/1‑251 (MQ=255)                                                                                                           cggcagggTAGCCCAGTCTCTTTGGGACCCCATCACGATG‑GCGACACGCCCCGGATTATTGCGGGAAGACATGCTTCTTAAACCTCCTGTGGTGCACACCTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGTTTGCGTGG‑CTGT‑GAAATCAGCAAAT‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCg                                                                                                                        <  1:98746/245‑1 (MQ=255)                                                                                                           cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGACAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCCCCATTGCACCTTCCg                                                                                                                        >  2:632098/1‑251 (MQ=255)                                                                                                            tGCATGGTAGCCCAGGCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGt                                                                                                                       >  1:425095/1‑251 (MQ=255)                                                                                                             tcATGGTAACCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCTGGAAGACATGCGTCTTAAACCTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGAACGAAAAGAGAGCA‑AGGAAAACGGGTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAGACGGAAAATGAATCAGCACAACGTCATCCGCCGGGACTTTCACCATTGAACCTTACGt                                                                                                                       >  2:70029/2‑250 (MQ=255)                                                                                                                   tAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGcc                                                                                                                  >  1:321455/1‑249 (MQ=255)                                                                                                                        cAGTCGCTTTTGGACCCCATCACGATG‑GCGACACCCGCAGGATTCTTGCGGGAAGACATGCGTCTTAAAACTCCTGCGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATAAATAAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGCATGCCAGGCAcc                                                                                                           >  1:255624/1‑251 (MQ=255)                                                                                                                         aGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGCAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGGATGCCAGGCAccc                                                                                                          >  2:516870/1‑251 (MQ=255)                                                                                                                            cGCTTTTGGCCCCCATCACG‑TGTGCGACGCGCGCGGGATTATTGCGAGAAGACGTGAGTCTAAAAACTCCGGTGGTGAAAACCTCGCGGCTTTAGGGGGCCCAGAAAATAGCACAAAAAGAGAGAA‑AGAAAAAGGGTGGCGGGT‑CTGT‑GAAATCAGAAAAG‑TTGCGGGTTT‑TTTAAACGGGAAAAGGAATCAGCAAAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGt                                                                                                       <  2:336823/251‑1 (MQ=255)                                                                                                                                    gACCCCATCCCGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACCCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAcgcg                                                                                               <  1:201399/251‑1 (MQ=255)                                                                                                                                               ggaggccgagatgtaccGGATTATGGCGGGAAAGCTAGCGCCTTAAAACCCCTTTGGTGTAGAACGCTCGGCTTTGGATGGCACGGAAAATAACACGAAAGGAGAGCA‑AGGAAACCGGTTGC‑TGGTCTGT‑AAAAGCAACAA‑G‑GTGGGGGTTTGTGTAAACCGAAAAATGAACCACCCCTCCTTCATCCGCCGTGACTTTCACCAGTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGt                                                                                     <  2:408313/236‑1 (MQ=255)                                                                                                                                                      gACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAAGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGAt                                                                              <  1:593324/251‑1 (MQ=255)                                                                                                                                                             gcCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGAt                                                                       >  1:400586/1‑251 (MQ=255)                                                                                                                                                                  attattGCGGGAAGACAGGCGTCTTAAAACCCCTGGGGGGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGGGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGtt                                                                  <  1:236959/251‑1 (MQ=255)                                                                                                                                                                         ggggAAGACAGGCGTCTTAAAACTCCTGTGGGGCACAACCCTCGGCTTTAGAGGGCCCAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CGTT‑GAAATCAGCAAAT‑TGGCGGGTTT‑TTTAAAACGGAAAAGGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCAc                                                           <  1:304551/250‑1 (MQ=255)                                                                                                                                                                                 cATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGc                                                   >  2:865215/1‑251 (MQ=255)                                                                                                                                                                                    gctgcTTAAACCCCCTGTGGTGCACAACTCCCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAAAGAGCA‑AGGAAAACGTTTGCGTGG‑‑TGTTAAAATCAGCAAAT‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGAt                                                <  1:293037/247‑1 (MQ=255)                                                                                                                                                                                        cTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGt                                            >  1:85044/1‑251 (MQ=255)                                                                                                                                                                                        cTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑GTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCGGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAAGTCTTGCACCGCCGGGCGAtttgt                                            >  1:198063/1‑247 (MQ=255)                                                                                                                                                                                            aaccccgtgtgttgaaaaatccccGGCTTGAGAGGGCACAGAGAATAGCACGAAAAGAGAGAA‑AGGAAGTCGGTTGCGTGG‑CTGG‑GAAATCAGCAAAG‑TTGCGGTTTT‑TTTAAACCGGAAAATGAACCGGCTAACCGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGc                                        <  2:57219/228‑1 (MQ=255)                                                                                                                                                                                               cTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCg                                     >  2:35189/1‑251 (MQ=255)                                                                                                                                                                                                   tgtgGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCCCCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGAGATAAGCTCATGCACCGCCGGGCGATGGGTATGCCCGTggg                                 >  2:277250/1‑250 (MQ=255)                                                                                                                                                                                                     tggtgCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATg                                                                                                                                                                       <  2:54419/115‑1 (MQ=255)                                                                                                                                                                                                     tggtgCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATg                                                                                                                                                                       >  1:54419/1‑115 (MQ=255)                                                                                                                                                                                                                  ccGGCTTTAGAGGGCACAGAGAATAGAACAAAAAGAGAGCA‑AGGAAAACGGTTGCTTGG‑CGGT‑GAAACCAGCAAAG‑TTGCGGTTTT‑TTAAAAACGGAAAATTAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCt                  <  1:692272/250‑1 (MQ=255)                                                                                                                                                                                                                    ggCTTTAGAGGGCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATGGCCCCTg                 >  1:646480/1‑250 (MQ=255)                                                                                                                                                                                                                     gCTTTAGAGGGCACAGAGAATAGCACGAAAATAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGAt               >  1:54317/1‑251 (MQ=255)                                                                                                                                                                                                                               gCACAGAGAATAGCACGAAAAGAGAGCA‑AGGAAACCGGTTGCGGGG‑CTGT‑GAAATCAGCAAAG‑TTGGGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCAt     <  1:865215/251‑1 (MQ=255)                                                                                                                                                                                                                                  aagagACTAGCACGAAAAGAGAGCA‑AGGAAACCGGTTGCGTGG‑CTGT‑GAAAACAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATAAATCAGCTCAACTTTATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCCCCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCTCCTGATGTTTTTCCATCGc  <  1:393862/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATG‑GCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCA‑AGGAAAACGGTTGCGTGG‑CTGT‑GAAATCAGCAAAG‑TTGCGGGTTT‑TTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAAAGCAGGTTGCTGATTGGCGATAAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGC  >  NC_000913/552923‑553397

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