BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,586,406	T→C	Y119C (TAC→TGC)	hsdR ←	type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,586,406	0	T	C	100.0%	1185.0 / NA	282	Y119C (TAC→TGC)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (141/141); total (141/141)

Only 100 of 295 total aligned reads displayed.

GTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCG  >  NC_000913/4586308‑4586503

gttTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcag                                                                                                 >  1:7116032/1‑101 (MQ=255)

gttTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcag                                                                                                 <  2:7115489/101‑1 (MQ=255)

 ttttAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcagt                                                                                                <  2:7115497/101‑1 (MQ=255)

 ttttAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcagt                                                                                                <  2:7115507/101‑1 (MQ=255)

 ttttAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcagt                                                                                                >  2:7116036/1‑101 (MQ=255)

  tttAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGgcagta                                                                                               >  1:7116037/1‑101 (MQ=255)

    taccGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTAcc                                                                                             >  2:7116041/4‑101 (MQ=255)

        gTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGAc                                                                                         >  1:7116044/1‑101 (MQ=255)

        gTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGAc                                                                                         <  1:7115643/101‑1 (MQ=255)

          cAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAg                                                                                       >  2:7116048/1‑101 (MQ=255)

           aGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGc                                                                                      <  1:7115023/101‑1 (MQ=255)

           aGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGc                                                                                      >  1:7116049/1‑101 (MQ=255)

             cACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCa                                                                                    >  1:7116052/1‑101 (MQ=255)

              aCTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAg                                                                                   >  2:7116055/1‑101 (MQ=255)

              aCTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAg                                                                                   <  2:7115893/101‑1 (MQ=255)

               cTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAgg                                                                                  <  2:7115607/101‑1 (MQ=255)

                ttCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGc                                                                                 <  1:7115580/101‑1 (MQ=255)

                ttCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGc                                                                                 >  2:7116058/1‑101 (MQ=255)

                   cTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGa                                                                              >  2:7116062/1‑101 (MQ=255)

                   cTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGa                                                                              <  1:7115695/101‑1 (MQ=255)

                     ggTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAAc                                                                            <  1:7115757/101‑1 (MQ=255)

                     ggTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAAc                                                                            <  1:7115987/101‑1 (MQ=255)

                       tGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAAccc                                                                          >  1:7116065/1‑101 (MQ=255)

                        gATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCg                                                                         <  1:7115524/101‑1 (MQ=255)

                         aTAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGa                                                                        >  2:7116068/1‑101 (MQ=255)

                          tAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAg                                                                       >  2:7116071/1‑101 (MQ=255)

                            gagGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTc                                                                     <  2:7115675/101‑1 (MQ=255)

                            gagGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTc                                                                     >  2:7116074/1‑101 (MQ=255)

                            gagGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTc                                                                     <  2:7115967/101‑1 (MQ=255)

                                ttttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAg                                                                 >  2:7116077/1‑101 (MQ=255)

                                ttttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAg                                                                 >  1:7116078/1‑101 (MQ=255)

                                ttttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAg                                                                 <  1:7115707/101‑1 (MQ=255)

                                 tttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAgg                                                                >  2:7116081/1‑101 (MQ=255)

                                  ttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGc                                                               <  1:7115630/101‑1 (MQ=255)

                                  ttCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGc                                                               <  1:7115798/101‑1 (MQ=255)

                                    caccacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGcac                                                             >  2:7116087/1‑101 (MQ=255)

                                     accacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGcaca                                                            >  1:7116089/1‑101 (MQ=255)

                                     accacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGcaca                                                            <  2:7115466/101‑1 (MQ=255)

                                      ccacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACAt                                                           <  2:7115405/101‑1 (MQ=255)

                                       cacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATc                                                          <  2:7115565/101‑1 (MQ=255)

                                        acGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCt                                                         <  1:7115650/101‑1 (MQ=255)

                                        acGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCt                                                         >  2:7116092/1‑101 (MQ=255)

                                        acGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCt                                                         <  2:7115654/101‑1 (MQ=255)

                                         cGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTg                                                        >  2:7116095/1‑101 (MQ=255)

                                          gTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTgg                                                       >  1:7116097/1‑101 (MQ=255)

                                           ttCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTggg                                                      >  2:7116101/1‑101 (MQ=255)

                                             ccGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCa                                                    >  1:7116105/1‑101 (MQ=255)

                                               ggCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATc                                                  <  1:7115832/101‑1 (MQ=255)

                                                 cAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGt                                                >  1:7116108/1‑101 (MQ=255)

                                                     acaAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTgag                                            >  2:7116112/1‑101 (MQ=255)

                                                     acaAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTgag                                            >  1:7116113/1‑101 (MQ=255)

                                                      caAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTgaga                                           <  1:7115354/101‑1 (MQ=255)

                                                       aaaCACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGAt                                          <  1:7115922/101‑1 (MQ=255)

                                                           accggcaccggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGtt                                      >  2:7116116/1‑101 (MQ=255)

                                                              ggcaccggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATg                                   >  2:7116120/1‑101 (MQ=255)

                                                               gcaccggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGa                                  >  1:7116122/1‑101 (MQ=255)

                                                                caccggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGat                                 <  2:7115831/101‑1 (MQ=255)

                                                                 accggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGata                                >  1:7116128/1‑101 (MQ=255)

                                                                 accggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGata                                >  2:7116127/1‑101 (MQ=255)

                                                                   cggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGa                                  >  1:7116134/1‑97 (MQ=255)

                                                                   cggGAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATAtt                              >  2:7116132/1‑101 (MQ=255)

                                                                    gggAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTc                             >  1:7116136/1‑101 (MQ=255)

                                                                     ggAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCg                            <  1:7115613/101‑1 (MQ=255)

                                                                     ggAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCg                            <  2:7115592/101‑1 (MQ=255)

                                                                      gAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGt                           >  1:7116143/1‑101 (MQ=255)

                                                                      gAAGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGt                           >  2:7116140/1‑101 (MQ=255)

                                                                        aGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGc                         <  1:7115683/101‑1 (MQ=255)

                                                                        aGTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGc                         >  2:7116146/1‑101 (MQ=255)

                                                                         gTCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCa                        >  1:7116149/1‑101 (MQ=255)

                                                                          tCATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCAc                       >  2:7116151/1‑101 (MQ=255)

                                                                           cATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCAcc                      <  2:7115783/101‑1 (MQ=255)

                                                                           cATAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCAcc                      >  1:7116155/1‑101 (MQ=255)

                                                                             tAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGc                    >  2:7116158/1‑101 (MQ=255)

                                                                             tAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGc                    <  1:7115725/101‑1 (MQ=255)

                                                                             tAATCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGc                    <  2:7115742/101‑1 (MQ=255)

                                                                               aTCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCt                  >  2:7116162/1‑101 (MQ=255)

                                                                                tCTTTAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTg                 <  1:7115699/101‑1 (MQ=255)

                                                                                  tttAGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGt               >  2:7116165/1‑101 (MQ=255)

                                                                                     aGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTAc            <  2:7115574/101‑1 (MQ=255)

                                                                                     aGTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTAc            <  1:7115673/101‑1 (MQ=255)

                                                                                      gTGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTAcc           <  1:7115570/101‑1 (MQ=255)

                                                                                       tGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCa          >  1:7116168/1‑101 (MQ=255)

                                                                                       tGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCa          <  2:7115690/101‑1 (MQ=255)

                                                                                       tGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCa          <  2:7115730/101‑1 (MQ=255)

                                                                                       tGACCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCa          <  2:7115888/101‑1 (MQ=255)

                                                                                         aCCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAAt        <  2:7115714/101‑1 (MQ=255)

                                                                                         aCCAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAAt        >  2:7116170/1‑101 (MQ=255)

                                                                                           cAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgc      <  1:7115598/101‑1 (MQ=255)

                                                                                           cAGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgc      >  1:7116172/1‑101 (MQ=255)

                                                                                            aGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcg     >  2:7116175/1‑101 (MQ=255)

                                                                                            aGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcg     <  2:7115799/101‑1 (MQ=255)

                                                                                            aGACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcg     <  1:7115460/101‑1 (MQ=255)

                                                                                             gACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcgg    <  2:7115625/101‑1 (MQ=255)

                                                                                             gACGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcgg    <  1:7115807/101‑1 (MQ=255)

                                                                                             gACGGCAGCACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCt                                                         <  2:7115347/48‑1 (MQ=39)

                                                                                              aCGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcggc   <  1:7115836/101‑1 (MQ=255)

                                                                                              aCGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcggc   <  1:7115738/101‑1 (MQ=255)

                                                                                              aCGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcggc   <  2:7115542/101‑1 (MQ=255)

                                                                                               cGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcggcg  <  1:7115475/101‑1 (MQ=255)

                                                                                               cGGCAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATgcggcg  >  1:7116180/1‑101 (MQ=255)

GTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCG  >  NC_000913/4586308‑4586503

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