Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,741,289	C→G	intergenic (+526/+227)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,741,289	0	C	G	100.0%	85.4 / NA	94	intergenic (+526/+227)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (49/45);  total (49/45)

TATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCT  >  AE015451/4741183‑4741388                                                                                                               |                                                                                                    tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                    >  1:5849181/1‑125 (MQ=255)   ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                                                  >  1:5849153/1‑125 (MQ=255)   ttGAACTGGGTATAGGACCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                                                  <  1:5848948/125‑1 (MQ=255)       aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                                              <  1:5848952/125‑1 (MQ=255)       aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                                              <  1:5848953/125‑1 (MQ=255)       aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                                              >  1:5849089/1‑125 (MQ=255)        cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                             <  2:5848955/125‑1 (MQ=255)        cTGGGTATAGGTCCACCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                             <  1:5848956/125‑1 (MQ=255)         tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                             <  1:5848957/124‑1 (MQ=255)         tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                            >  1:5849414/1‑125 (MQ=255)         tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                            >  1:5849231/1‑125 (MQ=255)          gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                                           <  2:5848958/125‑1 (MQ=255)          gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                                           >  2:5849168/1‑125 (MQ=255)          gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                                           >  2:5849214/1‑125 (MQ=255)           ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                             <  2:5848960/122‑1 (MQ=255)           ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTGGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                          <  1:5848959/125‑1 (MQ=255)            gTATCGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         <  2:5848964/125‑1 (MQ=255)            gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         <  1:5848961/125‑1 (MQ=255)            gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         <  2:5848963/125‑1 (MQ=255)            gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         >  1:5849352/1‑125 (MQ=255)            gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         >  1:5849312/1‑125 (MQ=255)             tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCTTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                                         >  1:5849349/1‑124 (MQ=255)              ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                       <  1:5848965/125‑1 (MQ=255)              ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                       <  2:5848966/125‑1 (MQ=255)              ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                        >  1:5849430/1‑124 (MQ=255)               taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGCACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                      <  1:5848967/125‑1 (MQ=255)                aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                     <  1:5848968/125‑1 (MQ=255)                aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGTCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                     <  2:5848969/125‑1 (MQ=255)                aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                     >  2:5849165/1‑125 (MQ=255)                 ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                                    <  2:5848971/125‑1 (MQ=255)                 ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                     <  2:5848970/124‑1 (MQ=255)                  gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGTCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                   >  1:5849377/1‑125 (MQ=255)                  gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                                    >  1:5849260/1‑124 (MQ=255)                  gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                   >  2:5848974/1‑125 (MQ=255)                   tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                                  >  2:5849319/1‑125 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                   <  2:5848996/123‑1 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                 <  1:5848973/125‑1 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                 <  1:5849061/125‑1 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                 <  2:5849003/125‑1 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                                  <  2:5848982/124‑1 (MQ=255)                    ccAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                                                   >  1:5848996/1‑123 (MQ=255)                     ccgCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                                               >  2:5849254/3‑110 (MQ=255)                     cAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCAGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                                >  1:5849370/1‑125 (MQ=255)                      aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca                                                               <  2:5849028/125‑1 (MQ=255)                      aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca                                                               <  1:5849059/125‑1 (MQ=255)                      aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                                >  2:5849393/1‑124 (MQ=255)                       gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                 <  1:5849056/122‑1 (MQ=255)                       gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                              <  1:5849051/125‑1 (MQ=255)                       gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                              >  2:5849322/1‑125 (MQ=255)                       gCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                              >  1:5849115/1‑125 (MQ=255)                        ctctGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                 <  1:5849005/121‑1 (MQ=255)                        ctctGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc                                                             >  1:5849193/1‑125 (MQ=255)                         tctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc                                                             >  2:5849327/1‑124 (MQ=255)                         tctGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc                                                             >  1:5849320/1‑124 (MQ=255)                          ctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                                           <  2:5849017/125‑1 (MQ=255)                          ctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                            <  1:5848989/124‑1 (MQ=255)                          ctGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                            <  1:5849057/124‑1 (MQ=255)                             gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                                                        <  1:5849032/125‑1 (MQ=255)                             gTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                                                        <  1:5849072/125‑1 (MQ=255)                              tCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                                                       <  1:5849088/125‑1 (MQ=255)                               cACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                                                      <  1:5848993/125‑1 (MQ=255)                               cACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGGACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                                                       <  1:5848984/124‑1 (MQ=255)                                aCCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAAGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg                                                          >  1:5849486/1‑120 (MQ=255)                                aCCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg                                                          >  1:5849387/1‑120 (MQ=255)                                 ccTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                                                    >  1:5849234/1‑125 (MQ=255)                                 ccTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCC‑GCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                                   >  1:5849459/1‑125 (MQ=255)                                  cTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                                   <  2:5849045/125‑1 (MQ=255)                                  cTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                                   >  2:5849178/1‑125 (MQ=255)                                                            acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc                            >  1:5849340/7‑124 (MQ=255)                                                            acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc                            >  1:5849407/7‑124 (MQ=255)                                                            acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCa                           >  2:5849157/7‑125 (MQ=255)                                                            acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCCAGTTAGGCCTTGTACTTGGACCCGGCCACCTc                            >  1:5849435/7‑124 (MQ=255)                                                             cttgcaTTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc                          <  2:5849083/120‑1 (MQ=255)                                                             cttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCAGGCCACCTc                            >  1:5849566/6‑123 (MQ=255)                                                              ttgcaTTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACt                         <  1:5849011/121‑1 (MQ=255)                                                              ttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACt                         >  2:5849240/5‑125 (MQ=255)                                                               tgcaTTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTg                        <  2:5848975/122‑1 (MQ=255)                                                                gcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGa                       >  1:5849039/3‑125 (MQ=255)                                                                  aTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAgg                     >  1:5849328/1‑125 (MQ=255)                                                                   ttCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAggg                    >  1:5849455/1‑125 (MQ=255)                                                                   ttCGTGAA‑TCCAGGCACATGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCATGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACACGGCCACCTCACTGAggg                    >  2:5849398/1‑125 (MQ=255)                                                                   ttCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAggg                    <  2:5849048/125‑1 (MQ=255)                                                                   ttCGTGAA‑TCCAGCTACAGGCCCCGGACACGTCCTCCTTCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTAGGACCCGGCCACCTCACt                         >  1:5849782/1‑120 (MQ=255)                                                                               aggcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCt              >  2:5849309/7‑120 (MQ=255)                                                                                ggcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATcc        >  2:5849275/6‑125 (MQ=255)                                                                                ggcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATc         <  1:5849064/119‑1 (MQ=255)                                                                                 gcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATccc       >  1:5849100/5‑125 (MQ=255)                                                                                   acaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATccc       <  1:5849047/121‑1 (MQ=255)                                                                                   acaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcg      >  2:5849258/3‑124 (MQ=255)                                                                                    caGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtc    >  1:5849718/2‑125 (MQ=255)                                                                                    caGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtc    >  1:5849497/2‑125 (MQ=255)                                                                                     aGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcc   >  2:5849256/1‑125 (MQ=255)                                                                                      ggCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct  >  1:5849403/1‑125 (MQ=255)                                                                                      ggCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct  <  1:5849042/125‑1 (MQ=255)                                                                                                               |                                                                                                    TATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCT  >  AE015451/4741183‑4741388

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