Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,004	T→A	100%	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%	172.8 / NA	60	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (29/31);  total (29/31)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116860‑4117147                                                                                                                                                   |                                                                                                                                                aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCGGACAA‑TAG‑TGGTTTCAAGCTGGTTCGTAATATCGTATGCTGCAATTTGAACGGa                                                                                                                                             >  2:529337/1‑149 (MQ=255) aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGa                                                                                                                                             >  1:716566/1‑149 (MQ=255)     cctCAGGTTCTAAATGGCCGGTGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTg                                                                                                                                         <  2:455698/147‑1 (MQ=255)       tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGat                                                                                                                                       >  1:514322/1‑149 (MQ=255)        cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGCCGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGata                                                                                                                                      >  1:239139/1‑149 (MQ=255)         aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGatat                                                                                                                                     <  2:302890/149‑1 (MQ=255)           gTTTTAAATGCTCGCAGATCTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGTCAA‑TCG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTTCAATTTGAACGGAACTGCTATCg                                                                                                                                   >  2:580772/1‑149 (MQ=255)                aaaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCttt                                                                                                                              <  2:745733/149‑1 (MQ=255)                  aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAACCCGACGGCGTCTCG‑‑GGTTCGTAATACCGATCGCTGCAATTTGAACGGAACTGGTATCGGCTTTCg                                                                                                                            >  1:411893/1‑149 (MQ=255)                  aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCg                                                                                                                            >  1:148874/1‑149 (MQ=255)                  aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCg                                                                                                                            >  1:48445/1‑149 (MQ=255)                  aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGGTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGGACGGGACTGCTATCGGCTTTCg                                                                                                                            >  1:433/1‑149 (MQ=255)                    gCTCCCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGcc                                                                                                                          <  1:346614/149‑1 (MQ=255)                     cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCa                                                                                                                         <  1:531024/149‑1 (MQ=255)                        gCATATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTAGAGATAGGCTTGCCGTTTTCTTTTTCCCAGACTA‑AAG‑TGGTTTCAGGCTGTTCCGTAATACCGCGAGCTGCCATTTGAACGGAACTGATATCGGCTTTCGACAGCa                                                                                                                      <  2:45475/149‑1 (MQ=255)                              ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTTTTTTTTTT‑TTG‑TGGTTTCACGCTGGTTCGTAATACCGGTGGCTGCAATTGGAACGGAACTGATATCGGCTTTCGCCTGTACTTCtc                                                                                                                >  2:581197/1‑148 (MQ=255)                                        cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                                                                                      >  2:20727/1‑149 (MQ=255)                                          acgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATGTCGGCTTTCGCCCGCACTTCTACCAGCGTGGAGc                                                                                                    >  2:570357/1‑149 (MQ=255)                                            ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGGGCtt                                                                                                  >  1:521009/1‑149 (MQ=255)                                                cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTggg                                                                                              <  2:580951/149‑1 (MQ=255)                                                 tGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑GGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGt                                                                                             >  1:354624/1‑149 (MQ=255)                                                   ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           >  1:160410/1‑149 (MQ=255)                                                   ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           <  1:158735/149‑1 (MQ=255)                                                   ccAGACAATGGCCTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                           <  2:486811/149‑1 (MQ=255)                                                    cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                                          >  1:658237/1‑149 (MQ=255)                                                       aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCTATAGCTGCAATTTGAACGGTACTGATATAGGACTAATTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCa                                                                                       <  1:206651/149‑1 (MQ=255)                                                       aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTCCCAGGCTGGAGCCTTGGGTGGCCCa                                                                                       >  2:97010/1‑149 (MQ=255)                                                         aagtgCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                                     <  2:126274/145‑1 (MQ=255)                                                         aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTATACCAGCGTGGAGCTGTGGGTGGCCCGGa                                                                                     >  1:261403/1‑149 (MQ=255)                                                                   tagatacGCTTGCCGGTTTCTTTTTCCCATACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                           <  1:616351/149‑1 (MQ=255)                                                                    agttagGCTTGGCGGTTTCTTTTTCCCAGACGA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGCTCTGATATCGGCTTTCGCCAGCACTTCTGCCAGCGTGGAGCTTTGGGTGGGCCAGATTTCCATTggg                                                                          <  2:395758/146‑1 (MQ=255)                                                                         ggCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCGGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg                                                                     <  2:283592/149‑1 (MQ=255)                                                                              gCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                                                <  1:61644/149‑1 (MQ=255)                                                                                 ggTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAACTGCAATTTGAAAGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGGGCCTTGGGTGGCCCAGATTTCCCTTGGGTCGGGTTCTAccc                                                             >  1:372069/1‑149 (MQ=255)                                                                                    ttCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTGGGTCGGGTTCTAccccgc                                                          >  2:195053/1‑146 (MQ=255)                                                                                     tCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAcc                                                         >  1:509876/1‑149 (MQ=255)                                                                                       ttttcCCCAGGCAA‑GAG‑GGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACGTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                       <  1:276460/149‑1 (MQ=255)                                                                                          ttCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                    >  1:359664/1‑149 (MQ=255)                                                                                              cAGTCAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                <  1:63438/149‑1 (MQ=255)                                                                                              cAGACAA‑TTG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTCCCATCGGGTCTTGTCCTACCCAACCTGCTTTTgg                                                <  2:716566/149‑1 (MQ=255)                                                                                              cAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                <  2:160410/149‑1 (MQ=255)                                                                                              cAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                >  2:802524/1‑149 (MQ=255)                                                                                              cAGACAA‑TAG‑TGGTTTCACGCCGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                <  1:81977/149‑1 (MQ=255)                                                                                                      ctgatGGTTTCACGATGGATTGTAATTCCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                          <  2:215813/145‑1 (MQ=255)                                                                                                        g‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTg                                       <  2:579637/149‑1 (MQ=255)                                                                                                                cACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGGGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCCCACAt                                >  2:427156/1‑149 (MQ=255)                                                                                                                    cTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc                            <  2:48445/149‑1 (MQ=255)                                                                                                                     ttgTTCGTAATACCGTTAGCTGCAATTTGAACGGAACTGGTTTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCAAATTcgcg                           <  2:683206/147‑1 (MQ=255)                                                                                                                     tGGTTCGTACGACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTCCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgcg                           <  2:411893/149‑1 (MQ=255)                                                                                                                       gTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt                         <  2:144067/149‑1 (MQ=255)                                                                                                                         tCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGAGCTTTGGGTGGCCGAGATTTCCATTGGGGCGTGTTCTACCCAGCCTGGGGTTGGATGGATTTGCTCAAATTCGGGCTGc                       >  2:305645/1‑149 (MQ=255)                                                                                                                          cGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTATACCAACTAATGTTTGGGGTAGATTTGCTCAAATTCGCGCTGCg                      <  1:847805/149‑1 (MQ=255)                                                                                                                           gTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                     <  2:354624/149‑1 (MQ=255)                                                                                                                            tAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTCGAATTCGCGCTGCGac                    >  1:723085/1‑149 (MQ=255)                                                                                                                                 ccGATAGCTGGAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGTTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt               <  1:476124/149‑1 (MQ=255)                                                                                                                                   gATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCGGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAGTTGCTCAAATCCGCGCTGCGCCACGCtta             >  2:344424/1‑147 (MQ=255)                                                                                                                                     tAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATg           >  1:57189/1‑149 (MQ=255)                                                                                                                                     tAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATg           <  2:131108/149‑1 (MQ=255)                                                                                                                                        cTGCAATTTGAACGGAACTGATATCGGCCTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata        <  2:106814/149‑1 (MQ=255)                                                                                                                                              tttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa  >  2:53508/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                AAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAA‑TAG‑TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCA  >  NC_000913/4116860‑4117147

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