Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	370,321	G→A	G15G (GGG→GGA)	mhpB →	3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	370,321	0	G	A	100.0%	156.0 / NA	46	G15G (GGG→GGA)	mhpB	3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/22);  total (24/22)

GTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGC  >  NC_000913/370075‑370560                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                gTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATg                                                                                                                                                                                                                                             >  2:136990/1‑251 (MQ=255)       cgcgTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGAcc                                                                                                                                                                                                                                       >  1:268760/1‑251 (MQ=255)             ggCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGgcgc                                                                                                                                                                                                                                 >  1:89476/1‑251 (MQ=255)               cGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCa                                                                                                                                                                                                                                >  1:75340/1‑250 (MQ=255)                           gccccgggggaGCTGGTGCGCGCAACAAAAGGCTTCGCTGGTGGTGATGCGCCCGGATCGTTTTGTTGCCGCCCCCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGAt                                                                                                                                                                                                                   <  2:292541/241‑1 (MQ=255)                                tGCGTAGCTGGTTCGCGCAACACAATGCTTCGCAGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGt                                                                                                                                                                                                              >  1:151189/1‑251 (MQ=255)                                  cGTAGCTGGTTCGCGCAACACAATGCTTCGCGGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCCCCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCa                                                                                                                                                                                                            <  1:271506/251‑1 (MQ=255)                                       cgggTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGc                                                                                                                                                                                                       <  1:319002/249‑1 (MQ=255)                                       cTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGACCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGCGCTCGATGAGGTCAATGGc                                                                                                                                                                                                       >  1:375305/1‑251 (MQ=255)                                         ggTTCGCGAAACAAAATGCTTCGCTGGGGGTGAGGCGCCCGGATCGCTTTGTGGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGt                                                                                                                                                                                                     <  1:122937/251‑1 (MQ=255)                                           ttCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGCGGTGCTCgcttatgtaaattgcgtg                                                                                                                                                                                                    >  2:197744/1‑233 (MQ=255)                                                             tCGCTGGGGGTGATGCGCCCGGATCTCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAg                                                                                                                                                                                  <  2:238247/250‑1 (MQ=255)                                                              cGCTGGGGGGGATGCGCCGGGATCGCTTTGTTGCCGCCCCCGCCATTCCGCAAACCCGGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGc                                                                                                                                                                                 <  2:188558/250‑1 (MQ=255)                                                                   gtggtgATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTg                                                                                                                                                                           >  2:20470/1‑251 (MQ=255)                                                                           gcgcCCGGATCGCTTGGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtc                                                                                                                                                                   <  1:63418/251‑1 (MQ=255)                                                                            cgcCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtct                                                                                                                                                                  <  1:9465/251‑1 (MQ=255)                                                                             cccgggatccttttGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATtctc                                                                                                                                                                 <  2:325601/240‑1 (MQ=255)                                                                               cgggATCGTTTTGTGTCCTCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTccc                                                                                                                                                               <  2:777/249‑1 (MQ=255)                                                                                        tttGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtgg                                                                                                                                                      >  1:21218/1‑251 (MQ=255)                                                                                         ttgttgCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACtggt                                                                                                                                                     >  1:146985/1‑251 (MQ=255)                                                                                               ccgccaccgccaTTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCt                                                                                                                                               >  2:116860/1‑251 (MQ=255)                                                                                                             ccGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGcc                                                                                                                                      >  1:9882/1‑246 (MQ=255)                                                                                                               gCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCAc                                                                                                                               <  1:286424/251‑1 (MQ=255)                                                                                                                aaaaCCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGGATTGCGGCATTCTCCCCTGCACTGGTGGTGCTGTTTGCGCCAGATCACt                                                                                                                              >  2:364065/2‑251 (MQ=255)                                                                                                                              gACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTc                                                                                                                <  1:350967/251‑1 (MQ=255)                                                                                                                                cccTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCt                                                                                                               >  1:3118/1‑250 (MQ=255)                                                                                                                                        aaaCTGGCGTCGGTGATGACTCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTGTGCGCCAGAGCCCTACAACGGCGTTTTCTATGACGTg                                                                                                       <  2:58325/250‑1 (MQ=255)                                                                                                                                        aaaCTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTg                                                                                                       <  1:85619/250‑1 (MQ=255)                                                                                                                                                    gTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCt                                                                                          <  1:178586/251‑1 (MQ=255)                                                                                                                                                    gTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCt                                                                                          <  1:240798/251‑1 (MQ=255)                                                                                                                                                              tGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGt                                                                                >  1:349834/1‑251 (MQ=255)                                                                                                                                                               gACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGtt                                                                               <  1:87529/251‑1 (MQ=255)                                                                                                                                                                          aTGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTTGTGGTGCTGTTTGCGCCAGATCTCTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGc                                                                    <  1:216290/251‑1 (MQ=255)                                                                                                                                                                                        gTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGAttt                                                        <  1:330274/249‑1 (MQ=255)                                                                                                                                                                                        gTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCg                                                      >  1:12948/1‑251 (MQ=255)                                                                                                                                                                                                   aGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACGCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTTCGGc                                            >  2:113102/1‑250 (MQ=255)                                                                                                                                                                                                    gCCTGATATGCACGCTTATCTTCACTGTCTTTCCAACTCGCCGCTGGTGGGATATGTAGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGc                                            >  1:347262/1‑249 (MQ=255)                                                                                                                                                                                                               aCGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTGCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGctc                               <  2:208446/251‑3 (MQ=255)                                                                                                                                                                                                                 gCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGt                             <  1:358723/251‑1 (MQ=255)                                                                                                                                                                                                                 gCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGt                             >  1:359795/1‑251 (MQ=255)                                                                                                                                                                                                                       cTTCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCtgt                       >  2:122700/1‑251 (MQ=255)                                                                                                                                                                                                                        ttCACTGTCTTTCCCACTCGCCGCTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCtgtg                      >  2:41224/1‑251 (MQ=255)                                                                                                                                                                                                                                          cgccgcTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGa                                                                            <  1:296291/179‑1 (MQ=255)                                                                                                                                                                                                                                          cgccgcTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGa                                                                            >  2:296291/1‑179 (MQ=255)                                                                                                                                                                                                                                          cgccgcTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCtgt    >  2:127242/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccgcTGGTGGGATATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTgc  >  2:157501/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGACCCGGCGCAAGAGGTGCTCGATGAGGTCAATGGCGTGATTGCCAGCGCCCGCGAGCGTATTGCGGCATTCTCCCCTGAACTGGTGGTGCTGTTTGCGCCAGATCACTACAACGGCTTTTTCTATGACGTGATGCCACCGTTCTGTTTAGGCGTTGGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGC  >  NC_000913/370075‑370560

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