BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,321,028	A→G	R228R (CGT→CGC)	mrsA ←	MrsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,321,028	0	A	G	100.0%	627.6 / NA	201	R228R (CGT→CGC)	mrsA	MrsA
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (126/75); total (126/75)

Only 100 of 217 total aligned reads displayed.

ATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGC  >  NC_000913/3320938‑3321117

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        >  1:1039779/1‑93 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        >  2:5617/1‑93 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        <  2:3889569/93‑1 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        <  2:333736/93‑1 (MQ=255)

 tCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGa                                                                                        >  1:75481/1‑93 (MQ=255)

    aCTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGcc                                                                                     >  1:1209487/1‑93 (MQ=255)

    aCTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGcc                                                                                     >  1:4261230/1‑93 (MQ=255)

     cTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCt                                                                                    <  1:3822696/93‑1 (MQ=255)

     cTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCt                                                                                    >  2:1242558/1‑93 (MQ=255)

      tttATTGGCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTg                                                                                   >  2:4457413/1‑93 (MQ=255)

      tttATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTg                                                                                   >  1:3240258/1‑93 (MQ=255)

      tttATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTg                                                                                   >  1:3653203/1‑93 (MQ=255)

         aTTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAg                                                                                >  2:2153792/1‑93 (MQ=255)

           tGCCTTCATTGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcg                                                                              <  2:4376036/93‑1 (MQ=255)

            gCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgc                                                                             >  1:3424488/1‑93 (MQ=255)

             ccTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgcg                                                                            <  2:2341384/93‑1 (MQ=255)

             ccTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgcg                                                                            >  1:2499219/1‑93 (MQ=255)

              cTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgcgc                                                                           >  2:4214401/1‑93 (MQ=255)

              cTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAgcgcgc                                                                           >  1:1182310/1‑93 (MQ=255)

               ttCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGcgcgcg                                                                          >  2:1122859/1‑93 (MQ=255)

               ttCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAACACGCGAGCCTGGAGcgcgcg                                                                          >  2:1521826/1‑93 (MQ=255)

                tCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGa                                                                         >  2:3230971/1‑93 (MQ=255)

                  aTGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAAc                                                                       <  2:2538973/93‑1 (MQ=255)

                   tGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGTGGCGCGCGAACg                                                                      >  2:3990765/1‑93 (MQ=255)

                   tGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACg                                                                      >  2:1961233/1‑93 (MQ=255)

                   tGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACg                                                                      >  2:1021267/1‑93 (MQ=255)

                    ggTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACgc                                                                     >  2:618969/1‑92 (MQ=255)

                    ggTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGt                                                                     <  1:4175437/93‑1 (MQ=255)

                    ggTCAACCATAATCACGCGATCGCCATCGCCGGCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGt                                                                     >  1:2814690/1‑93 (MQ=255)

                     gTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTc                                                                    <  1:2074280/93‑1 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTa                                                                <  2:1571025/93‑1 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTa                                                                >  2:2084365/1‑93 (MQ=255)

                         aCCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGACTGGAGCGCGCGAACGTCGGTa                                                                >  1:2302258/1‑93 (MQ=255)

                          ccATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAg                                                               >  1:2181675/1‑93 (MQ=255)

                           cATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGc                                                              >  2:2905453/1‑93 (MQ=255)

                           cATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGc                                                              >  2:2227773/1‑93 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccc                                                            >  1:3547716/1‑93 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccc                                                            >  2:3164621/1‑93 (MQ=255)

                             tAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGCAGccc                                                            >  2:480878/1‑93 (MQ=255)

                              aaTCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGcccc                                                           >  1:2688995/1‑93 (MQ=255)

                               aTCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccccc                                                          >  2:1488731/1‑93 (MQ=255)

                               aTCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGccccc                                                          >  1:4475141/1‑93 (MQ=255)

                                  ccgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACt                                                       <  1:1276332/92‑1 (MQ=255)

                                  aCGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACt                                                       >  1:3997676/1‑93 (MQ=255)

                                    gcaatcgccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTc                                                     <  2:2496152/90‑1 (MQ=255)

                                       atcgccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGc                                                  >  2:1772364/1‑93 (MQ=255)

                                        tccccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCa                                                 <  2:2707928/90‑1 (MQ=255)

                                          gccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCAtt                                               <  2:1635182/93‑1 (MQ=255)

                                          cccatcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCAtt                                               <  2:3799816/92‑1 (MQ=255)

                                            catcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGCAGCCCCCACTTCGGCATTGa                                             >  1:1785854/1‑93 (MQ=255)

                                             atcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGAt                                            >  1:1364658/1‑93 (MQ=255)

                                                cccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGtt                                         <  2:4015317/92‑1 (MQ=255)

                                                 ccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGttt                                        >  2:3909901/1‑93 (MQ=255)

                                                   gTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTac                                      <  2:412970/93‑1 (MQ=255)

                                                   gTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAAGCCCCACTTCGGCATTGATGTTTac                                      >  2:3942686/1‑93 (MQ=255)

                                                    tCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTaca                                     >  2:4347480/1‑93 (MQ=255)

                                                     cGAAGGCAATGCCGAGATCCGCTTTTACAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTacac                                    <  1:2233567/93‑1 (MQ=255)

                                                     cGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTacac                                    >  1:3142845/1‑93 (MQ=255)

                                                       aaGGCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCg                                  >  2:665447/1‑93 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                >  1:2318097/1‑93 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                >  1:2942758/1‑93 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                <  2:725535/93‑1 (MQ=255)

                                                         ggCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGtt                                <  1:775764/93‑1 (MQ=255)

                                                          gCAATACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGttt                               >  2:357012/1‑93 (MQ=255)

                                                              tACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCt                           <  1:3342829/93‑1 (MQ=255)

                                                              tACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCt                           >  2:3529963/1‑93 (MQ=255)

                                                              tACCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCt                           >  1:4586334/1‑93 (MQ=255)

                                                               aCCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTc                          <  2:3669821/93‑1 (MQ=255)

                                                               aCCGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTc                          >  1:200807/1‑93 (MQ=255)

                                                                ccGAGATCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTca                         <  2:50303/93‑1 (MQ=255)

                                                                  gagaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTcaca                       <  2:3384876/93‑1 (MQ=255)

                                                                   agaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACaa                      <  1:1193560/93‑1 (MQ=255)

                                                                    gaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAc                     >  2:1342098/1‑93 (MQ=255)

                                                                    gaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAc                     <  2:1046123/93‑1 (MQ=255)

                                                                    gaTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAc                     <  2:1902205/93‑1 (MQ=255)

                                                                     gtCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAcc                    <  1:170437/92‑1 (MQ=255)

                                                                     aTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAcc                    >  2:3816127/1‑93 (MQ=255)

                                                                     aTCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAAcc                    >  1:1692488/1‑93 (MQ=255)

                                                                      tCCGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCg                   >  2:4649839/1‑93 (MQ=255)

                                                                        cGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGAt                 >  1:4431814/1‑93 (MQ=255)

                                                                        cGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGAt                 >  1:4520030/1‑93 (MQ=255)

                                                                        cGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGAt                 >  1:4945/1‑93 (MQ=255)

                                                                        cGCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGAt                 >  2:1215561/1‑93 (MQ=255)

                                                                         gCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATa                <  2:110009/93‑1 (MQ=255)

                                                                         gCTTTTTCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGCTTACACCGTTTGGCTCACAACCGATa                >  1:2970532/1‑93 (MQ=255)

                                                                           tttttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGc              <  1:37115/93‑1 (MQ=255)

                                                                           tttttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGc              <  1:1572647/93‑1 (MQ=255)

                                                                           tttttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGc              <  1:3489234/93‑1 (MQ=255)

                                                                             tttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGa            <  2:4535213/93‑1 (MQ=255)

                                                                             tttCAGCCAGCACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGa            <  1:2600886/93‑1 (MQ=255)

                                                                                 agccagcGCGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAAc        <  2:2289101/93‑1 (MQ=255)

                                                                                 agccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAAc        >  1:3828488/1‑93 (MQ=255)

                                                                                  gccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACg       >  2:3237270/1‑93 (MQ=255)

                                                                                  gccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACg       >  1:245619/1‑93 (MQ=255)

                                                                                   ccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGt      <  1:886777/93‑1 (MQ=255)

                                                                                   ccagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGt      >  1:3883971/1‑93 (MQ=255)

                                                                                    cagcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGtt     >  2:2636397/1‑93 (MQ=255)

                                                                                    cagcAACCGAGCCTGGAGCGCACGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACTCCGTTTGGCTCAAAACCGATAGCGATAACGtt     <  2:4582971/93‑1 (MQ=255)

                                                                                     agcACGCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTc    >  2:129082/1‑93 (MQ=255)

                                                                                       cacgCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGc  >  2:820070/1‑93 (MQ=255)

ATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGC  >  NC_000913/3320938‑3321117

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