Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,004	T→A	D73V (GAT→GTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,004	0	T	A	100.0%	212.6 / NA	70	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/40);  total (30/40)

AACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116861‑4117147                                                                                                                                                |                                                                                                                                                aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTGAGACCGATAGCTGCAATTTGAACGGaa                                                                                                                                            >  1:540394/1‑149 (MQ=255)  aCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAAGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAAc                                                                                                                                           >  2:184107/1‑149 (MQ=255)    cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTg                                                                                                                                         <  2:391833/149‑1 (MQ=255)        aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                     <  1:1142428/149‑1 (MQ=255)          gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGCCAATAGGGGTTTCACGCTGGTTCGTACTATCGCGAGCTGGACTTTGAACGGACCTAATATCg                                                                                                                                   >  1:1154619/1‑149 (MQ=255)           ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCgg                                                                                                                                  >  1:1086627/1‑149 (MQ=255)                aaTGCTCGCAGATTTCTGCGGTAGGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTc                                                                                                                             <  1:30830/149‑1 (MQ=255)                aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAAGTGATATCGGCTTTc                                                                                                                             >  2:899126/1‑149 (MQ=255)                   gCTGGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGGTATAGATAGGCGTGCCGGTTTCTTTTTCCCCGACAATAGTGGTTTCACGCTGGTTCGTACTACCGATAGCTACAATTTGAACGGAACTGATATCGGCTTTCGcc                                                                                                                          <  2:31736/149‑1 (MQ=255)                   gCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGcc                                                                                                                          >  1:1185865/1‑149 (MQ=255)                     tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAg                                                                                                                        <  2:1110416/149‑1 (MQ=255)                            tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                                                                                 <  1:977516/149‑1 (MQ=255)                            tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                                                                                 <  2:721769/149‑1 (MQ=255)                                tGCGGTACGACGGCACTGCCAGACGATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGAGATCGGCTTTCGCCAGCACTTCTACCa                                                                                                             >  1:796650/1‑149 (MQ=255)                                 gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCttctgg                                                                                                            >  2:1048665/1‑144 (MQ=255)                                    gtgcgacgGCACTGCCAGACAATGGCTTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                                                                                         <  1:1007271/146‑1 (MQ=255)                                     tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCGTGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCTGCTTTCGCCAGCACTTCTACCAGCTTg                                                                                                        <  2:446240/149‑1 (MQ=255)                                        gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                                                                                     <  2:523322/149‑1 (MQ=255)                                        gacgGCACTGCAAGACAATGGCGTTATCGATAGGCTTGAAGGTTTCTTGTTCCCAGACAATAGTGGTTTAACACTTTTTCGAATTGCATTGAGCGGCAATTTGAACGGAACTGATCTCGGCTTTCGCCACGACTTCTACCAGCGTGGAg                                                                                                     <  1:911568/149‑1 (MQ=255)                                          cgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGCGCt                                                                                                   >  1:717000/1‑149 (MQ=255)                                                tGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGACTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGt                                                                                             >  2:350317/1‑149 (MQ=255)                                                  ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           >  1:995870/1‑149 (MQ=255)                                                      aCAATGGCGTTATAGTTAGGCTTGCCGGTTTCCTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTAACATCGTCTATCTTTGGGTGGCCCa                                                                                       <  2:148161/149‑1 (MQ=255)                                                         aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATCCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGAt                                                                                    >  1:742701/1‑149 (MQ=255)                                                           ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                     <  1:63567/98‑1 (MQ=255)                                                           ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                     >  2:63567/1‑98 (MQ=255)                                                                tataGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                             <  1:454059/149‑1 (MQ=255)                                                                  tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTgg                                                                           <  1:898152/149‑1 (MQ=255)                                                                  tagatagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                           <  1:470764/149‑1 (MQ=255)                                                                    gatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGGTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                         <  1:579389/149‑1 (MQ=255)                                                                       agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                                      <  1:1066017/149‑1 (MQ=255)                                                                         gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                    <  2:639393/149‑1 (MQ=255)                                                                         gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                    >  1:250358/1‑149 (MQ=255)                                                                               cGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAcc                                                              >  1:1063177/1‑149 (MQ=255)                                                                                ggTTGCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTGCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAccc                                                             <  1:378923/149‑1 (MQ=255)                                                                                   ttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                                          >  1:819424/1‑149 (MQ=255)                                                                                    tCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCGGATTTCCATTGGGTCGTGTTATACCCCAcc                                                         >  1:505553/1‑149 (MQ=255)                                                                                    tCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAcc                                                         <  1:65593/149‑1 (MQ=255)                                                                                      tttttCCCAGACACTAGTGGGTTCACGCTGGTTCGTAATGCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                       <  1:883149/149‑1 (MQ=255)                                                                                        tttCCCAGACCATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTATCAGCGTGGATCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGt                                                     <  1:969678/149‑1 (MQ=255)                                                                                         ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTATCCAACCTGGtt                                                    <  2:633869/149‑1 (MQ=255)                                                                                           cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                  <  1:1158841/149‑1 (MQ=255)                                                                                           cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                  <  2:1185865/149‑1 (MQ=255)                                                                                           cccAGACAATAGTGGTTGCACGCTGGTTCGTAATCCCGATAGCTGCAATTTGGACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                  <  1:237234/149‑1 (MQ=255)                                                                                            ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                 >  1:434046/1‑149 (MQ=255)                                                                                             ccgACAATTGTGGTTTCACGCTGGTTCGTAATTCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTgg                                                <  2:769397/147‑1 (MQ=255)                                                                                             cAGCCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCGGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                <  1:20001/149‑1 (MQ=255)                                                                                             cAGACTATCGCGGTTCCTCGCGTGTTCGTAATACCGGTAGCTGCAATTCGCACGGACCTGAGATCGGCTTGCCCCAGCACTTCTGCCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATGGGGTCTTGTTCTACCCAACCTGGTTTTgg                                                <  2:451998/149‑1 (MQ=255)                                                                                             cAGACAATAGTGGTTTCACGCTGGGTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTCTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                <  2:519109/149‑1 (MQ=255)                                                                                                   aTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                          >  2:999692/1‑149 (MQ=255)                                                                                                   aTAGGGGTTTCGCGCTGGGTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                          <  2:1116146/149‑1 (MQ=255)                                                                                                    tAGTGGTTTCACGCTGGTTCGTAATACCTATAGCTGCAATTTGAACGGAACTGATATCGCATTTAGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAACTGGGTGTGGGTAgcat                                         >  2:552895/1‑146 (MQ=255)                                                                                                         gTTTCACGCTGGTTCGTAATCCCGATAGATGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAATTCCCTTAGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGcgg                                    >  2:554530/1‑147 (MQ=255)                                                                                                          tttCACGCCGGTTCGTAATACCGATAGCTGCAATTTGAACGGGACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGCCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa                                   <  1:335956/149‑1 (MQ=255)                                                                                                              gcGCTGGTTCCTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCTGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTCTGGGTAGATTTGCTCAAAtt                               <  2:622875/148‑1 (MQ=255)                                                                                                              aCGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                               >  2:1122252/1‑149 (MQ=255)                                                                                                                     ttCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTGCTACCCAACCTGGTTTTGGGTCGATTTGCTCAAATTCGCGCTg                        >  2:642606/1‑149 (MQ=255)                                                                                                                     ttCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGACTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCATGTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                        <  1:350317/149‑1 (MQ=255)                                                                                                                         tAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                    <  1:125888/149‑1 (MQ=255)                                                                                                                         tAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCCCACATTCGCGCGTCGac                    >  2:633757/1‑149 (MQ=255)                                                                                                                          aaTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                   <  2:1003012/149‑1 (MQ=255)                                                                                                                           aTACCGAGAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTGGGTTTGCTCAAATTCTCGCTGCGaccc                  >  2:360814/1‑147 (MQ=255)                                                                                                                               cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGCCCCGCTa              >  1:23699/1‑149 (MQ=255)                                                                                                                               cGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTa              >  1:198482/1‑149 (MQ=255)                                                                                                                                 aTAGCTGCAATTTGAACGGAACTGATTTCGGCTTTCGCCAGCACTTCTACCAGCGTGGGGCTTTGGGTGGCCCAGCTTTCCGTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt            <  1:485469/149‑1 (MQ=255)                                                                                                                                      tGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTCGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat       >  1:144527/1‑149 (MQ=255)                                                                                                                                      tGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat       >  1:666030/1‑149 (MQ=255)                                                                                                                                       gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt      <  1:484746/149‑1 (MQ=255)                                                                                                                                        cAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTg     <  1:721102/149‑1 (MQ=255)                                                                                                                                           tttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATGTTGGCa  >  2:666654/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                AACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116861‑4117147

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