Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,123	G→A	R34C (CGC→TGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,123	0	G	A	100.0%	310.9 / NA	95	R34C (CGC→TGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (47/48);  total (47/48)

GATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATA  >  NC_000913/4116989‑4117244                                                                                                                                       |                                                                                                                          gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGc                                                                                                                       >  2:671585/1‑139 (MQ=255) gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGc                                                                                                                       >  2:430919/1‑139 (MQ=255)  aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                      >  1:782509/1‑139 (MQ=255)  aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                      >  1:662080/1‑139 (MQ=255)   tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     <  2:410188/139‑1 (MQ=255)   tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     <  2:1193619/139‑1 (MQ=255)      cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGacac                                                                                                                  <  2:782509/139‑1 (MQ=255)       tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                 <  1:1140931/139‑1 (MQ=255)         cAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCt                                                                                                               >  1:530186/1‑139 (MQ=255)          aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              <  1:430919/139‑1 (MQ=255)               gATCGGAACTGATATCGGCTTTCGCCGGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGat                                                                                                         <  1:414341/139‑1 (MQ=255)                aTCTGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGata                                                                                                        <  1:836011/139‑1 (MQ=255)                      aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  2:500231/1‑122 (MQ=255)                      aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  1:500231/122‑1 (MQ=255)                          atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:1074182/1‑118 (MQ=255)                          atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:1074182/118‑1 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                      <  2:984528/110‑1 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCg                                                                                                                      >  1:984528/1‑110 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:255505/1‑113 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:255505/113‑1 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                         >  2:1087864/1‑139 (MQ=255)                               ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                         >  2:280530/1‑139 (MQ=255)                                    tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     >  1:1229046/1‑106 (MQ=255)                                    tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     <  2:1229046/106‑1 (MQ=255)                                                cTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                                                        >  2:147606/1‑139 (MQ=255)                                                cTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                                                        >  2:1232003/1‑139 (MQ=255)                                                      gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                 <  1:1121222/92‑1 (MQ=255)                                                      gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                 >  2:1121222/1‑92 (MQ=255)                                                       cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                 <  1:147606/139‑1 (MQ=255)                                                       cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                 <  1:671585/139‑1 (MQ=255)                                                        gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTaa                                                                                                             >  2:609824/1‑94 (MQ=255)                                                        gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTaa                                                                                                             <  1:609824/94‑1 (MQ=255)                                                        gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            <  1:1178760/111‑1 (MQ=255)                                                        gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            >  2:1178760/1‑111 (MQ=255)                                                         tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:822645/1‑87 (MQ=255)                                                         tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:822645/87‑1 (MQ=255)                                                           gAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACGCAACCTGGTATTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  2:799062/1‑85 (MQ=255)                                                           gAGCTTTAGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  1:799062/85‑1 (MQ=255)                                                                ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTgg                                                                                                    >  1:880477/1‑95 (MQ=255)                                                                ttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTgg                                                                                                    <  2:880477/95‑1 (MQ=255)                                                                      ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                  >  2:1234461/1‑91 (MQ=255)                                                                      ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                  <  1:1234461/91‑1 (MQ=255)                                                                       gCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg                                                 <  1:280530/139‑1 (MQ=255)                                                                         ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                               <  1:1087864/139‑1 (MQ=255)                                                                         ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                               <  2:530186/139‑1 (MQ=255)                                                                                  cATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     <  2:96213/60‑1 (MQ=255)                                                                                  cATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGa                                                                                                                     >  1:96213/1‑60 (MQ=255)                                                                                  cATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                       >  1:166522/1‑106 (MQ=255)                                                                                  cATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                       <  2:166522/106‑1 (MQ=255)                                                                                       ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                 <  2:434205/59‑1 (MQ=255)                                                                                       ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACg                                                                                                                 >  1:434205/1‑59 (MQ=255)                                                                                       ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            <  2:662791/80‑1 (MQ=255)                                                                                       ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            >  1:662791/1‑80 (MQ=255)                                                                                         tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttttt                               >  2:1138147/1‑139 (MQ=255)                                                                                         tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttttt                               <  2:1153566/139‑1 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  2:698067/1‑52 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:1059570/52‑1 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  1:698067/52‑1 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:1059570/1‑52 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:1209856/52‑1 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:1209856/1‑52 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              <  2:65343/57‑1 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              >  1:65343/1‑57 (MQ=255)                                                                                            tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                            >  2:381359/1‑139 (MQ=255)                                                                                             gtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGt                           <  2:645036/139‑1 (MQ=255)                                                                                                  aCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGt                      >  2:794678/1‑139 (MQ=255)                                                                                                      aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              <  1:703181/47‑1 (MQ=255)                                                                                                      aaCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              >  2:703181/1‑47 (MQ=255)                                                                                                       aCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                       <  1:890999/85‑1 (MQ=255)                                                                                                       aCCTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                       >  2:890999/1‑85 (MQ=255)                                                                                                        ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   >  1:356495/1‑40 (MQ=255)                                                                                                        ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGaca                                                                                                                   <  2:356495/40‑1 (MQ=255)                                                                                                        ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              <  1:1291285/45‑1 (MQ=255)                                                                                                        ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTa                                                                                                              >  2:1291285/1‑45 (MQ=255)                                                                                                        ccTGGTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                >  1:246430/1‑139 (MQ=255)                                                                                                            gTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATc                                                                                                <  2:1282049/55‑1 (MQ=255)                                                                                                            gTTTTGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATc                                                                                                >  1:1282049/1‑55 (MQ=255)                                                                                                             ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                                                           >  1:1200187/1‑91 (MQ=255)                                                                                                             ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                                                           <  2:1200187/91‑1 (MQ=255)                                                                                                             ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                >  2:162445/1‑134 (MQ=255)                                                                                                             ttttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                <  1:162445/134‑1 (MQ=255)                                                                                                               ttGGGTAGATTTGCTCACATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGGCCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGTGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCga         <  1:794678/139‑2 (MQ=255)                                                                                                               ttGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt         <  2:662080/139‑1 (MQ=255)                                                                                                                tGGGTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa        >  2:356002/1‑139 (MQ=255)                                                                                                                  ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                  <  2:189970/47‑1 (MQ=255)                                                                                                                  ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCa                                                                                                  >  1:189970/1‑47 (MQ=255)                                                                                                                  ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                  >  2:903911/1‑127 (MQ=255)                                                                                                                  ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                  <  1:903911/127‑1 (MQ=255)                                                                                                                  ggTAGATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGt      <  1:381359/139‑1 (MQ=255)                                                                                                                      gATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCata  >  1:593965/1‑139 (MQ=255)                                                                                                                      gATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGATCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCata  >  2:1290898/1‑139 (MQ=255)                                                                                                                       aTTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            >  1:788986/1‑48 (MQ=255)                                                                                                                       aTTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGa                                                                                            <  2:788986/48‑1 (MQ=255)                                                                                                                                 aTTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa                                           >  1:382982/1‑87 (MQ=255)                                                                                                                                 aTTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa                                           <  2:382982/87‑1 (MQ=255)                                                                                                                                       |                                                                                                                          GATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATA  >  NC_000913/4116989‑4117244

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