Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,789,286	G→A	G101D (GGT→GAT)	glaH →	glutarate dioxygenase GlaH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,789,286	0	G	A	100.0%	170.2 / NA	53	G101D (GGT→GAT)	glaH	glutarate dioxygenase GlaH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/20);  total (33/20)

GCGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGGTGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACA  >  NC_000913/2789040‑2789531                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTc                                                                                                                                                                                                                                                    <  1:564047/251‑1 (MQ=255)    gATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCgat                                                                                                                                                                                                                                                 >  2:48623/1‑251 (MQ=255)       tcaccctcacccCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTCAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCgatgat                                                                                                                                                                                                                                              >  2:289683/1‑251 (MQ=255)           cctcacccCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGGTGCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAAGCGTTTCTGCGGTTGCGGGTAGCCAAAATTCGTGACGATCTATGTGCGAACCAGCTGCAACCATGGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGa                                                                                                                                                                                                                                          <  1:584211/251‑1 (MQ=255)            ctcacccCGACGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGGTGCCGAGGGGCCCGTGCAGGCGCTGGAGAACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCACCCATTGCTGTTGAAGACGTTGCTCAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCgatgagttga                                                                                                                                                                                                                                          >  2:477456/1‑245 (MQ=255)                 cccGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGaagcaag                                                                                                                                                                                                                                    <  1:407856/251‑1 (MQ=255)                        gcgcAATCCCCGCGTCCGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATga                                                                                                                                                                                                                             <  1:113214/251‑1 (MQ=255)                                     gTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTCAGCAGGTTGCCGAGTGGAACGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGACAAAATTCTTGACGATCTCTGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGACGCAAGCGGATGAGCTGGTGAAGCTg                                                                                                                                                                                                                >  1:252315/1‑251 (MQ=255)                                        tgctgGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCg                                                                                                                                                                                                             >  1:303441/1‑251 (MQ=255)                                           tgGAACTCCCCTT‑ACCAGAAAAGACGACCAAACACTTTCTTGAGCAGGTTGCCGCGTGGCCCGTGCAGGCGCTGGAGTCCAACTCGTTTCTGCGTTTTCGGGTAGCCAAACTTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGCTGTGAAGAAAGCGGTTGAGGTGGGGAAGCTGGCgccg                                                                                                                                                                                                          >  2:546368/1‑248 (MQ=255)                                                 tCACCTTCACC‑GAACAGACGACCAAACCGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGACGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGAAAGCGGATGAGATGGGGAAGCTGGCGGCGGCGGTg                                                                                                                                                                                                    >  2:27892/1‑251 (MQ=255)                                                  cACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGGTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGt                                                                                                                                                                                                     <  1:295772/249‑1 (MQ=255)                                                         aCC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACACATCGTTTCTGCGTTTTCGGGTAGCCACAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTTGCGCATCt                                                                                                                                                                                            >  2:222225/1‑251 (MQ=255)                                                                    cgacCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGGTGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGt                                                                                                                                                                                  >  2:230416/1‑251 (MQ=255)                                                                    cgacCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGGGGCGCATCTGATTGGCCgg                                                                                                                                                                                  >  2:430657/1‑250 (MQ=255)                                                                     gacCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGGGGCGCATCTGATTGGCCGtt                                                                                                                                                                                 >  1:17533/1‑251 (MQ=255)                                                                             aGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTc                                                                                                                                                                         <  1:103291/251‑1 (MQ=255)                                                                             aGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTc                                                                                                                                                                         >  1:322990/1‑251 (MQ=255)                                                                                ttCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTcgat                                                                                                                                                                      >  1:304258/1‑251 (MQ=255)                                                                                        cAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTcgatgcgatg                                                                                                                                                                >  1:31055/1‑249 (MQ=255)                                                                                        cAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAg                                                                                                                                                              <  2:264275/251‑1 (MQ=255)                                                                                         aGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGc                                                                                                                                                             <  1:48623/251‑1 (MQ=255)                                                                                              gCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCa                                                                                                                                                       <  2:467044/251‑1 (MQ=255)                                                                                                  tctGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCATGACGATTTATGTGCAAACCATCGGCAACCATTGCGGTGGAAGACGCTGCTAAACCGCCCTGAAGGTGCGTGGTGGACCAAGGCGGTGCAAGTCGATGATGTGAATCGAGCTGATGAGATGGAGAAGCTGGCGACGGCGGTGGCGCATCTGATTGACCGGTCCAATTGCCAAGCGATGAGCGT‑TCAGTa                                                                                                                                                    <  2:523212/248‑1 (MQ=255)                                                                                                     ggCCCGTGCAGGCGCTGGAGTACAACTCGTTTCTCCGCTTTCGGGTAGCCCCAATTCTTGAAGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTCGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGAGGTGGATGTCGCTGATGTGGCGCAAGCGGATGAGATGGTGAAGCTGGCAACGGCGGTGGCGCATCTGCTTGGCCGTTCCAATTTCGTTGAGATGAGCGG‑TCAGTa                                                                                                                                                    >  1:523212/1‑247 (MQ=255)                                                                                                             cAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCACAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGttt                                                                                                                                        >  1:251757/1‑251 (MQ=255)                                                                                                              aGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGtttt                                                                                                                                       <  1:27892/251‑1 (MQ=255)                                                                                                                             aaTCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTc                                                                                                                        >  1:479233/1‑251 (MQ=255)                                                                                                                                      tGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGACCCCGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGAGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGCGGCGAATCTGATTGGCAGTTCCAATTTCGATGCGATGAGCGG‑GCAGGATTACGCGAGTTTTGTAGTGAAAAATGTCGATAACTc                                                                                                                >  1:255014/1‑250 (MQ=255)                                                                                                                                             tCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGc                                                                                                         <  1:537894/250‑1 (MQ=255)                                                                                                                                                   aGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTgc                                                                                                  <  2:479233/251‑1 (MQ=255)                                                                                                                                                    gCCAAAATTCTTGACGATCTATGTGCGAACCATCTGAAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGAGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTTGCGCATCTGATAGGCCGTTCCAATTTCGATGCGATGAGAGG‑TCAGTATTACGGGCGTTTTGGAATGTAAAATGTCGATAACTCAGaacgctaactgcc                                                                                                 >  2:439153/1‑239 (MQ=255)                                                                                                                                                     ccAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGGGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTgcgc                                                                                                >  1:66713/1‑251 (MQ=255)                                                                                                                                                       aaaaTTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGACGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGcc                                                                                               >  1:394211/1‑250 (MQ=255)                                                                                                                                                          aTTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCCGTAGTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGcc                                                                                           >  2:157696/1‑251 (MQ=255)                                                                                                                                                                        aTGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATgg                                                                             <  2:10694/251‑1 (MQ=255)                                                                                                                                                                        aTGTGCGAACCAGCTGCAACCATTGCTGTGGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATg                                                                              <  1:506800/250‑1 (MQ=255)                                                                                                                                                                        aTGTGCGAACAAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGACGGTGCGTTGTTGATCAATGAGGTGGATGTCGCTGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGACGGGGGCGAATCTGATTGGACGTTCCAATTTCGATGAGATGGGCGGGTC‑GTATTACGCGCGGTTTGTAGTGAAACATGTCGTTAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGAGATgg                                                                             >  2:219696/1‑251 (MQ=255)                                                                                                                                                                             cGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTACACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACCGCTATCTGCGCCAGCCGCCCCGCGTGATGGAGCt                                                                         >  2:410302/1‑250 (MQ=255)                                                                                                                                                                                  cAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTacaaa                                                                    >  2:456906/1‑248 (MQ=255)                                                                                                                                                                                       gCAACCATTGATGTTGCAGACGTTGCTAAACCGCGCTGACGGTGCGGTGTTGATCAATGCGGTGGATGTCGACGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGGGGCGCCTCTGATTGGCCGTTCCACTTTCGCTGCGATGCGCGG‑TCAGGTTTACGCGCGTTTTGTAGGGAAAAATGTCGATAACTCAGACAGCGATCTGCGCCAGCCGCACAGCGTGATGGCGCTCCACAacgacg                                                              >  2:513621/1‑251 (MQ=255)                                                                                                                                                                                         aaCCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGc                                                            <  1:430657/251‑1 (MQ=255)                                                                                                                                                                                           ccATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCAc                                                          >  1:474529/1‑251 (MQ=255)                                                                                                                                                                                              ttGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGGGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATCACTCAGACAGCTATCTGCGCCAGCCGCCCCGCGTGATGGCGCTACACAACGACGGCACgg                                                        >  1:103046/1‑249 (MQ=255)                                                                                                                                                                                               tGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTa                                                       <  1:8681/250‑1 (MQ=255)                                                                                                                                                                                               ggCTGTTGAAGACGTTGCTAAACCGCGCGGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATGTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTAt                                                      <  2:17533/250‑1 (MQ=255)                                                                                                                                                                                                          aCGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGAt                                           <  1:289683/251‑1 (MQ=255)                                                                                                                                                                                                                  aaaCCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGAtt                                   >  1:60101/1‑251 (MQ=255)                                                                                                                                                                                                                     ccGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGTGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGGTTTGTAGTGAAAAATGACGATAACTAAGACAGCTATCTGCGCCCGCCGCACCGCGTGATGGAGCTACACCACTACGGCACGTATGTCGAAGAGATAACCGATTACg                                >  2:516881/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGCTGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTCTCTGCGCCCGCCGCACCGCGTGATGGAGCTTCACAAAGCCGGCCCGTATGTCGAAGAGATAACAGCTTACGTGCtgatgatga                    >  1:354751/1‑250 (MQ=255)                                                                                                                                                                                                                                  gCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCCGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCtgatgatga                    >  1:380101/1‑250 (MQ=255)                                                                                                                                                                                                                                      tgttgATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATc               >  1:264550/1‑251 (MQ=255)                                                                                                                                                                                                                                             cAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGc        <  1:23109/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ggtggatgTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACa  <  1:546368/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACC‑GAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGGTGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGG‑TCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACA  >  NC_000913/2789040‑2789531

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