Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,393	C→T	L389L (CTG→TTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,393	0	C	T	100.0%	134.8 / NA	44	L389L (CTG→TTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/28);  total (16/28)

TAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAAC‑ATATGAAACGTGCATTTATTATGGTGC  >  NC_000913/4619149‑4619627                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                tAAAGCGGCCGAAGGCGTTGGTCGAATGGGAACGGCCCAAAAAGGCCCCACCGACTTCGTTGAGAACGCAGCGAAGTATCAGCAGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGATGGGGATGGCAGTGGTTGCCATGGGCGGCGGTCTACTTCCGGCATCTGCCACCATTGTTTACAGCGTCGGCGTTACTGCTATGGCGCGTTTggggg                                                                                                                                                                                                                                           >  2:492891/1‑249 (MQ=255)                  ttGGCCGTCTGGTACCGGTACAAAAGTGCCCGACCGACTTCATGGAGACCTCCGCGAATCACCTGACGTCCGCCATGCTGACGAACGCAGTCTATGCTCATACCGAAGGTTTTGTCAGGGAATCGGTTACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCTTCGTCAGGCATCTGACGCCATTGATTACAGCGTCGGATTTACTGATATGGCGCGTTTGGGCGagcattgagccggtca                                                                                                                                                                                                                           <  2:62730/250‑16 (MQ=255)                   tGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAg                                                                                                                                                                                                                          <  1:274465/250‑1 (MQ=255)                               aGCGGCACAAAAAGGCCCAACCGGCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg                                                                                                                                                                                                              <  1:379803/250‑1 (MQ=255)                               aGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGTCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg                                                                                                                                                                                                              <  1:77781/250‑1 (MQ=255)                                  ggCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTa                                                                                                                                                                                                          <  1:19236/251‑1 (MQ=255)                                      caAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGCCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCATTACTGATATGGCGCGTTTGGGCGGCGAGGTAGACGGGCAGCGTCCGCTGGCGGTTaaccc                                                                                                                                                                                                      >  2:193677/1‑247 (MQ=255)                                           aGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTGCGCTGGGGGGTATctcaccgg                                                                                                                                                                                                 >  2:47189/1‑244 (MQ=255)                                            ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAcgcgc                                                                                                                                                                                                 >  1:476055/1‑249 (MQ=255)                                            ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGGTATCCAcgcgg                                                                                                                                                                                                 >  2:437241/1‑249 (MQ=255)                                                caaccgacTTCGTTGAGAACTACGCGAATTATCTGCCGACAGCGATGCGTACGAAAGCAGTCTATGCTGATACCGAATGTTTGGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTTCTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGAc                                                                                                                                                                                           <  2:406831/249‑1 (MQ=255)                                                            gtGAGAACTACGCGAAGTATCGGCCGACAGCGATGATGACGAAAGCAGTATATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGACAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCT‑AAAGACGAAAACAGCTgg                                                                                                                                                                               <  2:434412/250‑1 (MQ=255)                                                                 aaCTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAgg                                                                                                                                                                           >  1:214521/1‑250 (MQ=255)                                                                   cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCCCG‑AAAGACGAAAAAAGCTGGCAGGa                                                                                                                                                                          >  1:186712/1‑249 (MQ=255)                                                                     aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAg                                                                                                                                                                        <  1:495334/249‑1 (MQ=255)                                                                     aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAgcg                                                                                                                                                                      <  2:82640/251‑1 (MQ=255)                                                                            gtttcgtcctgcagtgttGCTGACGAAAGCAGTCTAGGCTTATGCAGAATGTTGTGGCAGTGAAATGGAGACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGGGGCGGGCGTCGTCAGGCATCTTCCTCCAGTGATTTCAGCTGCTGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGATATCCACGCG‑AAAGACTAAAACAGCTGGCAGGAAgcggcgcaa                                                                                                                                                                <  2:441375/233‑4 (MQ=255)                                                                                  gCCGACAGCGATGCGGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGt                                                                                                                                                           <  1:350867/249‑1 (MQ=255)                                                                                          agaggCTGACGAAAGCAGTTTCTGCTGAGACCGAAGGTTGTGGCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATGGATTACAGCGTCGGCTTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGAAGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGCGA‑CAGCGGCa                                                                                                                                                <  2:97849/247‑1 (MQ=255)                                                                                             tGCTGAAGAAAGCAGTCTATGCTGAAACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGAAGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTGACCACGCGGAAAG‑CGAAAACAGCTGGGAGGGAGCGGCGGAGACGGGGA‑AAGCAGCAAtt                                                                                                                                             >  1:181723/1‑251 (MQ=255)                                                                                                  aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACt                                                                                                                                        <  2:81775/251‑1 (MQ=255)                                                                                                            tCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGATATCCAAGCG‑ACAGCCCAAAACCACCGGCAGGCACCGGCGCAAGCGGGGAGAAG‑GGCAATTAAACTTTCCGAtca                                                                                                                               >  1:364682/1‑248 (MQ=255)                                                                                                                 gCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CAcc                                                                                                                        <  1:437241/251‑1 (MQ=255)                                                                                                                        ccGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGc                                                                                                                 <  1:502749/251‑1 (MQ=255)                                                                                                                            agtttttGTCAGGGAAAAGGATACCCGCGCGCTGGGGAGGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCaca‑c                                                                                                            <  2:366105/248‑1 (MQ=255)                                                                                                                              ggtttGGCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCa                                                                                                          <  1:496009/249‑1 (MQ=255)                                                                                                                                   gTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTg                                                                                                      <  1:199447/250‑1 (MQ=255)                                                                                                                                           aaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTGTCCACCCG‑GAAGACGAAAAAAGCTGGGAGGAAAAGGCGAAAACGGGGA‑AAGGCGCAATTAAAATTTCCGATTAAA‑CACCGGaaaagacacaaaattgctattccc                                                                                              >  1:503623/1‑230 (MQ=255)                                                                                                                                                                gATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGA‑AAGCGGCAAATAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAAATGTCTATCGCCGTATCAGCGAATAACggg                                                                           >  1:424322/1‑247 (MQ=255)                                                                                                                                                                                     cgtcggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAttt                                                   <  1:120460/241‑1 (MQ=255)                                                                                                                                                                                       gcgggcgtcGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACT‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGa                                                 <  1:43828/243‑1 (MQ=255)                                                                                                                                                                                        gggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGTTATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGTCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAGATTTGaa                                                <  2:86294/244‑1 (MQ=255)                                                                                                                                                                                        cggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGGGGGGA‑AAGAGGCAATTAAACCTGCCGATAAAA‑CACCGGAAAG‑ACAACCAACTGTCTATCGCCGCATCACAGAAACACGGGATACTGATCTGATCATTTAAAtttta                                                 >  2:73279/8‑248 (MQ=255)                                                                                                                                                                                            cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCa                                             <  1:427915/247‑1 (MQ=255)                                                                                                                                                                                            cgtcGTCAGGCAGCTGACACCATTGATTACAGCGTCGGCGTTACTGATAGGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACT‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCAc                                            <  2:205469/248‑1 (MQ=255)                                                                                                                                                                                                gtcgggCACCTGACACCATTGATTACAGCGTCGGCTGTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCCCTGAg                                        <  2:53276/247‑1 (MQ=255)                                                                                                                                                                                                           gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAAC‑a                            <  1:479376/251‑1 (MQ=255)                                                                                                                                                                                                              cccatttatgaCAGCGTCTGCTTTACTGATCTGGCGGTTTTGGGCGCCCAGGTAGACGGGCAGCGACCGCTGGCTGTTATCCACGCG‑AAAGACGAAAACTTCTGTCAGGAAGCGGCGAAAGAGGTGC‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCACCGTATCAGCGAATAACTGTATACTGATCTGATCATTTAAATGTGAAGCACTGCGTACGGAGCAC‑atat                         <  2:167956/241‑1 (MQ=255)                                                                                                                                                                                                               ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCAGCG‑AAAGACGAAAACAGCCGGCAGGAAGCGGCGGAAGCGGGGG‑ACGAGGCCATTAAACTTGCCGCTAAAA‑CACCCGGAAGCACA‑CCAACTGGCTATTGCCCTATCAGCGAAAAACGGGATACTGATCTGATTCTTTAAATTTTAAGCACTgtggaccgggaaaatatt                         >  2:175285/1‑234 (MQ=255)                                                                                                                                                                                                                      gtACAGCGTCGGTTTCACTGATATGGCGGGGTGGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAGCCACGCG‑AAAGACGAATACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CAGCGGAAAGCACA‑CCAACTGTCTATCGTCGTATCAGCGAATAACGGTATTCTGATCTGATCATTTAAATTTGAAGCACTGAGGACGTAGAAC‑ATCTGAAACGTg                 <  2:86449/250‑1 (MQ=255)                                                                                                                                                                                                                       tACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGAAGGTCAGCGTCCGCTGGCGGTTTTCCAAGCG‑AAAAACGAAAACAACTGGCAGGAAGCGGGGGAAGCGGGGA‑AAGAGGCAATTAAACTTTCCAATAAAA‑GACCGGAAAACCCA‑CCAACAGTCTATCGCCGTATAAGCGAATAACGGGCTACTGATCTTATCATTTAAATTTGTAGAACTgcggacagacaacataagtaaagggc                 >  1:373706/1‑226 (MQ=255)                                                                                                                                                                                                                            ggTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAAC‑ATATGAAACGTGCATtta           <  1:158341/250‑1 (MQ=255)                                                                                                                                                                                                                                   tttACTGCTATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGCTATCCACCCG‑AAAGACGGACACAGCTCGCAGGAAGCGGCGCAAGCGGTGA‑CCGCGGCCATTAAACATTCCCATAAAGACCCCCGAA‑GCACC‑CCAACTGTCTATCCCCGTTTTAGCGAATAACGGGGAACTGATTCGGTCATTTAAATTTTAAGCACCGACGACGGGGAAC‑AAATGAAAAAGGCATTTTTTATggg    >  2:46747/1‑250 (MQ=255)                                                                                                                                                                                                                                     tACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAAAAGCTGGCAGGGAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACCGGCTATCGCCGTATCAGCGAATAACGGGATACTGATCTGATCATTTTAATTTGAAGCAATGCGGACGGGGAACCTTAT‑AAAAGTGCATTTTTTATGGTgg  >  1:381693/1‑250 (MQ=255)                                                                                                                                                                                                                                                     cTGGGCGACCAGGTATACGGTCAGCGTCCGCTGGCGGTTAACCACGCG‑AAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTa                                                                                                                                            <  1:245225/100‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                TAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCG‑AAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGA‑AAGCGGCAATTAAACTTGCCGATAAAG‑CACCGGAAAGCACA‑CCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAAC‑ATATGAAACGTGCATTTATTATGGTGC  >  NC_000913/4619149‑4619627

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