BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,505	G→A	R462R (CGG→CGA)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,505	0	G	A	99.2%	815.4 / ‑5.6	236	R462R (CGG→CGA)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (122/111); total (123/113)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 247 total aligned reads displayed.

TGCCTGGCAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATGC  >  NC_000913/2208415‑2208594

tGCCTGGCAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGAcc                                                                                         >  1:3888832/1‑93 (MQ=255)

   cTGGCAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACctct                                                                                      >  2:5683688/1‑93 (MQ=255)

   cTGGCAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACcgct                                                                                      <  2:8116953/93‑1 (MQ=255)

     gtcAAATGTTGA‑CGCGCCCGTTAAGCGCGCCGGCGATTCTGGCTTTGCGGGAATGTTTGCCACAGAGTCCAACCGCCTGCATACGCCCTCTTa                                                                                    <  2:5408395/91‑1 (MQ=255)

       cAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAAt                                                                                  >  1:6027888/1‑93 (MQ=255)

        aaaTGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCTGCGTCCCCCCTCCTTCCTACCCCCTCTCAATc                                                                                 >  1:9692406/1‑93 (MQ=255)

          aTGTTGAGCGCG‑CCGTTACGCCATCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAg                                                                               <  2:7219638/93‑1 (MQ=255)

           tGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAgg                                                                              >  2:7435673/1‑93 (MQ=255)

            gTTGAGCGCG‑CGGTTACGCCAGCCGGGCATTGTGGCTTGGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             <  1:4604833/93‑1 (MQ=255)

            gTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             <  2:2365436/93‑1 (MQ=255)

            gTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGt                                                                             >  1:5535473/1‑93 (MQ=255)

             ttGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTc                                                                            <  2:9190735/93‑1 (MQ=255)

             ttGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTc                                                                            >  2:3755507/1‑93 (MQ=255)

               gAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGa                                                                          >  1:7573350/1‑93 (MQ=255)

                   gcg‑cCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAAc                                                                      >  1:6844657/1‑93 (MQ=255)

                    cg‑cCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACt                                                                     >  2:2213208/1‑93 (MQ=255)

                    cg‑cCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACt                                                                     >  1:162046/1‑93 (MQ=255)

                    cg‑cCGTTACGCCAGCCGGGCATTGTGGCTTGGCGTTATTATTTACGCCGGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACt                                                                     <  2:5451059/93‑1 (MQ=255)

                     g‑cCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtt                                                                    <  1:3805043/93‑1 (MQ=255)

                     g‑cCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtt                                                                    <  1:3856564/93‑1 (MQ=255)

                       ccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttt                                                                   >  2:8620936/1‑92 (MQ=255)

                       ccGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACtta                                                                   >  1:9916717/1‑93 (MQ=255)

                         gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 >  2:6501945/1‑93 (MQ=255)

                         gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 <  1:934382/93‑1 (MQ=255)

                         gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatt                                                                 <  2:2385077/93‑1 (MQ=255)

                         gTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCCTAATCAGGTCGATAACttatt                                                                 <  1:3094386/93‑1 (MQ=255)

                          ttACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattt                                                                >  1:3486075/1‑93 (MQ=255)

                          ttACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattt                                                                <  1:8777470/93‑1 (MQ=255)

                          ttACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattt                                                                >  1:7732856/1‑93 (MQ=255)

                          ttACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattt                                                                <  2:8780187/93‑1 (MQ=255)

                           tACGCCAGCCGGGCATTGTGGCTTTGCGTGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattta                                                               <  2:957062/93‑1 (MQ=255)

                           tACGCCAGCCGGGCATTGTGGCTTTGCGGGAATCTTTACGCCAGCGTCCCCCCGCCTGCATACGCCCTCTTAATCAGGTCGATAACttattta                                                               >  1:7282819/1‑93 (MQ=255)

                           tACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttattta                                                               <  2:5226110/93‑1 (MQ=255)

                            aCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACttatttat                                                              <  2:6015247/93‑1 (MQ=255)

                            aCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGAGAACttatttat                                                              >  1:6642088/1‑93 (MQ=255)

                            aCGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGAATACGACCTCTTAATCAGGTCGATAACttatttat                                                              <  1:1791302/93‑1 (MQ=255)

                             cGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTtatttatt                                                             <  1:6525033/93‑1 (MQ=255)

                              gCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTc                                                            <  1:4751615/93‑1 (MQ=255)

                              gCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTc                                                            <  2:2738901/93‑1 (MQ=255)

                              gCCAGCCGGGCATTGTGGCTGTGCGGGAATATTTACTCCAGCGACCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTc                                                            <  2:6780667/93‑1 (MQ=255)

                               ccAGCCGGGCATTGTGGCTTGGCGGGAATATTTACGCCAGCGCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCt                                                           <  2:7544260/93‑1 (MQ=255)

                                 aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                         >  2:4545353/1‑93 (MQ=255)

                                 aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                         >  2:9679595/1‑93 (MQ=255)

                                 aGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGc                                                         <  1:4401546/93‑1 (MQ=255)

                                  gCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGcc                                                        >  1:6270866/1‑93 (MQ=255)

                                    cGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg                                                      <  1:3760059/93‑1 (MQ=255)

                                       gCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATTCGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCg                                                   <  2:2850370/93‑1 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCTGGTCGCTg                                                >  1:9450740/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:7127992/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:9404593/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:5706215/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:5229027/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:1460450/1‑93 (MQ=255)

                                          ttGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTg                                                >  1:302117/1‑93 (MQ=255)

                                           tgtgGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGa                                               <  2:1121968/93‑1 (MQ=255)

                                            gtgGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAg                                              >  2:2600693/1‑93 (MQ=255)

                                                 tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                         <  2:4687307/93‑1 (MQ=255)

                                                 tttGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTAt                                         >  2:3079848/1‑93 (MQ=255)

                                                    gCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTAGCTGAGTGTATTTc                                      >  1:4423519/1‑93 (MQ=255)

                                                      gggAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGc                                    >  1:109678/1‑93 (MQ=255)

                                                        gAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTc                                  <  2:1516499/93‑1 (MQ=255)

                                                        gAATATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTc                                  >  1:5381783/1‑93 (MQ=255)

                                                         aaTATTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCg                                 >  1:2176958/1‑93 (MQ=255)

                                                           tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               >  2:729598/1‑93 (MQ=255)

                                                           tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               >  1:2662223/1‑93 (MQ=255)

                                                           tatTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGt                               <  2:5668922/93‑1 (MQ=255)

                                                            atTTACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGtt                              <  2:4189397/93‑1 (MQ=255)

                                                            atTTACGCCAGCGGCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGtt                              <  2:1792507/93‑1 (MQ=255)

                                                             tttACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTg                             <  2:141204/93‑1 (MQ=255)

                                                             tttACGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTg                             >  1:1323094/1‑93 (MQ=255)

                                                              ttccGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTgg                            <  2:2149835/90‑1 (MQ=255)

                                                               tACGCCAGCGTCCCCCCGCCTGCCTACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTggg                           <  1:8953760/93‑1 (MQ=255)

                                                               tACGCCAGCGCCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTggg                           <  2:1066914/93‑1 (MQ=255)

                                                                 cGCCAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAc                         <  1:7987918/93‑1 (MQ=255)

                                                                   ccAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcag                       >  2:8876248/1‑93 (MQ=255)

                                                                    cAGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagc                      <  1:7275358/93‑1 (MQ=255)

                                                                     aGCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagca                     >  1:5912774/1‑93 (MQ=255)

                                                                      gCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagcag                    >  1:2452391/1‑93 (MQ=255)

                                                                      gCGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGAcagcag                    >  2:8528653/1‑93 (MQ=255)

                                                                       cGTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagc                   <  2:4736519/93‑1 (MQ=255)

                                                                       cGTCCCCCCGCCTGCAGACAACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACagcagc                   <  2:6234541/93‑1 (MQ=255)

                                                                        gTCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgcc                  >  1:6737168/1‑93 (MQ=255)

                                                                         tCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccg                 >  2:7515252/1‑93 (MQ=255)

                                                                         tCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccg                 >  2:792983/1‑93 (MQ=255)

                                                                         tCCCCCCGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccg                 >  2:6654429/1‑93 (MQ=255)

                                                                          ccccccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAgccgc                >  2:5858657/1‑93 (MQ=255)

                                                                             cccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCaa             >  1:4157430/1‑92 (MQ=255)

                                                                              ccGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGa            <  2:686622/93‑1 (MQ=255)

                                                                               cGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGAc           >  2:8700272/1‑93 (MQ=255)

                                                                               cGCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGAc           >  1:101163/1‑93 (MQ=255)

                                                                                gCCTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGGATTTCGCTCGGTTGGGACAGCAGCCACCAGACg          >  2:9113082/1‑93 (MQ=255)

                                                                                  cTGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCt        >  2:6955988/1‑93 (MQ=255)

                                                                                   tGCATACTACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTg       <  1:6390031/93‑1 (MQ=255)

                                                                                   tGCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTg       >  1:1974994/1‑93 (MQ=255)

                                                                                    gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg      <  2:7826349/93‑1 (MQ=255)

                                                                                    gCATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTgg      <  1:3366256/93‑1 (MQ=255)

                                                                                     cATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGa     <  1:2252006/93‑1 (MQ=255)

                                                                                     cATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGa     <  1:8366109/93‑1 (MQ=255)

                                                                                     cATACGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCAGCGAGACGCTGGa     >  2:5071895/1‑93 (MQ=255)

                                                                                        aCGACCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATgc  <  2:1961405/93‑1 (MQ=255)

TGCCTGGCAAATGTTGAGCGCG‑CCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGGTCGCTGAGTGTATTTCGCTCGGTTGGGACAGCAGCCGCCAGACGCTGGATGC  >  NC_000913/2208415‑2208594

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