Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	113,695	G→A	W84* (TGG→TGA)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	113,695	0	G	A	100.0%	133.3 / NA	42	W84* (TGG→TGA)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/22);  total (20/22)

AAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAG  >  NC_000913/113454‑113941                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       aagaagATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACGGCTGTGCATAAACACTATTCTGGTGAAGAGTGACAATCGttt                                                                                                                                                                                                                                               <  2:520279/251‑1 (MQ=255)   gaagaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTAt                                                                                                                                                                                                                                             >  2:513784/1‑251 (MQ=255)     agaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATca                                                                                                                                                                                                                                           <  1:184981/251‑1 (MQ=255)     agaTCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATc                                                                                                                                                                                                                                            >  2:559168/1‑250 (MQ=255)                        aaaGACGTTCTCATCCGCCCTAAACGCTCCACTCGTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGTAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGTCTTCGTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGAtgt                                                                                                                                                                                                                        <  1:269795/251‑1 (MQ=255)                                         cccTAAACGCTCCACTCTTAAACGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATgg                                                                                                                                                                                                       >  1:161265/1‑251 (MQ=255)                                         cccTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATgg                                                                                                                                                                                                       <  1:158136/251‑1 (MQ=255)                                                      cctctTAAAAGCCGTTCCGATGTTGGACTTGAACGTCAAGTCACCTGCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTATGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACGCTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTAc                                                                                                                                                                                          <  2:15846/250‑1 (MQ=255)                                                                  cGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCgg                                                                                                                                                                               >  2:199344/1‑250 (MQ=255)                                                                   gTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGAt                                                                                                                                                                             >  2:8811/1‑251 (MQ=255)                                                                        gATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCAGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCg                                                                                                                                                                         >  2:147852/1‑250 (MQ=255)                                                                           gTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGaaa                                                                                                                                                                      <  1:564693/250‑1 (MQ=255)                                                                              gAACTGTAACGTCAATTCACCTTCAATCATTCAGGTCAGAGCTTGTACGGCGTGCCGATTATCTCCGCAAATATGGACTCCGTATGCACATTTTCTATGGCCTCTGCGCTGTCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCTTTTATCAACAATTCTTCCGCTTATGTGCTGAAACATGTGATTGTTTCAACCGGTACGTCTGATGCCGAtgtctaaagaa                                                                                                                                                                   <  2:85242/250‑11 (MQ=255)                                                                                        gTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGAtt                                                                                                                                                        >  2:341979/1‑251 (MQ=255)                                                                                        gTCAAATCACCCTCGAACAATCAGGTCAGGGCTGGTCCTGGGTTCCGATTAACGCCGCAAAAATTGACATCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTGTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATTAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTGTACCGGTACGTCTGATGCGGAGTTCGAAAAAACTAAACAGAtt                                                                                                                                                        <  2:34229/251‑1 (MQ=255)                                                                                                                   agagCTGGTCCGGCGTGCCGATTATCGCCGCTAATATGGACCCCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATCTTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTc                                                                                                                             >  2:383036/1‑251 (MQ=255)                                                                                                                         ggTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGttt                                                                                                                         <  1:513784/249‑1 (MQ=255)                                                                                                                             cGGCGTGCCGATTATCGACGCAAATATGGACACCGTAGGCAAATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGCTGTGCTGAAACATGTGATGGTTTCTAACGGTACGTCTGAGGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTg                                                                                                                   >  2:103430/1‑251 (MQ=255)                                                                                                                                    ccGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGc                                                                                                            <  1:424061/251‑1 (MQ=255)                                                                                                                                         tATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGactg                                                                                                       >  2:464013/1‑248 (MQ=255)                                                                                                                                         tATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATg                                                                                                       >  1:77862/1‑251 (MQ=255)                                                                                                                                               cacAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGAGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTAtt                                                                                                 >  1:192886/3‑251 (MQ=255)                                                                                                                                                gcAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTc                                                                                                <  2:424140/251‑1 (MQ=255)                                                                                                                                                 cAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGGTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGGTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTATCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTcc                                                                                               <  2:45931/251‑1 (MQ=255)                                                                                                                                                    atatGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGa                                                                                             >  1:264425/1‑250 (MQ=255)                                                                                                                                                          acacCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGGTATGCCGAACTCTTc                                                                                      <  2:218805/251‑1 (MQ=255)                                                                                                                                                              cGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGc                                                                                  <  1:559168/251‑1 (MQ=255)                                                                                                                                                                taggcccgttttctctggactaagcgcTGGCTTCTTTTGATATTTTGACTGCTGTGCATCAACACTATTCTGTCGAAGCGTGACAAGCGTTTATCAACAATTCTTCCGATGATGTGCTGCAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTCAACAGATTCTCGACCTGAACCCGGCATTAACCTTCGTTTGTATTGACGTGGCGACTGGATATTCCGAACACTTCGTGc                                                                                  >  2:240237/24‑249 (MQ=255)                                                                                                                                                                tAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGc                                                                                  >  1:364775/1‑249 (MQ=255)                                                                                                                                                                    cacaTTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGCAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTct                                                                            >  2:213128/1‑250 (MQ=255)                                                                                                                                                                              aTGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACTCTTCGTGCAGTTCGTTGCGAAAgc                                                                  <  1:431498/251‑1 (MQ=255)                                                                                                                                                                                ggCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCGGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAgcgc                                                                <  1:501936/251‑1 (MQ=255)                                                                                                                                                                                             gCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGcc                                                   <  1:282636/251‑1 (MQ=255)                                                                                                                                                                                                 cttTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGAcc                                               >  1:174379/1‑251 (MQ=255)                                                                                                                                                                                                      gATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAAc                                          <  1:464013/251‑1 (MQ=255)                                                                                                                                                                                                                ctgctgTGCATAAACACTATTCTGTCGAAGAGAGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCt                                >  1:445732/1‑251 (MQ=255)                                                                                                                                                                                                                        cATAAACACTATTCTGTCGAAGGGTGACAAGCGTTTATCAACAATTCTTCCGCTGCTGTGCTGAAACATGTGTTGGTTTCTACCGGTACGTCTGTTGAGGCTTTCGAAAAAGCTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGATTGTATTGACGTGGCGAATGGTTATTCAGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAACACCATTTGTGCTGGTAACGt                        >  2:312470/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTg                <  1:505018/249‑1 (MQ=255)                                                                                                                                                                                                                                    tCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTg             <  1:362729/250‑1 (MQ=255)                                                                                                                                                                                                                                     cTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGa            >  1:585112/1‑250 (MQ=255)                                                                                                                                                                                                                                           aaGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgtgt     <  2:174379/251‑1 (MQ=255)                                                                                                                                                                                                                                              agTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTgag  <  2:163724/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       AAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAG  >  NC_000913/113454‑113941

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