BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,246,188	A→G	E213E (GAA→GAG)	lamB →	phage lambda receptor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,246,188	0	A	G	100.0%	888.0 / NA	273	E213E (GAA→GAG)	lamB	phage lambda receptor protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (133/140); total (133/140)

Only 100 of 280 total aligned reads displayed.

CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277

caaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCt                                                                                         >  1:5462410/1‑93 (MQ=255)

 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        <  1:3013673/93‑1 (MQ=255)

 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCAGGGCGGCACATTAGAGCTg                                                                                        >  2:2081811/1‑93 (MQ=255)

   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      <  2:5379974/93‑1 (MQ=255)

   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      >  1:5616731/1‑93 (MQ=255)

   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      <  1:3308869/93‑1 (MQ=255)

   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      >  2:4002136/1‑93 (MQ=255)

    atTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgt                                                                                     <  1:4548222/93‑1 (MQ=255)

      ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   <  2:543086/93‑1 (MQ=255)

       tATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTc                                                                                  >  2:1950654/1‑93 (MQ=255)

        aTGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCg                                                                                 <  2:3011569/93‑1 (MQ=255)

          gACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAc                                                                               >  2:4947482/1‑93 (MQ=255)

           aCTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACt                                                                              >  1:1118454/1‑93 (MQ=255)

             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            >  2:2150601/1‑93 (MQ=255)

              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                           >  2:4024427/1‑93 (MQ=255)

               taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                          >  2:631519/1‑93 (MQ=255)

                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                         <  1:2574675/93‑1 (MQ=255)

                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                         <  2:5267457/93‑1 (MQ=255)

                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                         <  1:5298037/93‑1 (MQ=255)

                 ccAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        <  2:1623189/93‑1 (MQ=255)

                 ccAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        >  1:973835/1‑93 (MQ=255)

                 ccAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        >  1:4377083/1‑93 (MQ=255)

                 acAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        <  1:607753/92‑1 (MQ=255)

                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      >  2:4883605/1‑93 (MQ=255)

                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      >  1:648674/1‑93 (MQ=255)

                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                      >  2:1362368/1‑93 (MQ=255)

                    aCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     >  1:4889588/1‑93 (MQ=255)

                      gAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                   >  1:3146960/1‑93 (MQ=255)

                       aaaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCa                                                                  >  1:600541/1‑93 (MQ=255)

                        aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCac                                                                 >  1:4166432/1‑92 (MQ=255)

                        aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCaa                                                                 >  1:765466/1‑93 (MQ=255)

                         aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                >  2:723194/1‑93 (MQ=255)

                         aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGGGTCGAATACGGTCGTGcacac                                                                >  2:1552752/1‑89 (MQ=255)

                          ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               <  2:4904753/93‑1 (MQ=255)

                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             >  1:3784673/1‑93 (MQ=255)

                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             >  2:3125999/1‑93 (MQ=255)

                             cgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGc                                                            <  1:4769162/93‑1 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  2:4524186/1‑93 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          <  2:4438889/93‑1 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          <  2:5479231/93‑1 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  2:6253720/1‑93 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  1:4005136/1‑93 (MQ=255)

                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          <  1:5054229/93‑1 (MQ=255)

                                acgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTg                                                         >  1:1637354/1‑93 (MQ=255)

                                  gacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGAt                                                       <  1:2164206/93‑1 (MQ=255)

                                  gacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGAt                                                       <  1:6208824/93‑1 (MQ=255)

                                  gacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGAt                                                       >  2:410312/1‑93 (MQ=255)

                                    cgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATaa                                                     >  2:2352285/1‑93 (MQ=255)

                                     gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    >  1:27069/1‑93 (MQ=255)

                                      ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                   >  1:3574170/1‑93 (MQ=255)

                                      ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                   >  2:3723882/1‑93 (MQ=255)

                                        ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 <  2:3498155/93‑1 (MQ=255)

                                        ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 <  2:838204/93‑1 (MQ=255)

                                         tCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATc                                                <  2:2866213/93‑1 (MQ=255)

                                         tCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATc                                                >  2:5383605/1‑93 (MQ=255)

                                         tCGATGGGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATc                                                <  2:4716798/93‑1 (MQ=255)

                                           gATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGt                                              >  2:5189780/1‑93 (MQ=255)

                                            gtgtgCGTTTAGCGCAGATGGAAATCAGCCCGGGCGGCAGATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:854577/92‑1 (MQ=255)

                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:6132113/93‑1 (MQ=255)

                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGGCACCTTGCGTGATAACTATCGTCt                                            >  2:2601569/1‑93 (MQ=255)

                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                            <  1:1367805/93‑1 (MQ=255)

                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                            >  1:1870393/1‑93 (MQ=255)

                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                            >  2:224794/1‑93 (MQ=255)

                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                          <  1:141406/93‑1 (MQ=255)

                                                gCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGt                                         >  2:950024/1‑93 (MQ=255)

                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAGTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  2:2513811/1‑93 (MQ=255)

                                                  gTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTg                                       >  1:2271396/1‑93 (MQ=255)

                                                    ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  2:5776594/93‑1 (MQ=255)

                                                    ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     >  2:6093448/1‑93 (MQ=255)

                                                      aGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGATTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATgg                                   >  2:3213873/1‑93 (MQ=255)

                                                       gcgcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgc                                  <  2:4554427/93‑1 (MQ=255)

                                                        cgcAGATGGAAATCAACCCGGGCGGCAAATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcg                                 <  2:5154300/93‑1 (MQ=255)

                                                         gcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcgc                                >  1:5120851/1‑93 (MQ=255)

                                                           aGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCAt                              >  1:373041/1‑93 (MQ=255)

                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             <  1:3903573/93‑1 (MQ=255)

                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             >  1:1012650/1‑93 (MQ=255)

                                                            gATCGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTCCGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             <  1:2623083/93‑1 (MQ=255)

                                                               ggAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCgta                          >  2:563395/1‑91 (MQ=255)

                                                               ggAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaa                          >  2:4143504/1‑93 (MQ=255)

                                                                  aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  1:1323632/93‑1 (MQ=255)

                                                                  aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  2:1183219/93‑1 (MQ=255)

                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  1:6237005/93‑1 (MQ=255)

                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  1:1678345/93‑1 (MQ=255)

                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                      <  2:6003962/93‑1 (MQ=255)

                                                                    tCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACg                     <  1:3058105/93‑1 (MQ=255)

                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   <  1:5486233/93‑1 (MQ=255)

                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                   >  2:3353763/1‑93 (MQ=255)

                                                                       aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                  >  1:4139509/1‑93 (MQ=255)

                                                                       aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                  <  2:4229377/93‑1 (MQ=255)

                                                                        cccGGGCGGTACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTg                 <  2:3852028/93‑1 (MQ=255)

                                                                        cccGGGCGGCACATTAGAGCTGGGTGTCCACTACGGTCGTGCCAACTTGCGTGGTAACTATCGTCTGGGTGATGGCGCATCGCAAAACGGCTg                 >  2:5857936/1‑93 (MQ=255)

                                                                         ccGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTgg                <  1:4866334/93‑1 (MQ=255)

                                                                           gggCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGtt              <  2:5014629/93‑1 (MQ=255)

                                                                             gcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTAt            <  2:1775067/93‑1 (MQ=255)

                                                                                gcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCa         <  1:2474609/93‑1 (MQ=255)

                                                                                 cacaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAc        <  1:2377605/93‑1 (MQ=255)

                                                                                 cacaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAc        <  1:3685791/93‑1 (MQ=255)

                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  2:2636539/93‑1 (MQ=255)

                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       <  2:1110916/93‑1 (MQ=255)

                                                                                   caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg      <  1:5811021/93‑1 (MQ=255)

CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246098‑4246277

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