Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,018	A→G	T264A (ACG→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,018	0	A	G	100.0%	238.3 / NA	64	T264A (ACG→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/37);  total (27/37)

GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAG  >  NC_000913/4618871‑4619159                                                                                                                                                    |                                                                                                                                              gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgg                                                                                                                                            <  1:980086/151‑1 (MQ=255) gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgg                                                                                                                                            >  1:767299/1‑151 (MQ=255)     tctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgggcg                                                                                                                                         >  1:44624/1‑148 (MQ=255)      ctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTg                                                                                                                                         <  1:824969/149‑1 (MQ=255)      ctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTg                                                                                                                                         <  2:575075/149‑1 (MQ=255)      ctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTAGCGTGCGCACAACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGAAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGt                                                                                                                                          >  1:16208/1‑148 (MQ=255)          gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATc                                                                                                                                   <  2:224434/151‑1 (MQ=255)           aaGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAata                                                                                                                                     >  2:846710/1‑148 (MQ=255)                   gCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgc                                                                                                                          <  1:382578/151‑1 (MQ=255)                     cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGCACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAgtatacttaccc                                                                                                                            >  2:622093/1‑139 (MQ=255)                       aaGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCt                                                                                                                      <  1:221051/151‑1 (MQ=255)                       aaGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCt                                                                                                                      <  1:125163/151‑1 (MQ=255)                         gCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGt                                                                                                                    <  2:943252/151‑1 (MQ=255)                          cGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGtt                                                                                                                   >  2:871481/1‑151 (MQ=255)                               gttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGAt                                                                                                               <  2:694413/150‑1 (MQ=255)                               gttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGGAACCCGCGTCTGTTTGAt                                                                                                               >  2:193841/1‑150 (MQ=255)                                ttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGAt                                                                                                               >  1:518687/1‑149 (MQ=255)                                       ggCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGAt                                                                                                       <  2:524337/150‑1 (MQ=255)                                       ggCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATg                                                                                                      >  1:100251/1‑151 (MQ=255)                                       ggCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATg                                                                                                      <  1:433743/151‑1 (MQ=255)                                        gCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATgg                                                                                                     >  2:357011/1‑151 (MQ=255)                                         cTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATgg                                                                                                     >  1:252976/1‑150 (MQ=255)                                          tAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGGAACCCGCGTCTGTTTGATGTCACGATGgcg                                                                                                   >  2:339229/1‑151 (MQ=255)                                                   tGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCtgtg                                                                                              >  1:687557/1‑147 (MQ=255)                                                       gTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGgaga                                                                                      <  2:136500/151‑1 (MQ=255)                                                         gcgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATg                                                                                    <  2:492255/151‑1 (MQ=255)                                                            caccacCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTTTGGTTGATGGCACGATGGGCCTGGGCGTGGAGATg                                                                                    >  2:799513/1‑148 (MQ=255)                                                            caccacCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATg                                                                                    >  2:469878/1‑148 (MQ=255)                                                              ccaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGAt                                                                               <  1:67448/151‑1 (MQ=255)                                                              ccaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGAt                                                                               <  1:426687/151‑1 (MQ=255)                                                                   gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCg                                                                          <  1:138434/151‑1 (MQ=255)                                                                       tgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCa                                                                       >  1:502055/1‑150 (MQ=255)                                                                         ctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCaaa                                                                     >  2:329978/1‑150 (MQ=255)                                                                          tCACCAACATGAATTAGGTACTGGCTTCCAGTGCAGGTAACGCGGTTGAAGTTCTTGAAGCGGTGCAGTTCCTGGCGGGTGAATCTCGTAGCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCa                                                                       <  2:892562/147‑1 (MQ=255)                                                                           cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACt                                                                   <  2:999835/150‑1 (MQ=255)                                                                           cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTg                                                                  <  2:65345/151‑1 (MQ=255)                                                                            aCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTgg                                                                 <  2:586095/151‑1 (MQ=255)                                                                             ccGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACt                                                                   <  1:989814/148‑1 (MQ=255)                                                                             ccGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGc                                                                <  1:118668/151‑1 (MQ=255)                                                                                 cATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGaaa                                                            <  1:117733/151‑1 (MQ=255)                                                                                   tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGaa                                                             >  1:879531/1‑148 (MQ=255)                                                                                     aaTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATg                                                        <  2:767214/151‑1 (MQ=255)                                                                                        cAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACg                                                     >  1:361666/1‑151 (MQ=255)                                                                                        cAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACg                                                     <  1:548959/151‑1 (MQ=255)                                                                                           gTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                  <  2:76331/151‑1 (MQ=255)                                                                                            tACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                  >  1:294476/1‑150 (MQ=255)                                                                                                ggCCTCCAGTGCGGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAg                                               <  2:566733/149‑1 (MQ=255)                                                                                                 gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgc                                            <  2:939446/151‑1 (MQ=255)                                                                                                    tCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgc                                          >  2:501789/1‑150 (MQ=255)                                                                                                     ccAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgcgc                                        >  1:578555/1‑151 (MQ=255)                                                                                                     ccAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgc                                          <  1:962570/149‑1 (MQ=255)                                                                                                         tGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGcgcgcgcg                                       <  1:447114/148‑1 (MQ=255)                                                                                                                aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAg                              <  1:109502/150‑1 (MQ=255)                                                                                                                 aCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGc                            <  1:526598/151‑1 (MQ=255)                                                                                                                    cgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCgg                          >  2:438713/1‑150 (MQ=255)                                                                                                                    cgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCggg                         >  2:933710/1‑150 (MQ=255)                                                                                                                     ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTg                        <  2:442060/151‑1 (MQ=255)                                                                                                                      gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGTGAAATTGCAGGCgg                          >  1:618393/1‑148 (MQ=255)                                                                                                                          aaGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGa                   <  1:958439/151‑1 (MQ=255)                                                                                                                            gTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACa                 <  2:6453/151‑1 (MQ=255)                                                                                                                                 tGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAAAgg             >  1:416702/1‑150 (MQ=255)                                                                                                                                    aGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTaaa         <  2:44622/151‑1 (MQ=255)                                                                                                                                           cAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGTACAACGGTAAAGCGGCAg  <  2:16205/151‑1 (MQ=255)                                                                                                                                           cAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCgg     >  1:360234/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                              GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAG  >  NC_000913/4618871‑4619159

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