Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,197	G→A	P323P (CCG→CCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,197	0	G	A	100.0%	162.4 / NA	52	P323P (CCG→CCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/33);  total (19/33)

ACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAA  >  NC_000913/4618950‑4619442                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aCATGAATCAGGTACTGGCCTCCAGTGCCGGTAGCGCGGTAGAAGTTCGTGAAGCGGTGCAGTTACTGACGGGTGAGGTTCGTAGCCCGCGTCTTTTTGATTACACGATGGCGCTGTGCGTGGAGATGCGTATCTCCGGCAAACTGGCTAAATATGACGCCGAAGCGCGCGCGAAATTGCAAGCATTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGTCACAAAAAGGCccaacc                                                                                                                                                                                                                                                     <  1:54896/251‑1 (MQ=255)    ggaaccAGGTACTGGCCTCCAGTGCAGTTAACGGGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGAATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTG‑AGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACt                                                                                                                                                                                                                                                 <  1:29872/246‑1 (MQ=255)      aaTCAGGTACGGGTCTCCATTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGCCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACTCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGTCAGAAGTCTTTGGTCGTATTGTAGCGGCACAAAAAGGACCAACCGACtt                                                                                                                                                                                                                                                <  2:163296/251‑1 (MQ=255)       aTCAGGTACTGGCCTCCATTGCAGGTAACGCGGTTGAAGTGCGTGAAGCGGTGCAGTGCCTGACGGGGTAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTc                                                                                                                                                                                                                                               <  2:167958/251‑1 (MQ=255)              aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgaga                                                                                                                                                                                                                                        <  2:202468/251‑1 (MQ=255)              aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgag                                                                                                                                                                                                                                         <  1:40009/250‑1 (MQ=255)               cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGAACAAAAAGGCACAACCGACTTCGTTgaga                                                                                                                                                                                                                                        >  1:52639/1‑250 (MQ=255)                         gTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGa                                                                                                                                                                                                                               >  1:92664/1‑249 (MQ=255)                          ggCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                                                                                                                                                                                                            <  2:77050/250‑1 (MQ=255)                             aGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTGCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                                                                                                                                                                                                         <  1:104550/251‑1 (MQ=255)                                 agcgcgGTTGAAGTGCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGAGGGCGCTGTGCGTGGAGATGCGTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACTGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTATCTTCTCAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                                                                                                                                                                                                     <  1:250517/249‑1 (MQ=255)                                 aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTGGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTGGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                                                                                                                                                                                                     <  2:18693/251‑1 (MQ=255)                                       gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAACGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGc                                                                                                                                                                                                               >  1:9/1‑251 (MQ=255)                                          gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGa                                                                                                                                                                                                             <  1:33209/250‑1 (MQ=255)                                               tCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                                                                                                                                                        <  1:251714/250‑1 (MQ=255)                                                 gTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACg                                                                                                                                                                                                     >  1:68623/1‑251 (MQ=255)                                                        ggTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGc                                                                                                                                                                                                >  1:5692/1‑249 (MQ=255)                                                          tGCAGATCCTGACGGGTGAGTATCGAAACAAGAGTCTGTTTGCTGTTACGATGGAGCTGTGCGTGGAGATGCTTATCCCCGGCAAACTGACGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGAAGAGCTGGACAACGGTAAAGCGGCAGCAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGACCAACCGACTTCGTTGAGAACTCAGCGAAGCATCTGCCGACAGCGATGCGGACGAAAGCCGTc                                                                                                                                                                                            >  2:129810/1‑251 (MQ=255)                                                            caggtccggACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCTCGAAATTGCAAGCAGTGCTGGATAACGGTAAAGCGGCAGAAGGCTTTGGTCTTATGGTAGCGTCACAAAAAGGCCCAACCGACTTCGTTGAGGACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTa                                                                                                                                                                                          <  1:172523/243‑1 (MQ=255)                                                             aGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCTAAATTGCAAGCAGTGCTGGACCACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt                                                                                                                                                                                         <  2:68695/251‑1 (MQ=255)                                                              gTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTGGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATg                                                                                                                                                                                        <  1:138219/251‑1 (MQ=255)                                                                tCCTGACGGGTGCGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCt                                                                                                                                                                                      >  2:133100/1‑251 (MQ=255)                                                                 ccTGACGGGTGAGTAGCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTTTCTCCGGCAAACTGGCGAAAGATTACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCTCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTg                                                                                                                                                                                     <  1:141649/251‑1 (MQ=255)                                                                    gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGa                                                                                                                                                                                    <  2:261985/249‑1 (MQ=255)                                                                       gggTGAGTATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGTGGAGATGCTTATCTCAGGCAAACTGGCGAAAGATGACGCCGCAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCg                                                                                                                                                                               >  2:182724/1‑251 (MQ=255)                                                                                   tAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTc                                                                                                                                                                    >  1:237023/1‑250 (MQ=255)                                                                                       ccGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGa                                                                                                                                                               <  1:44294/251‑1 (MQ=255)                                                                                            gcTGTTTGATGTCCCGATGGCGCGGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCTGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGTTTTTGTCAGTGaaaggg                                                                                                                                                          <  2:68810/250‑4 (MQ=255)                                                                                                      gTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgcgc                                                                                                                                                >  2:1020/1‑251 (MQ=255)                                                                                                             tGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGAGCGCTgg‑gg                                                                                                                                        >  2:76085/1‑251 (MQ=255)                                                                                                              ggCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTgg‑gg                                                                                                                                        >  1:83814/1‑250 (MQ=255)                                                                                                               gcgcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGAt                                                                                                                                      >  1:146923/1‑251 (MQ=255)                                                                                                                  cTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGACTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGc                                                                                                                                   >  1:249418/1‑251 (MQ=255)                                                                                                                    gtgcgtgaagcaggttaTCTCCGGCAAACCGGCGAAAGATGACGCCGAAGAGCGCGCGAAATTGCAAGCAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGAAGTCTATGCTGATACAGAAGGTTTTGTCAGTGAAAAGGATACCCTCGCGCTGG‑GGAGGGcgg                                                                                                                                 >  2:19650/17‑249 (MQ=255)                                                                                                                        gTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGAGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGt                                                                                                                             >  2:16774/1‑251 (MQ=255)                                                                                                                                        ccGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGggcggc                                                                                                                <  1:41727/248‑1 (MQ=255)                                                                                                                                         cGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGcggcggtc                                                                                                            <  1:152338/251‑3 (MQ=255)                                                                                                                                            cAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCAAGCAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACGTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGGCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGcggcggtcgtc                                                                                                         <  2:9/251‑6 (MQ=255)                                                                                                                                             aaaCTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGcggcggtcgtcg                                                                                                        <  2:68623/251‑7 (MQ=255)                                                                                                                                              aaCTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATTCCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGcggcggtcgtcgt                                                                                                       <  1:141156/251‑8 (MQ=255)                                                                                                                                                    cGAAAGATGACGCCGACGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGACATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGAAGCGGTCGTCGTCAGGc                                                                                                  >  1:51015/1‑250 (MQ=255)                                                                                                                                                                ccGAAGTGCGCGACAAATTGCAAGCAGTGTTGGGCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGCAGCGGCTAAAAAAGGCCCATCCGACCTCGCGGAGAACTGCGCGAAGTATCTGCCGCCCGCGATGCTGACGAAAGACGTCGACGCTGATACCGAAGGTTTTGGCAGTGAAATGGATACCCGCGC‑CTGGTGGATGGCGGTGGTTGCAAGGGGCGGCGGTCGTCGACAGGCacatgacatgagt                                                                                     <  2:106300/251‑13 (MQ=255)                                                                                                                                                                     gcgcgcgcGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCGGATACCGAAGGTTTTGTCAGGGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTa                                                                                <  2:17319/251‑1 (MQ=255)                                                                                                                                                                          cgcgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGAGGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCg                                                                           >  1:152487/1‑251 (MQ=255)                                                                                                                                                                             gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                        <  2:92664/251‑1 (MQ=255)                                                                                                                                                                               aaTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGc                                                                      <  1:93857/251‑1 (MQ=255)                                                                                                                                                                                aTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCAACCGACGTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTGTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGAGGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCt                                                                     <  2:13513/251‑1 (MQ=255)                                                                                                                                                                                                    aCGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtttg                                                 <  2:208980/251‑4 (MQ=255)                                                                                                                                                                                                    aCGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGtttg                                                 <  1:230908/251‑4 (MQ=255)                                                                                                                                                                                                                   aaGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGGTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAGACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGAc                                  <  2:237023/251‑1 (MQ=255)                                                                                                                                                                                                                         atGGTCGTATGGTAGCGGCACAAAAAGGCCCAGCCGTCGTCGTTGAGAACTACGCGAAGGATCGGCCGACAGCGATGCTGACTAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAGCTGACACCAGTGATGACAGCGTCGGCTTGTCTGATAAGGCGCGTTTGGGCGACCAGGTAGACGGTc                              <  2:249418/248‑1 (MQ=255)                                                                                                                                                                                                                          tGGTCGTAAGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGCACTACGCGAAGTATCTGCCGCCAGCGATGCTGACGAACGCAGTCTATGCTGCAACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGAAATAGGCGGCGGTCGTCGTCCGACATCTGACACCATTGATTACAGCGTCGGCTTTACTGCTATGGCCCGTTTGGGCGACCAGATAGACggac                              >  1:172634/1‑246 (MQ=255)                                                                                                                                                                                                                                                   ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGaa  <  1:133100/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGG‑GGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAA  >  NC_000913/4618950‑4619442

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