BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	100.0%	401.9 / NA	128	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (54/74); total (54/74)

Only 100 of 149 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT  >  NC_000913/3553996‑3554173

cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                                         >  1:1054102/1‑93 (MQ=255)

 ggCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATACAATCCG‑a                                                                                       >  1:2158082/1‑93 (MQ=255)

  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       <  2:2069232/93‑3 (MQ=255)

  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                       >  1:4576262/1‑91 (MQ=255)

    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGGGGATAAAATCCgcagg                                                                                     >  2:428194/1‑89 (MQ=255)

      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                   <  1:1797099/93‑1 (MQ=255)

       gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                                  <  2:4170605/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                 <  1:3873226/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                 <  1:1095835/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                 <  2:2175031/93‑1 (MQ=255)

         aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                >  1:3469756/1‑93 (MQ=255)

            atccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                             <  2:118443/91‑1 (MQ=255)

              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                           <  1:262325/93‑1 (MQ=255)

               cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                          <  1:2838140/93‑1 (MQ=255)

               cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                          <  1:2026583/93‑1 (MQ=255)

                gTGGCGGAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                         <  1:2345617/93‑1 (MQ=255)

                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                        <  1:3783596/93‑1 (MQ=255)

                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                        >  2:1871512/1‑93 (MQ=255)

                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACGCAGGCGGCGCCGACggg                                                                        >  1:2088208/1‑93 (MQ=255)

                    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                                     <  1:3624893/93‑1 (MQ=255)

                    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                                     >  1:4181941/1‑93 (MQ=255)

                     gTAACAGCGTGGTTTGCTGATTAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                    >  1:850619/1‑93 (MQ=255)

                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                   >  2:4069735/1‑93 (MQ=255)

                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                   <  2:3555355/93‑1 (MQ=255)

                       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                  >  1:252449/1‑93 (MQ=255)

                        aCAGCGTGGTTTGCTGATGGATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCAGCTGGCGGCGCCGACGGGTTATctg                                                                 >  2:4056528/1‑91 (MQ=255)

                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                 >  2:919895/1‑93 (MQ=255)

                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCgg                                                                 >  2:922595/1‑93 (MQ=255)

                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCggg                                                                >  2:1038126/1‑92 (MQ=255)

                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                                >  2:4202339/1‑93 (MQ=255)

                          aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACGCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                               <  1:3585710/93‑2 (MQ=255)

                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCggggc                                                              >  1:4118909/1‑90 (MQ=255)

                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                              >  2:1621241/1‑91 (MQ=255)

                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                              >  1:3209630/1‑91 (MQ=255)

                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                         <  2:1531334/93‑1 (MQ=255)

                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                        <  2:4439807/93‑2 (MQ=255)

                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                        >  2:4063874/1‑92 (MQ=255)

                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                        >  2:220201/1‑92 (MQ=255)

                                 tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                                      <  2:1247371/93‑4 (MQ=255)

                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                  <  1:750/93‑1 (MQ=255)

                                       gattaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCAGT‑GCCGCCGCTCACtt                                                 <  1:4043451/93‑1 (MQ=255)

                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                 <  1:2766717/93‑1 (MQ=255)

                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                 <  1:2766718/93‑1 (MQ=255)

                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                <  1:2160305/93‑1 (MQ=255)

                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGCGTATCGGTGCCGG‑CGCTCACTTc                                                >  2:3435127/1‑93 (MQ=255)

                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                              >  1:3119911/1‑93 (MQ=255)

                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                              <  1:3828318/93‑1 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                            >  2:3283508/1‑93 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                            <  2:2404314/93‑1 (MQ=255)

                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                           <  1:2435959/93‑1 (MQ=255)

                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                          <  2:222069/93‑1 (MQ=255)

                                               ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCCCCCTTCACCAcgc                                         >  2:2450606/1‑93 (MQ=255)

                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                        <  1:4229768/93‑1 (MQ=255)

                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                        >  2:3314139/1‑93 (MQ=255)

                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                       <  1:4340818/93‑1 (MQ=255)

                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCCCACTTCAACACGCGa                                       >  1:1242478/1‑93 (MQ=255)

                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      <  2:3443757/93‑1 (MQ=255)

                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                      <  1:1810139/93‑1 (MQ=255)

                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                     >  2:3047164/1‑93 (MQ=255)

                                                    tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                    <  1:2282234/93‑1 (MQ=255)

                                                    tttCGCCAGCAGCCGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                    <  1:2282232/93‑1 (MQ=255)

                                                     ttCGCCAGCAGCGGCGCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCTCTCACTTCAACACGCGAAGGt                                   <  2:1312088/93‑1 (MQ=255)

                                                      tCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGTGATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                  <  2:423878/93‑1 (MQ=255)

                                                       cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                 <  1:995531/93‑1 (MQ=255)

                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                <  2:1442330/93‑1 (MQ=255)

                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGCCCa                                >  2:4429037/1‑93 (MQ=255)

                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                               <  1:161103/93‑1 (MQ=255)

                                                         ccAGCAGCGGCTCCAGCACCCGGAGCACAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                               <  2:2725390/93‑1 (MQ=255)

                                                         ccAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGTGCCAt                               >  1:1005292/1‑93 (MQ=255)

                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                             >  1:4213923/1‑93 (MQ=255)

                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                             >  1:3755294/1‑93 (MQ=255)

                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                            <  1:4345312/93‑1 (MQ=255)

                                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                           <  2:1326308/93‑1 (MQ=255)

                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                         >  2:4642465/1‑93 (MQ=255)

                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                         >  2:72531/1‑93 (MQ=255)

                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                         <  1:2525006/93‑1 (MQ=255)

                                                                cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                        <  1:3014300/93‑1 (MQ=255)

                                                                 ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                       >  2:641783/1‑93 (MQ=255)

                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                      <  2:4336504/93‑1 (MQ=255)

                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                      <  2:3824765/93‑1 (MQ=255)

                                                                   cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                     <  1:2979699/93‑1 (MQ=255)

                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                   >  1:1388452/1‑93 (MQ=255)

                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                   <  2:3629497/93‑1 (MQ=255)

                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                   >  2:2256637/1‑93 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                  <  1:282591/93‑1 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                  >  1:4560097/1‑93 (MQ=255)

                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc               <  1:3731424/93‑1 (MQ=255)

                                                                          aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg              >  1:3513648/1‑93 (MQ=255)

                                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAACCCACAGCGCgg             <  2:4676865/93‑1 (MQ=255)

                                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag         <  2:2660219/93‑1 (MQ=255)

                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACCGCGCGGgcagca       >  2:3127015/1‑93 (MQ=255)

                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagcc       >  2:1433318/1‑92 (MQ=255)

                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCCC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc      <  1:1685437/93‑1 (MQ=255)

                                                                                    tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg    >  1:3893443/1‑93 (MQ=255)

                                                                                    tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg    >  1:837997/1‑93 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGCCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt   <  2:3583466/93‑1 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt   <  1:2722032/93‑1 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt   <  1:4089410/93‑1 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt   <  2:4030535/93‑1 (MQ=255)

                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt  <  2:4350820/93‑1 (MQ=255)

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTT  >  NC_000913/3553996‑3554173

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