BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,888,375	A→G	L516P (CTT→CCT)	cysJ ←	sulfite reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,888,375	0	A	G	100.0%	876.5 / NA	250	L516P (CTT→CCT)	cysJ	sulfite reductase (NADPH)
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (127/123); total (127/123)

Only 100 of 260 total aligned reads displayed.

CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464

ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGa                                                                                         <  2:6005968/93‑1 (MQ=255)

 catcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGAt                                                                                        >  1:5503016/1‑93 (MQ=255)

 catcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGAt                                                                                        <  1:4125129/93‑1 (MQ=255)

 catcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGAt                                                                                        >  1:3219640/1‑93 (MQ=255)

 catcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGAt                                                                                        >  2:3109010/1‑93 (MQ=255)

  atcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATc                                                                                       <  1:3344881/93‑1 (MQ=255)

   tcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCg                                                                                      <  2:6261458/93‑1 (MQ=255)

   tcatTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCg                                                                                      <  2:1880946/93‑1 (MQ=255)

      ttGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATa                                                                                   >  2:5014237/1‑93 (MQ=255)

      ttGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATa                                                                                   <  1:6228316/93‑1 (MQ=255)

      ttGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATa                                                                                   >  2:2746886/1‑93 (MQ=255)

         aTCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgt                                                                                >  2:775166/1‑93 (MQ=255)

         aTCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgt                                                                                >  1:5077076/1‑93 (MQ=255)

           ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtgt                                                                              >  1:5589790/1‑93 (MQ=255)

           ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtgt                                                                              >  1:342892/1‑93 (MQ=255)

            cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTc                                                                             <  2:5114838/93‑1 (MQ=255)

            cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTc                                                                             >  1:3470462/1‑93 (MQ=255)

             aGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCa                                                                            <  2:1473076/93‑1 (MQ=255)

              gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAg                                                                           >  1:239657/1‑93 (MQ=255)

                 ccACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCAc                                                                        >  2:3461202/1‑93 (MQ=255)

                   acaGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGc                                                                      >  2:13795/1‑93 (MQ=255)

                   acaGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGc                                                                      <  1:1849074/93‑1 (MQ=255)

                    caGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGGAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGGAAGcc                                                                     >  1:2607906/1‑93 (MQ=255)

                     aGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCa                                                                    >  1:3721049/1‑93 (MQ=255)

                     aGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCa                                                                    >  1:2318408/1‑93 (MQ=255)

                      gCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCAt                                                                   <  1:1990502/93‑1 (MQ=255)

                      gCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCAt                                                                   <  1:6343743/93‑1 (MQ=255)

                         ccGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAATCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCtt                                                                <  2:3404447/93‑1 (MQ=255)

                         ccGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCtt                                                                >  1:5405657/1‑93 (MQ=255)

                         ccGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCtt                                                                <  2:2220831/93‑1 (MQ=255)

                          cgcgCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCttt                                                               <  1:783628/93‑1 (MQ=255)

                             gcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGAc                                                            >  2:1530078/1‑93 (MQ=255)

                             gcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGAc                                                            <  1:6445347/93‑1 (MQ=255)

                               ccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACga                                                          >  2:2120003/1‑92 (MQ=255)

                               ccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGt                                                          <  2:3625235/93‑1 (MQ=255)

                               ccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGt                                                          <  2:4714205/93‑1 (MQ=255)

                               ccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGt                                                          <  2:6206089/93‑1 (MQ=255)

                                 tGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAg                                                        >  1:4523746/1‑93 (MQ=255)

                                 tGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAg                                                        <  1:133062/93‑1 (MQ=255)

                                  gTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgc                                                       <  1:6377223/93‑1 (MQ=255)

                                  gTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgc                                                       <  2:2455/93‑1 (MQ=255)

                                    tCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgcgc                                                     <  2:2591074/93‑1 (MQ=255)

                                     cgcgCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCt                                                    <  1:4288749/93‑1 (MQ=255)

                                      gcgcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTg                                                   >  2:4158629/1‑93 (MQ=255)

                                       cgcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGc                                                  >  1:235638/1‑93 (MQ=255)

                                       cgcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGc                                                  >  1:144691/1‑93 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  2:3984024/1‑93 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  2:5244853/1‑93 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  2:113928/1‑93 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               >  1:1147906/1‑93 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               <  2:5726793/93‑1 (MQ=255)

                                          aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCAc                                               <  1:5885504/93‑1 (MQ=255)

                                           gTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACt                                              >  1:2106547/1‑93 (MQ=255)

                                           gTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACt                                              <  1:406132/93‑1 (MQ=255)

                                              tGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCa                                           <  2:3143196/93‑1 (MQ=255)

                                                 cTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCt                                        <  1:1573967/93‑1 (MQ=255)

                                                 cTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCt                                        >  1:2475780/1‑93 (MQ=255)

                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       <  1:2970869/93‑1 (MQ=255)

                                                       cGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACa                                  >  2:3228250/1‑93 (MQ=255)

                                                         tAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAgg                                >  1:2838053/1‑93 (MQ=255)

                                                          aaaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGa                               >  1:2229753/1‑93 (MQ=255)

                                                          aaaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGa                               >  2:4572572/1‑93 (MQ=255)

                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              <  2:5589277/93‑1 (MQ=255)

                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              >  2:2961480/1‑93 (MQ=255)

                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              <  2:4267494/93‑1 (MQ=255)

                                                           aaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGaa                              <  1:2510707/93‑1 (MQ=255)

                                                             ctttttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGt                            >  2:1967016/1‑93 (MQ=255)

                                                               ttttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCt                          <  1:4428024/93‑1 (MQ=255)

                                                                  tctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTcc                       <  2:6316239/93‑1 (MQ=255)

                                                                  tctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTcc                       <  1:415163/93‑1 (MQ=255)

                                                                    tttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGt                     >  1:5160151/1‑93 (MQ=255)

                                                                     ttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTa                    <  1:2727317/93‑1 (MQ=255)

                                                                      tttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaa                   <  2:446772/93‑1 (MQ=255)

                                                                       ttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaaa                  <  1:1021108/93‑1 (MQ=255)

                                                                       ttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaaa                  <  1:3804764/93‑1 (MQ=255)

                                                                       ttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaaa                  >  1:498286/1‑93 (MQ=255)

                                                                       ttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTaaa                  <  1:4956052/93‑1 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAgg                >  2:2289519/1‑92 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGt                <  2:5424934/93‑1 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGt                >  2:2306546/1‑93 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGt                <  1:1768188/93‑1 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGt                >  1:4193149/1‑93 (MQ=255)

                                                                         gATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGGCTTCCGTAAAgg                >  1:6139472/1‑92 (MQ=255)

                                                                          aTCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTg               >  2:2374583/1‑93 (MQ=255)

                                                                          aTCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTg               >  2:1765260/1‑93 (MQ=255)

                                                                           tCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGc              >  1:4789554/1‑93 (MQ=255)

                                                                            cgcgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCg             <  2:1332339/93‑1 (MQ=255)

                                                                             gcgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCgg            >  2:3390458/1‑93 (MQ=255)

                                                                              cgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGa           <  2:333021/93‑1 (MQ=255)

                                                                                cgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGAtt         >  2:1638577/1‑93 (MQ=255)

                                                                                cgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGAtt         <  1:5120802/93‑1 (MQ=255)

                                                                                 gACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTCCAGGAAGTCTTCCGTAAAGTGCGGATTa        >  2:1054222/1‑93 (MQ=255)

                                                                                 gACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTa        >  1:4712207/1‑93 (MQ=255)

                                                                                  aCCAGGCAGGATCGATACGTGTCAGCACGCCCTCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTAc       >  1:531395/1‑93 (MQ=255)

                                                                                  aCCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTAc       <  2:4748175/93‑1 (MQ=255)

                                                                                    caggcaggATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCa     >  2:4098353/1‑93 (MQ=255)

                                                                                     aGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaa    >  1:1536219/1‑93 (MQ=255)

                                                                                     aGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaa    <  1:664166/93‑1 (MQ=255)

                                                                                      ggCAGGATAGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaac   >  1:6290997/1‑92 (MQ=255)

                                                                                       gCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAg  >  2:2629319/1‑93 (MQ=255)

CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464

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