Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,199	G→A	Y190Y (TAC→TAT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,199	0	G	A	100.0%	255.1 / NA	70	Y190Y (TAC→TAT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/32);  total (38/32)

GAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGC  >  NC_000913/3814956‑3815443                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     aaaccccggctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TGGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCaa                                                                                                                                                                                                                                               <  1:213675/243‑1 (MQ=255)  aaaCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAAt                                                                                                                                                                                                                                              >  2:101046/1‑251 (MQ=255)      ccTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATc                                                                                                                                                                                                                                             >  1:438700/1‑248 (MQ=255)         gctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCt                                                                                                                                                                                                                                       >  1:212786/1‑251 (MQ=255)           tgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGCGTTTCTTTAAACGCCAAACTCTTCGAGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATACGCAATCAGATCttt                                                                                                                                                                                                                                     >  2:383990/1‑251 (MQ=255)           tgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCttt                                                                                                                                                                                                                                     >  1:112147/1‑251 (MQ=255)            ggCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTc                                                                                                                                                                                                                                    <  1:123596/251‑1 (MQ=255)            ggCAAGCCAGATCCTGGCCTTCGGTGAAGGTATGCGGGCACGTTTTGTCCGCAGCGAATTTAAATCCCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCtt                                                                                                                                                                                                                                      <  2:456528/249‑1 (MQ=255)             gCAAGCCAGATCCTGGCCTTCTGGGAAGGTAGGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGCCTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCa                                                                                                                                                                                                                                   <  1:126747/251‑1 (MQ=255)               caGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACCGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTccgg                                                                                                                                                                                                                                 >  2:196571/2‑248 (MQ=255)                       cccctgcctttagtaaaGGTATGCTGTAACTTTTTGTCCGCAGCAAATTTAAATACCGCCAGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAAGATGACGCCGGATGCCGTTTCACCCGGCGAGTTTCTTTAAACCCCAAACTCTTCGCGATAGCCCTTAACCGCCGCCAGATGTTCCGCCAGTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTgatgat                                                                                                                                                                                                                         <  2:267734/236‑1 (MQ=255)                         gtGGCCTTCGGTGAAGGTATGCTGTCACGTGTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTATTATGACGCCGGATGACTTTTCATCCGGCGAGTGTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACTGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATag                                                                                                                                                                                                                       <  1:67131/250‑1 (MQ=255)                           ggCCTTCTGTGACGGCCTGCTGTCACGTTTTGTCCGCAGCGAATATAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTATTATGACGCCGGATGACTTTTCATCCGGCGCGTTTCTTTAAACGCCAACCTCTTCGCGATAGGCCTTAACCCCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCATATAGGCAATCAGATCCTTCAGGGTGATGATag                                                                                                                                                                                                                       >  2:31687/1‑249 (MQ=255)                            gctttttgtgAGGGTATGCTGTCACGTTTTGGCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGCCTCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATTTGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAGCGCCAAACTCTTCGGGATGGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCTGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAgaga                                                                                                                                                                                                                    <  1:270436/245‑1 (MQ=255)                             ccTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGAt                                                                                                                                                                                                                   <  1:80890/251‑1 (MQ=255)                                tCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCa                                                                                                                                                                                                                 <  2:112147/250‑1 (MQ=255)                                     gAAGGTAGGCTTTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGTTT‑TGGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCTTATGACTTTTCCTCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGAGAGAGATCCCTTTg                                                                                                                                                                                                            <  2:70688/250‑1 (MQ=255)                                               tGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATc                                                                                                                                                                                                 >  1:61416/1‑251 (MQ=255)                                               tGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGAAAGAGATCCCTTTGCAGTTGTAATc                                                                                                                                                                                                 >  2:308877/1‑251 (MQ=255)                                                  cACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCAc                                                                                                                                                                                               <  1:82360/250‑1 (MQ=255)                                                             cGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCAc                                                                                                                                                                                                                >  1:157575/1‑222 (MQ=255)                                                                  ggaagttAAATCCCGCCCGACGGATGATT‑TTGACGCCCGCGGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATc                                                                                                                                                                              <  2:52861/246‑1 (MQ=255)                                                                          aatgccGCCCGACG‑AAGATTGTTGACTCCCGCTGGCTGCTAATTGCGGCACTGTTCAAGTAATATGACGCCGGATGGCTTTTCCTCCGGCGAGTGTCTTTAACCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAAGCCCCGAAAt                                                                                                                                                                      <  2:34639/247‑1 (MQ=255)                                                                          aatcccGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACTCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAAt                                                                                                                                                                      <  1:54035/247‑1 (MQ=255)                                                                             cccGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTGGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGAGGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAAtctc                                                                                                                                                                   <  1:204050/250‑1 (MQ=255)                                                                                          tgatT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGt                                                                                                                                                      >  2:154836/1‑251 (MQ=255)                                                                                                     gCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCNGa                                                                                                                                                                         >  2:228969/1‑222 (MQ=255)                                                                                                     gCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGa                                                                                                                                                                         <  1:228969/222‑1 (MQ=255)                                                                                                             gCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGa                                                                                                                                     >  1:3981/1‑250 (MQ=255)                                                                                                                   aaTTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGTCGCGCCCCCGTTCCTGACGATCAAGCGAAATCaat                                                                                                                              <  1:495415/251‑2 (MQ=255)                                                                                                                   aaTTGCGGCAAAGTTGCAGTAATATGACGCCGGATGCCTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGCTATGGCTTAACCGCCGCCAGCTGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCACTCAGATCTTTCAGGGTGATGATCGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGACTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAACTCaat                                                                                                                              >  2:438345/1‑250 (MQ=255)                                                                                                                       gCGGCACAGTTGCAGTAATATGACGCCGGATGCCTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACg                                                                                                                           <  2:333746/250‑1 (MQ=255)                                                                                                                        cGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTCCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgcc                                                                                                                                                                >  2:8323/1‑212 (MQ=255)                                                                                                                        cGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTCCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgcc                                                                                                                                                                <  1:8323/212‑1 (MQ=255)                                                                                                                          gCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGcc                                                                                                                         >  2:39885/1‑249 (MQ=255)                                                                                                                           cacaGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGc                                                                                                                      >  2:283982/1‑251 (MQ=255)                                                                                                                                             aCGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCAttagc                                                                                                    <  1:39885/251‑4 (MQ=255)                                                                                                                                               gCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCCCTTTGCAGTTGTACTCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGAAAGCGTCGCGCCATTAGCCt                                                                                                  >  2:314314/1‑251 (MQ=255)                                                                                                                                               gCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgccg                                                                                                                                                               <  1:229913/190‑1 (MQ=255)                                                                                                                                               gCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgccg                                                                                                                                                               >  2:229913/1‑190 (MQ=255)                                                                                                                                                ccGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTg                                                                                                 >  2:110792/1‑251 (MQ=255)                                                                                                                                                   gaggccTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGaat                                                                                              <  2:93769/246‑1 (MQ=255)                                                                                                                                                            tCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCAt                                                                                     >  1:435929/1‑251 (MQ=255)                                                                                                                                                             cATCCGGCGAGTTGCTTTAAACGCCAACCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGTCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCagc                                                                                    <  1:346901/251‑3 (MQ=255)                                                                                                                                                                ccGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGat                                                                                 <  1:139189/251‑2 (MQ=255)                                                                                                                                                                 cGGCGAGTTTCTTTAACCGCCAAACTCTTCGNGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGAGCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGatt                                                                                <  2:219980/251‑3 (MQ=255)                                                                                                                                                                    cGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTcgc                                                                             <  1:283982/251‑1 (MQ=255)                                                                                                                                                                      aGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTcgcg                                                                            >  2:104180/1‑250 (MQ=255)                                                                                                                                                                         ttCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATc                                                                        >  1:160135/1‑251 (MQ=255)                                                                                                                                                                             ttAAGCGCCAAACTCTTCGCGAGAGGCCTTAACCCCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGAGAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGcc                                                                     <  2:162532/250‑1 (MQ=255)                                                                                                                                                                                aaCGCAAAACTCGTCGGGATAGGCCTTAACCGCCGCCAGAAGTTCCGCCATTTACGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATCCAGATCACTATGCAGTTGTAATAACGCTCAACTTCCTGAATCGCCGAAATCTCGCCAAGCCCCCGTTCCTGACGATAAAGCGAAATCAATACGCCAGCAAGCGCCGAGCCATTCGCCTGAATAAACTCCATCGATTCGCCAATCGCCCTTc                                                                 >  2:493590/1‑251 (MQ=255)                                                                                                                                                                                    ccAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATANTCTCCATCGATTCGCGAATCGCCGTTCcggc                                                             >  2:219368/1‑251 (MQ=255)                                                                                                                                                                                           ctTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATc                                                      >  1:204904/1‑251 (MQ=255)                                                                                                                                                                                             tCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGAt                                                       <  1:101046/248‑1 (MQ=255)                                                                                                                                                                                              cgcgATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGACATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCATGAATAATCTCCATCGATTCGCGACTCGCCGTTCCGGCGGTGATcacc                                                   >  2:398117/1‑250 (MQ=255)                                                                                                                                                                                                 gATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCAcatc                                                 >  1:98382/1‑250 (MQ=255)                                                                                                                                                                                                  aTAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACACTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCAcatcac                                               >  1:455165/1‑250 (MQ=255)                                                                                                                                                                                                   tAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGNGa                                                                           >  2:245325/1‑222 (MQ=255)                                                                                                                                                                                                   tAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGa                                                                           <  1:245325/222‑1 (MQ=255)                                                                                                                                                                                                    aGGCCTTAACCGCCGCCAGATGGTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAAACAGATCTTTCAGGGTGATGATAGCGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAACTCTCGCCGCGCCCCCGTTCCTGACGATCCAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGCATCGCCGTTCCGGCGGTGATCCCCTCATCt                                             >  2:310817/1‑251 (MQ=255)                                                                                                                                                                                                          aaaCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGc                                       <  2:98382/250‑1 (MQ=255)                                                                                                                                                                                                            cccgccgCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCAt                                     <  2:3981/250‑1 (MQ=255)                                                                                                                                                                                                               gccgccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcg                                                               <  1:408139/222‑1 (MQ=255)                                                                                                                                                                                                               gccgccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCCAATCGCCGTTCcg                                                               >  2:408139/1‑222 (MQ=255)                                                                                                                                                                                                                       aTGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCCGCATTACGCGTCCtt                          >  2:190797/1‑251 (MQ=255)                                                                                                                                                                                                                         gTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGt                        >  1:324840/1‑251 (MQ=255)                                                                                                                                                                                                                            ccGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAAc                     >  2:300534/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTTTCCGGCTTCTCTACCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCCGCATTACCCGCCCTTGTAAcgggcc                >  2:414834/1‑247 (MQ=255)                                                                                                                                                                                                                                          ttCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTaccaacc         >  1:378564/1‑249 (MQ=255)                                                                                                                                                                                                                                               ttCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACAAGCATTACGCGTCCTTGTAACGCGCTACCAACCCGATTgc  >  1:266723/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGC  >  NC_000913/3814956‑3815443

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