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%	1159.7 / NA	335	L516P (CTT→CCT)	cysJ	sulfite reductase (NADPH)
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (172/163); total (172/163)

Only 100 of 352 total aligned reads displayed.

CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464

ccATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGa                                                                                         >  1:7864286/1‑93 (MQ=255)

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

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

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

    catTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGa                                                                                     <  1:5701032/93‑1 (MQ=255)

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

       tGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATAc                                                                                  <  1:8336222/93‑1 (MQ=255)

       tGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATAc                                                                                  >  1:2448973/1‑93 (MQ=255)

        gATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACg                                                                                 >  2:6301409/1‑93 (MQ=255)

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

          tCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtt                                                                               >  2:7241246/1‑92 (MQ=255)

           ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACgtgt                                                                              <  1:8551252/93‑1 (MQ=255)

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

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

             aGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGGCa                                                                            >  1:7986694/1‑93 (MQ=255)

              gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAg                                                                           <  2:7631998/93‑1 (MQ=255)

              gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAg                                                                           <  2:3395194/93‑1 (MQ=255)

                 ccACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATGGCGCGACCAGGCAGGATCGATACGTGTCAGCAc                                                                        >  1:6095934/1‑93 (MQ=255)

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

                   acaGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGAGCGACCAGGCAGGATCGATACGTGTCAGCACGc                                                                      >  2:7575032/1‑93 (MQ=255)

                    caGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATAACAGTCAGCACGcc                                                                     >  1:2221226/1‑93 (MQ=255)

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

                      gCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCAt                                                                   >  1:5057170/1‑93 (MQ=255)

                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 <  2:2238108/93‑1 (MQ=255)

                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 >  1:7617032/1‑93 (MQ=255)

                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 <  1:2820469/93‑1 (MQ=255)

                        tCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCt                                                                 <  2:7744082/93‑1 (MQ=255)

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

                          cgcgcgCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCttt                                                               <  2:5040114/93‑1 (MQ=255)

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

                           gcgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTg                                                              <  1:6211883/93‑1 (MQ=255)

                           gcgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTg                                                              <  1:7472363/93‑1 (MQ=255)

                           gcgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCCACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTg                                                              <  2:1731758/93‑1 (MQ=255)

                            cgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGa                                                             <  1:7008169/93‑1 (MQ=255)

                              cccTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACg                                                           >  1:4166242/1‑93 (MQ=255)

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

                                cTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTa                                                         <  1:2712765/93‑1 (MQ=255)

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

                                   ttCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAgcg                                                      >  2:300591/1‑93 (MQ=255)

                                    tCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCCCGCCATCTTTGACGTAgcgc                                                     <  2:7179074/93‑1 (MQ=255)

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

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

                                     cgcgCAGTTTGTCTTGTTCGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCt                                                    <  2:2596333/93‑1 (MQ=255)

                                        gcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGcc                                                 >  2:200902/1‑93 (MQ=255)

                                        gcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGcc                                                 >  2:5355362/1‑93 (MQ=255)

                                        gcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGcc                                                 <  2:6880203/93‑1 (MQ=255)

                                        gcAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGcc                                                 >  2:1824218/1‑93 (MQ=255)

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

                                           gTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATAAGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACt                                              >  1:4009250/1‑93 (MQ=255)

                                            tttgttttgtACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTc                                             <  2:8010445/93‑1 (MQ=255)

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

                                              tGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCa                                           >  2:1517341/1‑93 (MQ=255)

                                                tCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCAcc                                         <  1:8768347/93‑1 (MQ=255)

                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCCCCTg                                       >  1:3117632/1‑93 (MQ=255)

                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       >  2:4918615/1‑93 (MQ=255)

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

                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTg                                       >  1:1961891/1‑93 (MQ=255)

                                                  ttGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTCACGTAGCGCTGCCACTCCACCTg                                       >  2:3286143/1‑93 (MQ=255)

                                                   tGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTgg                                      >  2:3959040/1‑93 (MQ=255)

                                                   tGTACGTAAACTTTTTCTTTTTGATCGCGCGAACAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTgg                                      >  2:216249/1‑93 (MQ=255)

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

                                                        gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAg                                 >  1:3499738/1‑93 (MQ=255)

                                                        gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAg                                 >  1:3365694/1‑93 (MQ=255)

                                                        gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAg                                 <  1:6442295/93‑1 (MQ=255)

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

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

                                                          aaaCTTTTTCTTTTTGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGa                               <  1:5306141/93‑1 (MQ=255)

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

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

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

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

                                                               ttttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCt                          >  1:1616139/1‑93 (MQ=255)

                                                                tttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCtt                         >  1:2505077/1‑93 (MQ=255)

                                                                 ttctttttGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTc                        <  1:7700252/93‑1 (MQ=255)

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

                                                                   ctttttGATCGCGCGCCCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCg                      <  2:4247654/93‑1 (MQ=255)

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

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

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

                                                                        tGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAg                 >  1:5930925/1‑93 (MQ=255)

                                                                        tGATCGCGCGACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAg                 >  1:1151290/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                            cgcgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGGGCg             >  1:6748849/1‑93 (MQ=255)

                                                                              cgcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGa           >  2:3605180/1‑93 (MQ=255)

                                                                               gcgACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGAt          >  2:8444658/1‑93 (MQ=255)

                                                                                 gACCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTa        <  2:380417/93‑1 (MQ=255)

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

                                                                                  aCCAGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTAc       >  1:4490126/1‑93 (MQ=255)

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

                                                                                    caggcaggATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACTGGAAGTCTTCCGTAAAGTGCGGATTACCa     >  2:7954813/1‑93 (MQ=255)

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

                                                                                     aGGCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCCGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCTGATTACCaa    >  1:4792139/1‑93 (MQ=255)

                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   <  1:8789511/93‑1 (MQ=255)

                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   >  1:2159490/1‑93 (MQ=255)

                                                                                      ggCAGGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCaaa   >  1:1265581/1‑93 (MQ=255)

CCATCATTGATCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAG  >  NC_000913/2888285‑2888464

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