Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,068,287	G→A	F330F (TTC→TTT)	yihQ ←	sulfoquinovosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,068,287	0	G	A	100.0%	153.1 / NA	45	F330F (TTC→TTT)	yihQ	sulfoquinovosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

TGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGGAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTC  >  NC_000913/4068043‑4068514                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    ggggcgactACTCGCCGAATTCAGCCACCCGGCCGCCGCAGCCGAGTTCAATCATTTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGc                                                                                                                                                                                                                               <  2:114301/243‑1 (MQ=255)  gggCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACACGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCAGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGc                                                                                                                                                                                                                               >  1:157721/1‑250 (MQ=255)      aGATCCTCGCCGAACTCCGCCATCCAGCCGCCGCAGCCGAGTTAAATCATTTTCTTTTTGATCACTTCCTAGACCCGGGCGTAGGCTTCTGGATTAGTGAGTTCGCAAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGcc                                                                                                                                                                                                                          <  2:240427/251‑1 (MQ=255)        aTACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCcct                                                                                                                                                                                                                        >  1:91437/1‑251 (MQ=255)              gCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGAAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGAt                                                                                                                                                                                                                  >  1:224837/1‑251 (MQ=255)               ccGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGAtt                                                                                                                                                                                                                 <  1:124145/251‑1 (MQ=255)               ccGAAGCCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGAtt                                                                                                                                                                                                                 <  1:103775/251‑1 (MQ=255)                 aatgccgcccagccAGCCCCCGCAGCCGAGTTCACTCTTTTTTTTTTTGGTCACTGCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTGCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTcc                                                                                                                                                                                                               <  2:219187/239‑1 (MQ=255)                       aGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCt                                                                                                                                                                                                         >  1:382709/1‑251 (MQ=255)                                  cgccgcAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGa                                                                                                                                                                                               >  1:67219/1‑250 (MQ=255)                                      gcAGCCGAGTTCAATCATGTTCTTTTTTATCACTTCCTTGAACAAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCa                                                                                                                                                                                          <  2:328221/251‑1 (MQ=255)                                                                      cTTCCTTGAACCAGGCGTAGTCTTATAGCTTAGGGAAATCAACACCGCCGCCGTAAAACTCGCCACACTCCACCAGCTCGTCACCGCCCGCGCCATCTCTTGCCAGATAGCCGTGTTGTGCAGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGCTATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTCGTTTTGGCTGTTCAACTTcc                                                                                                                                                            >  2:167480/1‑249 (MQ=255)                                                                                ccAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATcc                                                                                                                                                >  1:282613/1‑250 (MQ=255)                                                                                 cgggCGTAGGCTTGTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGGCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAAGGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCAc                                                                                                                                               <  1:367397/249‑1 (MQ=255)                                                                                        aGGCTTCTGGATTAGTGAGATCGCCCACGCCGCCGTAAAACTCGCCAAACTCCACCAGCTAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATTCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTg                                                                                                                                        <  2:130105/251‑1 (MQ=255)                                                                                        aGGCTTCTGGATTAGAGAGATCGCAAACGCCGCCGTAAAACTAGGCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTCGCCAGCTAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGAAACCTCCGGGTTGATATAGGCCAGAAACTGCACGAGCTCCGGATTCCACTGCTTAATGAGTGAATCCAGTTGCGGGGAGTTTTCGTTATTCCACTTCCAGTTCGACATCACGCGtttt                                                                                                                                        >  2:344444/1‑250 (MQ=255)                                                                                          gCTTCTGGATTGGTGAGATCGACAACGCCTCCGTAAAACTCGCCAAACTCAACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATCGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTggc                                                                                                                                      >  2:133430/1‑249 (MQ=255)                                                                                          gCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTGTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGcc                                                                                                                                      <  1:145890/251‑1 (MQ=255)                                                                                               tGGATTAGTGAGATCGCCAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTCTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCaaa                                                                                                                                   >  2:4132/1‑249 (MQ=255)                                                                                                       tGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTg                                                                                                                                                                      <  2:9786/206‑1 (MQ=255)                                                                                                       tGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTg                                                                                                                                                                      >  1:9786/1‑206 (MQ=255)                                                                                                          gagcGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGGCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGa                                                                                                                      <  2:339553/248‑1 (MQ=255)                                                                                                          gaTCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGa                                                                                                                      <  1:200310/251‑1 (MQ=255)                                                                                                            ccGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTTGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAAt                                                                                                                    <  2:34086/250‑1 (MQ=255)                                                                                                                   ccccgccgTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGAc                                                                                                             <  2:185901/249‑1 (MQ=255)                                                                                                                       gccgTAAAACTCGCAGAACGCCACCAGATAGTCACCGCCATAGGCATCTTTTGCCAGATAGCCGGGTTGTGCCGCTTCTTCTCAGAGATCGTGATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCGCTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCCTTTGCCAAAAGAGCTTATACGAATACCGGACCAGt                                                                                                         <  2:252282/251‑1 (MQ=255)                                                                                                                          gTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCt                                                                                                      >  1:288575/1‑251 (MQ=255)                                                                                                                            aaaaCTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTgc                                                                                                    >  2:124210/1‑251 (MQ=255)                                                                                                                             aaaCTCGCCAAACTCCACCAGATAGTCACCGCCAGAAGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCCGTCCTgcg                                                                                                   >  2:325748/1‑251 (MQ=255)                                                                                                                               aCTCGCCAAACTCCACCAGAGAGTCACCGCCAGAGGCATCTTTTGCCATATAGCCGGGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCACCATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGGTTTCGCTGTTCCACTTCCAGTTCCGCATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGcc                                                                                                 <  2:118605/251‑1 (MQ=255)                                                                                                                                         cTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGCGGTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGTCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCACAAGAGGTCATACGAATACCGGGCCAGTCCTGCGCCCAGATGCCGt                                                                                       >  2:314880/1‑251 (MQ=255)                                                                                                                                                aGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCtt                                                                                <  1:299236/251‑1 (MQ=255)                                                                                                                                                    aGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACg                                                                            <  1:94258/251‑1 (MQ=255)                                                                                                                                                          cGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCcgcc                                                                      >  2:121890/1‑250 (MQ=255)                                                                                                                                                                agGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGc                                                                >  1:47931/1‑251 (MQ=255)                                                                                                                                                                 ggCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGCTCCACATCACTCGTTTGCCAAAAGAGGTCATACGAATTCCCGACCAGTCCTCCGCCCAGATGCCGTTGACCTTCACTCCCTCGTTACGCa                                                               <  2:242381/251‑1 (MQ=255)                                                                                                                                                                                    ccGTGTTGTGTCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGGTTGCCAAAAGAGGTCAGACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACTCATGGTGTCCAGTTCCTTCTg                                            <  2:45493/251‑1 (MQ=255)                                                                                                                                                                                    ccGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCCACATAAGGGTTGATCTAGGCCAGAAACTGCACGCTCTCCTGATTCCACTGCTTAATGCGTGCATCCTGTTGCGGGCAGCTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTTCCAAAAGAGGACATCCGAATACCGGACCCGTCCTGCCCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCCGTTtcttctt                                            >  1:17813/1‑250 (MQ=255)                                                                                                                                                                                      gtgtTGTGACGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCAACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCAACTTCAAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGCATACCGGACCAGTCCCGCGCCCAGATGCCGTTGCCCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTCCTGGc                                          >  2:118704/1‑251 (MQ=255)                                                                                                                                                                                          tgtgCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCAAGTTGCGGGGAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCacac                                      >  1:158068/1‑251 (MQ=255)                                                                                                                                                                                          tgtgCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACGGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACAGCGCGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCATATTCCGTTGACCTTCACGCCCGCGTTACGCATGTTGGCCAGTTGCTTTTGGcaca                                       <  2:17813/250‑1 (MQ=255)                                                                                                                                                                                                      tcGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTccc                             <  1:213732/248‑1 (MQ=255)                                                                                                                                                                                                      tcGCAGAGATCTTTATCGCTGGCAACATACGGGTTGAGATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCTCTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCAGGGTGTCCAGTTTCTTCTGGCCCACTTCCGTccc                             <  2:91437/248‑1 (MQ=255)                                                                                                                                                                                                       cgaagagATCTTTATCGATGGCAACATACAGGTTGAAATAGGTCAGAAACTGCACACTCTCCTGATTCCACTGCTTAATGCGTAAAGCCAGTTGCGGGTAGTTTACGCTGTTCCACTTCAAGTTCCACATCACGCATTTGCCAAAAGAGGTCAGACGAATACCGGACCTGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTCCGCATGGTGTCCAGTGTCTTCTGGCACGCTTCCATCCCGccc                         <  2:157721/248‑1 (MQ=255)                                                                                                                                                                                                                              aCATACGGGTTGATATAGGCCAGAAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTc  <  1:90704/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    TGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCGCCGTAAAACTCGCCAAACTCCACCAGATAGTCACCGCCAGAGGCATCTTTTGCCAGATAGCCGTGTTGTGCCGCTTCTTCGCAGAGATCTTTATCGCTGGCAACATACGGGTTGATATAGGCCAGGAACTGCACGCCCTCCTGATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTC  >  NC_000913/4068043‑4068514

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