Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,937	T→G	100%	K96Q (AAG→CAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,937	0	T	G	100.0%	339.7 / NA	93	K96Q (AAG→CAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (47/46);  total (47/46)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTC  >  NC_000913/4116798‑4117078                                                                                                                                            |                                                                                                                                              acttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTc                                                                                                                                      >  1:609668/1‑149 (MQ=255)    ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttt                                                                                                                                   <  3:705588/149‑1 (MQ=255)      ggcgCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttttt                                                                                                                                 <  1:622819/146‑1 (MQ=255)      ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttttt                                                                                                                                 <  2:321717/149‑1 (MQ=255)      ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttttt                                                                                                                                 <  3:744828/149‑1 (MQ=255)      ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttttt                                                                                                                                 <  1:599700/149‑1 (MQ=255)      ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCttttt                                                                                                                                 <  4:4106/149‑1 (MQ=255)       gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTc                                                                                                                                <  4:379382/149‑1 (MQ=255)        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTcc                                                                                                                               <  4:501097/149‑1 (MQ=255)          ccAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCCGACACTGGCGTTATCGATAGGCTGGCCGGTTTCTTTTTCCCa                                                                                                                             >  3:511971/1‑149 (MQ=255)            aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGa                                                                                                                           <  4:272793/149‑1 (MQ=255)                    aCGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTg                                                                                                                   <  4:453005/149‑1 (MQ=255)                    aCGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTg                                                                                                                   <  3:487568/149‑1 (MQ=255)                      gggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                 >  4:381621/1‑149 (MQ=255)                      gggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                 >  4:299595/1‑149 (MQ=255)                      gggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                 >  3:23642/1‑149 (MQ=255)                         tCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                                                                                              >  1:37857/1‑149 (MQ=255)                           aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                            >  1:17704/1‑149 (MQ=255)                           aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                                                                                            <  3:267030/149‑1 (MQ=255)                             tCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGc                                                                                                          <  2:477285/149‑1 (MQ=255)                              cACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCt                                                                                                         >  3:720944/1‑149 (MQ=255)                              cACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCt                                                                                                         >  3:371113/1‑149 (MQ=255)                               aCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTg                                                                                                        <  3:381621/149‑1 (MQ=255)                                  aGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGtt                                                                                                     >  1:676749/1‑149 (MQ=255)                                   gACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                                                                                                    >  3:723716/1‑149 (MQ=255)                                   gACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                                                                                                    <  4:499567/149‑1 (MQ=255)                                     ccGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                                                                                                  <  2:326181/149‑1 (MQ=255)                                        gTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAt                                                                                               >  3:625246/1‑149 (MQ=255)                                         tATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                                                                                              <  4:20397/149‑1 (MQ=255)                                           ttGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAcc                                                                                            >  4:599796/1‑149 (MQ=255)                                             gctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGAAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTTGTTTCACGCTGGTTCGTAATACCGa                                                                                          <  2:129848/149‑1 (MQ=255)                                              ctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAt                                                                                         >  3:715073/1‑149 (MQ=255)                                               tgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa                                                                                        <  4:85786/149‑1 (MQ=255)                                                gcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAg                                                                                       <  3:299595/149‑1 (MQ=255)                                                  ggATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCt                                                                                     >  4:483110/1‑149 (MQ=255)                                                        aaTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGAGAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                               >  4:518598/1‑149 (MQ=255)                                                        aaTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                               >  1:498580/1‑149 (MQ=255)                                                        aaTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                               >  4:564449/1‑149 (MQ=255)                                                         aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                              >  2:427906/1‑149 (MQ=255)                                                          tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTGATAGATAGGCTGGCCGGTATCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGTTGCAATTTg                                                                             >  3:355204/1‑149 (MQ=255)                                                            ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAt                                                                                         >  4:432457/1‑135 (MQ=255)                                                            ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAt                                                                                         <  3:432457/135‑1 (MQ=255)                                                             tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                          >  3:137117/1‑149 (MQ=255)                                                               tAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCgg                                                                        >  2:641929/1‑149 (MQ=255)                                                                aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGa                                                                       >  3:271151/1‑149 (MQ=255)                                                                 aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                      <  4:14130/149‑1 (MQ=255)                                                                 aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                      >  1:231293/1‑149 (MQ=255)                                                                   ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                    <  4:685379/149‑1 (MQ=255)                                                                    cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTg                                                                   <  2:498580/149‑1 (MQ=255)                                                                    cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTg                                                                   <  4:271151/149‑1 (MQ=255)                                                                      tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                 >  3:751234/1‑149 (MQ=255)                                                                        aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGatat                                                               >  2:102751/1‑149 (MQ=255)                                                                        aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGATTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCACTTTGATCGGAACTGatat                                                               >  3:485444/1‑149 (MQ=255)                                                                          gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCTCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                             >  2:296736/1‑149 (MQ=255)                                                                          gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                             >  1:468613/1‑149 (MQ=255)                                                                           ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                            >  4:755282/1‑149 (MQ=255)                                                                                aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                       <  2:137126/149‑1 (MQ=255)                                                                                aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                       <  2:487139/149‑1 (MQ=255)                                                                                        cAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTC‑TTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGCACTGATATCGGCTTTCGCCAGCACt                                              <  4:403765/149‑1 (MQ=255)                                                                                                  cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                     <  3:755282/149‑1 (MQ=255)                                                                                                        gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                               >  1:73343/1‑149 (MQ=255)                                                                                                            gCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACGATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                           <  1:134182/149‑1 (MQ=255)                                                                                                            gCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                           <  2:683657/149‑1 (MQ=255)                                                                                                             cACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                          <  4:715073/149‑1 (MQ=255)                                                                                                              aCTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                         <  2:200901/149‑1 (MQ=255)                                                                                                               cTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                        <  3:672062/149‑1 (MQ=255)                                                                                                                tGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGt                       >  3:7709/1‑149 (MQ=255)                                                                                                                 gCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                      <  2:676749/149‑1 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                               <  3:243749/91‑1 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt                                                                               >  4:243749/1‑91 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                      <  3:191908/116‑1 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                      >  4:191908/1‑116 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                     >  2:304081/1‑149 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                     >  2:384493/1‑149 (MQ=255)                                                                                                                  ccAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                     >  3:391491/1‑149 (MQ=255)                                                                                                                   cAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                    >  1:752814/1‑149 (MQ=255)                                                                                                                   cAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                    >  4:48776/1‑149 (MQ=255)                                                                                                                    aGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                   >  1:144935/1‑149 (MQ=255)                                                                                                                     gACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                  >  4:702570/1‑149 (MQ=255)                                                                                                                        aaTGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa               >  1:537633/1‑149 (MQ=255)                                                                                                                        aaTGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa               <  4:16299/149‑1 (MQ=255)                                                                                                                            gCGTTATAGATAGGCTGGACGGTTTCTTTTTCCCAGACAATAGTGGGTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc           <  4:720944/149‑1 (MQ=255)                                                                                                                               ttATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                              <  2:317217/127‑1 (MQ=255)                                                                                                                               ttATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa               >  1:319556/1‑142 (MQ=255)                                                                                                                               ttATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa               <  2:319556/142‑1 (MQ=255)                                                                                                                               ttATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAAAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                              >  1:317217/1‑127 (MQ=255)                                                                                                                                 atagatagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg      <  1:641929/149‑1 (MQ=255)                                                                                                                                  tagatagGCTGGCCGGTTTCTTTTTCTCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg     <  1:296736/149‑1 (MQ=255)                                                                                                                                  tagatagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg     <  2:537633/149‑1 (MQ=255)                                                                                                                                   agatagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg    <  2:73343/149‑1 (MQ=255)                                                                                                                                     atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc  >  1:466295/1‑149 (MQ=255)                                                                                                                                     atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc  >  1:748648/1‑149 (MQ=255)                                                                                                                                     atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTc  <  3:599796/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                              ACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTC  >  NC_000913/4116798‑4117078

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