BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,062	C→A	79.8%	W54L (TGG→TTG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,062	0	C	A	79.8%	256.8 / 82.5	119	W54L (TGG→TTG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (8/16); new base A (48/47); total (56/63)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

Only 100 of 120 total aligned reads displayed.

AATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAT  >  NC_000913/4116917‑4117208

aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGa                                                                                                                                                 <  1:334948/149‑1 (MQ=255)

     cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTcc                                                                                                                                            <  2:788122/149‑1 (MQ=255)

     cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTcc                                                                                                                                            <  3:315630/149‑1 (MQ=255)

      gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCa                                                                                                                                           <  2:81779/149‑1 (MQ=255)

       ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCAt                                                                                                                                          >  3:770309/1‑149 (MQ=255)

         atagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTg                                                                                                                                        <  1:738258/149‑1 (MQ=255)

          tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATAgg                                                                                                                                       <  3:703087/149‑1 (MQ=255)

           agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTggg                                                                                                                                      >  2:322573/1‑149 (MQ=255)

           agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTggg                                                                                                                                      >  3:865561/1‑149 (MQ=255)

            gatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGt                                                                                                                                     >  2:535283/1‑149 (MQ=255)

             atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTc                                                                                                                                    <  3:47035/149‑1 (MQ=255)

             atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCCGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCAAGATTTCCATTGGTTc                                                                                                                                    >  2:417127/1‑149 (MQ=255)

                ggattGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg                                                                                                                                 <  2:40999/146‑1 (MQ=255)

                ggCTTGCCGGTTTC‑TTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCgtgt                                                                                                                                <  2:892405/149‑1 (MQ=255)

                  cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGtt                                                                                                                               <  2:78082/149‑1 (MQ=255)

                   ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTc                                                                                                                              >  1:252055/1‑149 (MQ=255)

                         gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCa                                                                                                                        >  1:430489/1‑149 (MQ=255)

                         gTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCa                                                                                                                        >  4:354488/1‑149 (MQ=255)

                          tttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                                                                                                       >  4:802802/1‑149 (MQ=255)

                           ttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAAc                                                                                                                      >  3:74809/1‑149 (MQ=255)

                              tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                   >  1:641550/1‑149 (MQ=255)

                              tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                   >  3:143425/1‑149 (MQ=255)

                              tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                   >  1:156902/1‑149 (MQ=255)

                              tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                   >  4:837337/1‑149 (MQ=255)

                                 ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                                                                                <  3:145008/149‑1 (MQ=255)

                                   cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                                              <  1:322573/149‑1 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            <  4:114849/149‑1 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            <  1:13175/149‑1 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            >  3:801867/1‑149 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            >  3:283866/1‑149 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            <  4:648559/149‑1 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            >  2:921878/1‑149 (MQ=255)

                                     cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            <  1:359503/149‑1 (MQ=255)

                                      aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                           >  3:102533/1‑149 (MQ=255)

                                        aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                                                                                         >  1:910328/1‑149 (MQ=255)

                                          aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                                                                                                       >  2:221774/1‑149 (MQ=255)

                                          aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                                                                                                       >  4:400340/1‑149 (MQ=255)

                                           aTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                                                                                                      <  4:865561/149‑1 (MQ=255)

                                             aGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                                                                                                    <  4:283866/149‑1 (MQ=255)

                                               tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                                                                                  <  1:19576/149‑1 (MQ=255)

                                                 gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTc                                                                                                >  3:874506/1‑149 (MQ=255)

                                                  tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa                                                                                               >  2:757407/1‑149 (MQ=255)

                                                    tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                                                                             <  2:181449/149‑1 (MQ=255)

                                                        gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                                                                                  >  2:240182/1‑140 (MQ=255)

                                                        gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                                                                                  <  1:240182/140‑1 (MQ=255)

                                                        gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcg                                                                                         >  3:575997/1‑149 (MQ=255)

                                                         cTGGTTCGTAATACCGATAGCTACAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCGCAAATTcgc                                                                                        <  1:803513/149‑1 (MQ=255)

                                                             ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                                                                                    >  2:932597/1‑149 (MQ=255)

                                                             ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                                                                                    >  3:4979/1‑149 (MQ=255)

                                                              tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                   <  2:389134/149‑1 (MQ=255)

                                                              tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                   <  4:634013/149‑1 (MQ=255)

                                                              tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                   <  1:932597/149‑1 (MQ=255)

                                                                gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            <  1:625083/122‑1 (MQ=255)

                                                                gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                            >  2:625083/1‑122 (MQ=255)

                                                                  aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                                                                               <  1:393344/149‑1 (MQ=255)

                                                                      ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCt                                                                           <  2:749911/149‑1 (MQ=255)

                                                                        gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                         <  3:721906/149‑1 (MQ=255)

                                                                        gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                         >  4:43341/1‑149 (MQ=255)

                                                                           aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGa                                                                      <  3:202613/149‑1 (MQ=255)

                                                                              tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat                                                                   >  2:653583/1‑149 (MQ=255)

                                                                                 aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                  <  1:9753/147‑1 (MQ=255)

                                                                                 aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                  >  2:9753/1‑147 (MQ=255)

                                                                                    ttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATAGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCAt                                                             <  4:211104/149‑1 (MQ=255)

                                                                                     tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc                                                            <  1:412140/149‑1 (MQ=255)

                                                                                     tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc                                                            <  4:459877/149‑1 (MQ=255)

                                                                                       aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                          >  3:470444/1‑149 (MQ=255)

                                                                                        tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTg                                                         >  4:177750/1‑149 (MQ=255)

                                                                                         cgcaaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGa                                                        <  2:564739/146‑1 (MQ=255)

                                                                                          ggAACTGATATCTGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTTCATTGTATCGGGTTCTACCCAACCCGGTTTTGGGTAGGTTTGCCCAAAATCGCGCTGCGAGACGCTAATGATATTGGCATCgggct                                                       <  2:562469/149‑5 (MQ=255)

                                                                                           gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTa                                                                                                                            >  2:631441/1‑79 (MQ=255)

                                                                                             aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCa                                                    <  3:465181/149‑1 (MQ=255)

                                                                                              cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAt                                                   <  4:801867/149‑1 (MQ=255)

                                                                                                gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                 <  1:59460/149‑1 (MQ=255)

                                                                                                gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                 <  1:535283/149‑1 (MQ=255)

                                                                                                    tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                                                                               <  1:814706/115‑1 (MQ=255)

                                                                                                    tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                                                                               >  2:814706/1‑115 (MQ=255)

                                                                                                    tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGAc                                             >  4:281624/1‑149 (MQ=255)

                                                                                                     cGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTTGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACTCGCTAATGATATTGGCATCGTGATCCATTACGAcc                                            <  4:149456/149‑1 (MQ=255)

                                                                                                      ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcg                                           >  2:213620/1‑149 (MQ=255)

                                                                                                       gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgc                                          >  3:125674/1‑149 (MQ=255)

                                                                                                              ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                   >  3:705092/1‑149 (MQ=255)

                                                                                                              ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                   >  4:575016/1‑149 (MQ=255)

                                                                                                                  cACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtgg                               <  1:453270/149‑1 (MQ=255)

                                                                                                                  cACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtgg                               <  1:684266/149‑1 (MQ=255)

                                                                                                                   aCTTCTACCAGCGTGGAGCGTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggt                              <  4:364999/149‑1 (MQ=255)

                                                                                                                    cTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                             <  3:739717/149‑1 (MQ=255)

                                                                                                                    cTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                             <  4:4979/149‑1 (MQ=255)

                                                                                                                    cTTCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                             <  2:430489/149‑1 (MQ=255)

                                                                                                                     ttCTACCAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                            >  3:473688/1‑149 (MQ=255)

                                                                                                                          ccAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                       >  1:910860/1‑149 (MQ=255)

                                                                                                                           cAGCGTGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                      >  1:916188/1‑149 (MQ=255)

                                                                                                                                tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                                                                             <  2:421197/57‑1 (MQ=255)

                                                                                                                                tGGAGCTTTGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                 >  3:739678/1‑149 (MQ=255)

                                                                                                                                   aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc              <  4:874506/149‑1 (MQ=255)

                                                                                                                                        tGGGTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCATGGTCGAGCGCAAc         <  4:102533/149‑1 (MQ=255)

                                                                                                                                          ggTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCggg                                     <  4:545125/119‑1 (MQ=255)

                                                                                                                                          ggTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa       <  2:489999/149‑1 (MQ=255)

                                                                                                                                              gCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGTTTTGCTCAAATTCGCGCTGCGACACTCTAATGAGATTGGCATCGTGATCCATTACGACCTCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatata   >  3:90/1‑149 (MQ=255)

                                                                                                                                              gCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatata   >  2:579934/1‑149 (MQ=255)

                                                                                                                                               cccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatatat  >  3:698749/1‑149 (MQ=255)

AATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAT  >  NC_000913/4116917‑4117208

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