Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,364	C→T	E135E (GAG→GAA)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,364	0	C	T	100.0%	196.2 / NA	58	E135E (GAG→GAA)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/34);  total (24/34)

GAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTC  >  NC_000913/3815119‑3815606                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCATCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTcgcg                                                                                                                                                                                                                                                <  1:126951/251‑1 (MQ=255)        ttaactccacaacctctTCGCGATGTGTCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTACTCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCCATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGc                                                                                                                                                                                                                                          <  2:251352/237‑1 (MQ=255)          caacccccaaacccttccgaaTAGGCCTTACCCGCCGCCAGAGGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCACCTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCg                                                                                                                                                                                                                                        <  2:160803/230‑1 (MQ=255)             aCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGAGCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTcc                                                                                                                                                                                                                                    <  2:95855/251‑1 (MQ=255)               cccAAACTCTACGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCcgg                                                                                                                                                                                                                                  <  1:363546/250‑1 (MQ=255)                  aatctctTCGCGAGAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCGCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAAGCGCCGAAATCGCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑AGCGATTCGCGAATCGCCGTTCcggcgg                                                                                                                                                                                                                               <  2:97337/248‑1 (MQ=255)                     ctctTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTg                                                                                                                                                                                                                             >  1:356393/1‑250 (MQ=255)                         tCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATcac                                                                                                                                                                                                                        >  2:152678/1‑251 (MQ=255)                                  gCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTAc                                                                                                                                                                                                               <  1:22246/251‑1 (MQ=255)                                    cTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATc                                                                                                                                                                                                                          <  1:113953/238‑1 (MQ=255)                                    cTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGGGATc                                                                                                                                                                                                                          >  2:113953/1‑238 (MQ=255)                                           gccgccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACCTCATCTACCAGCAttcc                                                                                                                                                                                                      >  2:74729/1‑249 (MQ=255)                                           gccgccAGAGGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTAc                                                                                                                                                                                                      <  2:130700/251‑1 (MQ=255)                                               ccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGt                                                                                                                                                                                                  <  1:481731/251‑1 (MQ=255)                                                      ttCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGt                                                                                                                                                                                            >  1:168690/1‑250 (MQ=255)                                                         cgctaTTTCCGGGTTCTCTTCCAGATACGCAATCAGATCTTTCAGGGTGATGATAGACAACACTTTGCAGTTGTAATCACGCTCAACTTGCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAACCCAATACTCCAGCAACCGTCGCGCCATTAGCCTGACTCATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGCGGACACATCATATACCCGCATTACCCGTccacgtaa                                                                                                                                                                                          >  2:482581/5‑243 (MQ=255)                                                           ccATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAAcg‑cc                                                                                                                                                                                     >  1:158919/1‑250 (MQ=255)                                                               ttCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTa                                                                                                                                                                                  <  1:146978/250‑1 (MQ=255)                                                                    gCTTCTCTTCCAGATAGCCAATCAGCTCTTTAAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCAAGCTCAACTTCCTGAATCGACGAAATCTCGCCGCGCACCCGTTCCTGACGATCAAGCGAACTCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCCACGTTCCGGCGGTGATCCCCTCATCTAACAGCATTACGCCTCCTTGTAACG‑CGCTaccaacc                                                                                                                                                                            >  2:274004/1‑251 (MQ=255)                                                                                  aaggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAATCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTc                                                                                                                                                              <  2:146079/248‑1 (MQ=255)                                                                                   aggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTcc                                                                                                                                                             >  2:258214/3‑250 (MQ=255)                                                                                        aTCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCATCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTg                                                                                                                                                        <  1:73052/251‑1 (MQ=255)                                                                                          cagatCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCACCTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTgg                                                                                                                                                       <  2:421247/246‑1 (MQ=255)                                                                                               cTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTAACAACCAGATTGCCGCCTTCACCGTGGTCtttt                                                                                                                                                 >  1:350526/1‑251 (MQ=255)                                                                                                    aGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATAATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctc                                                                                                                                            >  1:390435/1‑248 (MQ=255)                                                                                                       gTGATGATAGAGATCACTTTGCAGTTGGAAGCACGCTCAACTTCCTGAATCGCCGAAAGCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCGTTTGCCTCttt                                                                                                                                         <  2:36264/251‑1 (MQ=255)                                                                                                            gatAGAGATCACTTTGCAGTTGTAATCCCGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGt                                                                                                                                    <  1:122296/251‑1 (MQ=255)                                                                                                            gatAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCTCCGTGGTCTTTTGCCTCTTTGCGGt                                                                                                                                    <  1:87967/251‑1 (MQ=255)                                                                                                            gatAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGt                                                                                                                                    >  1:435376/1‑251 (MQ=255)                                                                                                            gatAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGt                                                                                                                                    <  1:468255/251‑1 (MQ=255)                                                                                                                gagaTCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACCTCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCCGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTaa                                                                                                                                 >  1:283666/1‑250 (MQ=255)                                                                                                                   aTCACTTTGCAGTTGTAATCACGCTCAACTTCCTGACTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGGGGTCTTTTGCCTCTTTGCGGTTAAAGc                                                                                                                              >  1:360322/1‑250 (MQ=255)                                                                                                                        ttttcAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCAGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAACCTGAATAATCTCC‑ATCGATTAGCGAATCGCCGTTCCGGCGGCGATCACAGCCTCTACCAGCTATACGCGACCTTCTTACG‑CGCTACCAACCAGAATGCCGCCTTCACCGTGGGCTTTTCCCTCttttctgtttactcagtacg                                                                                                                        >  2:296712/5‑234 (MQ=255)                                                                                                                            cAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGc                                                                                                                      >  1:345057/1‑249 (MQ=255)                                                                                                                              gTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                                                                                          >  1:470121/1‑227 (MQ=255)                                                                                                                              gTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                                                                                          <  2:470121/227‑1 (MQ=255)                                                                                                                                  tAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTTCGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGGGGGCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGccgggcaa                                                                                                               >  2:83703/1‑243 (MQ=255)                                                                                                                                       ccGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGcaggtcaaggtcat                                                                                                         <  1:152678/250‑3 (MQ=255)                                                                                                                                             aaCTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCt                                                                                                   <  1:17763/251‑1 (MQ=255)                                                                                                                                                 tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCACGCGCCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGCTCACATCATCTACCCGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTTAAGCAGTACGGCAGGTCAAGGTCATGATCCTCCGc                                                                                               >  1:439511/1‑251 (MQ=255)                                                                                                                                                  ccTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCCGATTGCCGCCTTCACCGGGGTCTTTTGCCTCTTTGCGGGTAAAGAAGTACGGcagggtaagggcctggatgcacccc                                                                                              >  2:60538/1‑230 (MQ=255)                                                                                                                                                      aaTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAgcg                                                                                          >  1:349876/1‑249 (MQ=255)                                                                                                                                                              aaaTCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGGCAGGTCATGATGCTCCGCCCGCGCcacgggg                                                                                  >  1:228067/1‑247 (MQ=255)                                                                                                                                                              aaaTCTCGCCGCGCCCCCGGTCCTGACGAGCAAGGGAAATCAATACGCCAGCAAGGGTCGCGCGATGGGGCTGAGTAGTCGCC‑AGCGATGCGCGAATCGCCGTTCGGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAGCG‑CGCTACCAGCCAGATGGCCGCCTTCACCGGGGTCTTGTGCCTCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGgtggtggtggtcggggggggcggcgggg                                                                                  <  2:158919/251‑23 (MQ=255)                                                                                                                                                                                        cGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCttt                                                        <  1:100923/251‑1 (MQ=255)                                                                                                                                                                                          atcaaGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGt                                                      <  1:36864/247‑1 (MQ=255)                                                                                                                                                                                                cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAg                                                <  1:61057/251‑1 (MQ=255)                                                                                                                                                                                                cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTCTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCGCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAg                                                <  1:397270/251‑1 (MQ=255)                                                                                                                                                                                                   aatcaatccGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCGTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGc                                             <  2:317140/243‑1 (MQ=255)                                                                                                                                                                                                      caatacGCCAGCAAGCGTCGCGCCATTAGCCTGAATAACCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCaa                                          <  2:177370/247‑1 (MQ=255)                                                                                                                                                                                                          ccGCCAGCAAGCGTCGCGCCATGAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGTTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTTAGCAGGGCCAAAca                                       <  2:429597/249‑1 (MQ=255)                                                                                                                                                                                                                   aGCGTCGCGCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTGCCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTGTGTAAGCAGGGCCAAACAGCAGATCGaa                             <  1:396111/251‑1 (MQ=255)                                                                                                                                                                                                                        cgcgCCATTAGCCTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCaa                        <  1:149029/251‑1 (MQ=255)                                                                                                                                                                                                                         ccgccatgagCCTGAATAA‑CTCCCATCGATTCTCGAATCTCCGTGCCGGCGGTGGTCACACCATCTACCAGCATTCCGCGTCCTT‑TTACGTCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCGTTGCGGTTAAAGCAGTACGTCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATGGGGATCCCTTTGTAAGCATGGCCAAACTGCAGAGCGAACTCAAt                       <  2:475419/242‑1 (MQ=255)                                                                                                                                                                                                                                    cTGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCCCCGCCAGCGCCACGGCGGGAGGGGGCGCCAACGGGGTCCCTTTGTAAGCAGGGGCAAAACGCAGATAGAACTccattgcagaatccc            >  1:171302/1‑237 (MQ=255)                                                                                                                                                                                                                                     tGAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGaa                             <  2:356827/233‑1 (MQ=255)                                                                                                                                                                                                                                      gAATAATCTCC‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGTCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCTGTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATTCCAGAATcctcc          <  2:257230/251‑4 (MQ=255)                                                                                                                                                                                                                                               cc‑ATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCATGATGCTCCGCCAGCGCCACGGCGGGAGTGGGCGCCATCGGGATCCCCTTTTAAACAGGGCCAAACAGCAGAACGAACTCAACTCCAGAATTCACCAAAGCTTc  >  1:377798/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCC‑ATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACG‑CGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTC  >  NC_000913/3815119‑3815606

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