Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,211	G→A	D186D (GAC→GAT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,211	0	G	A	100.0%	247.0 / NA	71	D186D (GAC→GAT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/34);  total (37/34)

GCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCG  >  NC_000913/3814964‑3815454                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    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                                                                                                                                                                           <  1:8323/212‑1 (MQ=255)                                                                                                                cGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTCCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgcc                                                                                                                                                                           >  2:8323/1‑212 (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)                                                                                                                                                                                                                                               gggCAATCAGATTTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCa     <  2:438700/249‑1 (MQ=255)                                                                                                                                                                                                                                               aggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTTCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCa     <  1:443858/249‑1 (MQ=255)                                                                                                                                                                                                                                               aggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTc      <  2:324840/248‑1 (MQ=255)                                                                                                                                                                                                                                                  cAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCAccc  >  2:154842/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATT‑TTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCG  >  NC_000913/3814964‑3815454

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: ‑