BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,777,301	C→T	L63L (CTG→TTG)	yfjZ →	CP4‑57 prophage; putative antitoxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,777,301	0	C	T	100.0%	1342.8 / NA	312	L63L (CTG→TTG)	yfjZ	CP4‑57 prophage; putative antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (194/118); total (194/118)

Only 100 of 312 total aligned reads displayed.

TGGCTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTA  >  NC_000913/2777203‑2777398

tGGCTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTg                                                                                                 <  1:5049350/101‑1 (MQ=255)

 ggCTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGa                                                                                                <  1:5049411/101‑1 (MQ=255)

 ggCTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGa                                                                                                <  1:5049618/101‑1 (MQ=255)

   cTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGAcc                                                                                              >  1:5049945/1‑101 (MQ=255)

      aCCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACt                                                                                           <  2:5049281/101‑1 (MQ=255)

       ccGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTg                                                                                          <  2:5049539/101‑1 (MQ=255)

         gggCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGAc                                                                                               >  2:5049951/1‑94 (MQ=255)

         gggCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGt                                                                                        >  1:5049950/1‑101 (MQ=255)

            ccAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGaa                                                                                     <  2:5049231/101‑1 (MQ=255)

            ccAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGaa                                                                                     <  2:5049533/101‑1 (MQ=255)

             cAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAAc                                                                                    <  1:5049061/101‑1 (MQ=255)

             cAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAAc                                                                                    >  2:5049954/1‑101 (MQ=255)

              aGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACt                                                                                   <  1:5049303/101‑1 (MQ=255)

                tATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGa                                                                                 <  1:5049518/101‑1 (MQ=255)

                 aTCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGaa                                                                                >  1:5049958/1‑101 (MQ=255)

                   cACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATc                                                                              <  1:5049160/101‑1 (MQ=255)

                   cACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATc                                                                              <  1:5049204/101‑1 (MQ=255)

                    aCCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATcc                                                                             <  1:5049510/101‑1 (MQ=255)

                    aCCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATcc                                                                             >  2:5049961/1‑101 (MQ=255)

                     ccGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATccc                                                                            <  1:5049609/101‑1 (MQ=255)

                      cGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATcccc                                                                           >  1:5049964/1‑101 (MQ=255)

                        gTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGc                                                                         <  2:5049638/101‑1 (MQ=255)

                         tAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAAt                                                                               <  1:5049454/94‑1 (MQ=255)

                         tAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGcc                                                                        <  1:5049573/101‑1 (MQ=255)

                         tAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGcc                                                                        >  2:5049968/1‑101 (MQ=255)

                          aaGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCa                                                                       >  1:5049970/1‑101 (MQ=255)

                           aGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCAt                                                                      <  2:5049432/101‑1 (MQ=255)

                           aGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCAt                                                                      >  2:5049974/1‑101 (MQ=255)

                            gTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATg                                                                     >  1:5049976/1‑101 (MQ=255)

                            gTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATg                                                                     >  1:5049979/1‑101 (MQ=255)

                            gTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATg                                                                     <  1:5049447/101‑1 (MQ=255)

                              ttAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGcc                                                                   >  2:5049983/1‑101 (MQ=255)

                               tAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGccc                                                                  >  1:5049987/1‑101 (MQ=255)

                                aGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCa                                                                 <  1:5049779/101‑1 (MQ=255)

                                 gTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCaa                                                                >  2:5049991/1‑101 (MQ=255)

                                  tGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAAt                                                               >  1:5049992/1‑101 (MQ=255)

                                   gACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATg                                                              <  1:5049293/101‑1 (MQ=255)

                                    aCGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGc                                                             <  1:5049470/101‑1 (MQ=255)

                                    aCGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGc                                                             <  2:5049302/101‑1 (MQ=255)

                                      gCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGt                                                           >  1:5049996/1‑101 (MQ=255)

                                       ccGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTc                                                          >  2:5050001/1‑101 (MQ=255)

                                       ccGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTc                                                          >  1:5050002/1‑101 (MQ=255)

                                        cGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCa                                                         >  2:5050006/1‑101 (MQ=255)

                                         gAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCAc                                                        >  1:5050009/1‑101 (MQ=255)

                                          aaTGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCAcc                                                       >  1:5050012/1‑101 (MQ=255)

                                          aaTGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCAcc                                                       >  2:5050014/1‑101 (MQ=255)

                                            tGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCt                                                     <  1:5049555/101‑1 (MQ=255)

                                            tGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCt                                                     >  2:5050018/1‑101 (MQ=255)

                                            tGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCt                                                     <  1:5049318/101‑1 (MQ=255)

                                             gTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTg                                                    <  2:5049571/101‑1 (MQ=255)

                                              tCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGt                                                   >  2:5050022/1‑101 (MQ=255)

                                               ccTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTa                                                  >  1:5050023/1‑101 (MQ=255)

                                                 tAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTAcc                                                >  1:5050026/1‑101 (MQ=255)

                                                  aaGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTAc                                                 >  2:5050031/1‑99 (MQ=255)

                                                    gCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACcaca                                             >  1:5050035/1‑101 (MQ=255)

                                                     cTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACaa                                            >  2:5050036/1‑101 (MQ=255)

                                                     cTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACaa                                            >  1:5050039/1‑101 (MQ=255)

                                                      tGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAAc                                           <  2:5049697/101‑1 (MQ=255)

                                                      tGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAAc                                           >  1:5050043/1‑101 (MQ=255)

                                                      tGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAAc                                           >  2:5050045/1‑101 (MQ=255)

                                                       ggATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACg                                          <  2:5049132/101‑1 (MQ=255)

                                                       ggATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACg                                          >  1:5050050/1‑101 (MQ=255)

                                                        gATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACgg                                         >  1:5050055/1‑101 (MQ=255)

                                                        gATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACgg                                         <  2:5049394/101‑1 (MQ=255)

                                                         aTGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGt                                        >  2:5050057/1‑101 (MQ=255)

                                                           gtgGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAAc                                           <  2:5049883/96‑1 (MQ=255)

                                                           gtgGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGttt                                      >  2:5050061/1‑101 (MQ=255)

                                                            tgGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGtttt                                     >  2:5050064/1‑101 (MQ=255)

                                                             ggTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTa                                    <  1:5049489/101‑1 (MQ=255)

                                                             ggTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTa                                    <  1:5049420/101‑1 (MQ=255)

                                                              gTATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTAc                                   >  2:5050067/1‑101 (MQ=255)

                                                               tATTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTAcc                                  >  2:5050070/1‑101 (MQ=255)

                                                                aTTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                 <  1:5049581/101‑1 (MQ=255)

                                                                aTTTCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                 >  1:5050073/1‑101 (MQ=255)

                                                                 tttCCACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTg                                >  2:5050074/1‑101 (MQ=255)

                                                                    ccACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGa                             >  2:5050080/1‑101 (MQ=255)

                                                                    ccACATTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGa                             >  1:5050079/1‑101 (MQ=255)

                                                                       caTTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGc                          <  2:5049599/101‑1 (MQ=255)

                                                                       caTTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGc                          >  2:5050084/1‑101 (MQ=255)

                                                                        aTTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCt                                 <  2:5049497/93‑1 (MQ=255)

                                                                        aTTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGcc                         >  1:5050088/1‑101 (MQ=255)

                                                                        aTTTTATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGcc                         <  1:5049680/101‑1 (MQ=255)

                                                                          tttATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGa                       >  2:5050091/1‑101 (MQ=255)

                                                                          tttATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGa                       <  1:5049442/101‑1 (MQ=255)

                                                                            tATCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATa                     >  1:5050094/1‑101 (MQ=255)

                                                                              tCAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATAct                   >  1:5050098/1‑101 (MQ=255)

                                                                               cAGCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATActc                  >  2:5050100/1‑101 (MQ=255)

                                                                                 gCCAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCg                        >  2:5050103/1‑93 (MQ=255)

                                                                                   cAGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTgg              >  1:5050106/1‑101 (MQ=255)

                                                                                    aGATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGt             >  2:5050110/1‑101 (MQ=255)

                                                                                     gATAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTa            >  1:5050112/1‑101 (MQ=255)

                                                                                       tAGAGTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGt          >  2:5050116/1‑101 (MQ=255)

                                                                                        agagTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGtt         >  1:5050118/1‑101 (MQ=255)

                                                                                         gagTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTg        <  1:5049881/101‑1 (MQ=255)

                                                                                         gagTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTg        >  2:5050122/1‑101 (MQ=255)

                                                                                          agTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGc       >  1:5050125/1‑101 (MQ=255)

                                                                                           gTCGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCg      >  2:5050128/1‑101 (MQ=255)

                                                                                             cGATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGc    >  1:5050131/1‑101 (MQ=255)

                                                                                              gATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTa                                                  <  2:5049477/54‑1 (MQ=255)

                                                                                              gATGTTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCt   >  2:5050133/1‑101 (MQ=255)

TGGCTGACCGGGCCAGTATCACCGGTAAGTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTA  >  NC_000913/2777203‑2777398

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