BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,141,722	A→G	L153P (CTG→CCG)	hybB ←	putative cytochrome Ni/Fe component of hydrogenase‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,141,722	0	A	G	100.0%	666.4 / NA	206	L153P (CTG→CCG)	hybB	putative cytochrome Ni/Fe component of hydrogenase‑2
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (110/96); total (110/96)

Only 100 of 215 total aligned reads displayed.

CATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGA  >  NC_000913/3141631‑3141811

cATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACACCAGACGTTCAAACAg                                                                                          >  1:5118555/1‑93 (MQ=255)

 aTGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGt                                                                                         >  2:1500500/1‑93 (MQ=255)

  tGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTg                                                                                        <  1:4989599/93‑1 (MQ=255)

    gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGcc                                                                                      >  2:6046453/1‑93 (MQ=255)

       ggCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGa                                                                                   <  2:6271501/93‑1 (MQ=255)

           gcagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaa                                                                               >  1:6076041/1‑93 (MQ=255)

            cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCaaa                                                                              <  2:4844287/93‑1 (MQ=255)

             agcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAAc                                                                             <  2:3756183/93‑1 (MQ=255)

             agcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAAc                                                                             <  2:5268955/93‑1 (MQ=255)

              gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACt                                                                            >  1:3046807/1‑93 (MQ=255)

              gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACt                                                                            >  2:3668763/1‑93 (MQ=255)

               cgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTc                                                                           <  1:1341990/93‑1 (MQ=255)

                gcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTcc                                                                          <  2:4173748/93‑1 (MQ=255)

                 cACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCa                                                                         <  2:4883803/93‑1 (MQ=255)

                 cACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCa                                                                         <  2:5716398/93‑1 (MQ=255)

                 cACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCa                                                                         <  1:819575/93‑1 (MQ=255)

                 cACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCa                                                                         <  2:2046103/93‑1 (MQ=255)

                  aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        <  1:2576112/93‑1 (MQ=255)

                  aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        >  1:4710252/1‑93 (MQ=255)

                  aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAg                                                                        >  2:242061/1‑93 (MQ=255)

                   ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCAGCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAgg                                                                       >  1:66636/1‑92 (MQ=255)

                   ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                       >  2:5517076/1‑93 (MQ=255)

                   ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                       <  2:6208512/93‑1 (MQ=255)

                   ccGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGt                                                                       >  2:2870218/1‑93 (MQ=255)

                    cGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTg                                                                      >  1:4142700/1‑93 (MQ=255)

                    cGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTg                                                                      >  2:1339475/1‑93 (MQ=255)

                      aGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGcc                                                                    >  1:1520021/1‑93 (MQ=255)

                         gcgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATc                                                                 >  1:906986/1‑93 (MQ=255)

                          cgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCa                                                                <  1:156377/93‑1 (MQ=255)

                          cgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCa                                                                <  2:5580586/93‑1 (MQ=255)

                           gatgatGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCAc                                                               <  1:2079913/93‑1 (MQ=255)

                           gatgatAAAGAACATCACCTTGTTAAGTCGCTGTAGCGACGCCTTCCTCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCAc                                                               <  1:3296662/93‑1 (MQ=255)

                            atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACg                                                              >  2:2665150/1‑93 (MQ=255)

                             tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             <  1:546634/93‑1 (MQ=255)

                             tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             >  2:5897777/1‑93 (MQ=255)

                             tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             >  1:1215830/1‑93 (MQ=255)

                             tgatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGc                                                             >  1:2763117/1‑93 (MQ=255)

                              gatgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGcc                                                            >  2:4330601/1‑93 (MQ=255)

                                tgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGa                                                          <  1:4259238/93‑1 (MQ=255)

                                 gaataaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         <  1:494280/93‑1 (MQ=255)

                                 gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         <  2:2642075/93‑1 (MQ=255)

                                 gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         <  1:5659627/93‑1 (MQ=255)

                                 gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         <  2:2211313/93‑1 (MQ=255)

                                 gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                         >  1:6175883/1‑93 (MQ=255)

                                    gaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGAtata                                                      >  1:3762035/1‑93 (MQ=255)

                                      aCATCGCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                    <  1:4849374/93‑1 (MQ=255)

                                      aCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                    <  2:2566752/93‑1 (MQ=255)

                                      aCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                    >  1:1232769/1‑93 (MQ=255)

                                        aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                  <  2:4247286/93‑1 (MQ=255)

                                           aCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTc                                               <  1:3711345/93‑1 (MQ=255)

                                           aCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTc                                               >  1:2438096/1‑93 (MQ=255)

                                            ccTTGTTTAGTCGCTGTAGCGACTCCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCa                                              <  2:1758989/93‑1 (MQ=255)

                                             cTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCAt                                             >  2:2228044/1‑93 (MQ=255)

                                              ttGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTGCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                            <  1:1151128/93‑1 (MQ=255)

                                               tGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATAc                                           <  1:5947240/93‑1 (MQ=255)

                                                gTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACa                                          >  1:4625377/1‑93 (MQ=255)

                                                gTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACa                                          <  2:610420/93‑1 (MQ=255)

                                                 tttAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAg                                         >  1:2068957/1‑93 (MQ=255)

                                                   tAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                       <  2:2814415/93‑1 (MQ=255)

                                                   tAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                       >  1:5529021/1‑93 (MQ=255)

                                                    aGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAcc                                      <  1:3483350/93‑1 (MQ=255)

                                                    aGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAcc                                      <  2:540777/93‑1 (MQ=255)

                                                     gTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccg                                     >  1:1711668/1‑93 (MQ=255)

                                                     gTCGCTGTACAGACCACTTCCACCCCAGACGTTCAAACGGTGCCGGAGCACACTCCAGTGCCATCACGCCGATATAGATGGTCATACAAAcac                                     >  2:1891234/1‑91 (MQ=255)

                                                      tCGCTGTAGCGACCCCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgc                                    <  2:748590/93‑1 (MQ=255)

                                                      tCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgc                                    >  2:1612299/1‑93 (MQ=255)

                                                        gCTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgccg                                  >  1:5435218/1‑93 (MQ=255)

                                                         cTGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGt                                 >  1:3669049/1‑93 (MQ=255)

                                                          tGTAGCGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtc                                >  1:1861711/1‑93 (MQ=255)

                                                               cGACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGaa                           >  2:2113173/1‑93 (MQ=255)

                                                                gACACCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAAc                          >  1:4366717/1‑93 (MQ=255)

                                                                  cacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAATAg                        <  2:3476260/93‑1 (MQ=255)

                                                                  cacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAg                        >  1:4643601/1‑93 (MQ=255)

                                                                  cacCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAg                        >  1:2687573/1‑93 (MQ=255)

                                                                   acCTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGt                       <  2:567165/93‑1 (MQ=255)

                                                                   acCTCCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGt                       <  2:5511338/93‑1 (MQ=255)

                                                                     cTTCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTAc                     <  1:1407464/93‑1 (MQ=255)

                                                                      ttCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTAcc                    >  2:309242/1‑93 (MQ=255)

                                                                       tcctccCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGTCCGCCGTCTCGAACAGTACCg                   <  2:5118555/93‑1 (MQ=255)

                                                                       tCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCg                   >  2:5042555/1‑93 (MQ=255)

                                                                       tCCACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCg                   <  2:5822876/93‑1 (MQ=255)

                                                                         cACCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAg                 <  2:708726/93‑1 (MQ=255)

                                                                          cccccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                <  1:3908003/92‑1 (MQ=255)

                                                                          aCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                <  2:4284426/93‑1 (MQ=255)

                                                                          aCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                >  1:1663319/1‑93 (MQ=255)

                                                                          aCCCCAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                >  1:4898722/1‑93 (MQ=255)

                                                                           ccccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGtt               >  2:1944384/1‑93 (MQ=255)

                                                                            cccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTc              >  2:3179579/1‑93 (MQ=255)

                                                                            cccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTc              >  2:2471072/1‑93 (MQ=255)

                                                                             ccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCa             <  2:3915934/93‑1 (MQ=255)

                                                                             ccAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCa             <  1:320476/93‑1 (MQ=255)

                                                                              cAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCAc            >  1:3540649/1‑93 (MQ=255)

                                                                              cAGACGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCAc            >  1:5569621/1‑93 (MQ=255)

                                                                                 aCGTTCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGtt         <  1:2252955/93‑1 (MQ=255)

                                                                                   ggccAAACGGTGCGGGAGCAATCTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGa       <  2:5309348/90‑1 (MQ=255)

                                                                                   gTTCAAACGGTGCCGGAGCAAGCTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGa       <  1:1959673/93‑1 (MQ=255)

                                                                                    ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      <  1:5695748/93‑1 (MQ=255)

                                                                                    ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      >  1:3836067/1‑93 (MQ=255)

                                                                                    ttCAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa      >  2:3361967/1‑93 (MQ=255)

                                                                                      cAAACGGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGt    <  1:497277/93‑1 (MQ=255)

CATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGA  >  NC_000913/3141631‑3141811

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