BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,538,045	G→A	P857L (CCG→CTG)	narZ ←	cryptic nitrate reductase 2 alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,538,045	0	G	A	100.0%	830.6 / NA	213	P857L (CCG→CTG)	narZ	cryptic nitrate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (116/97); total (116/97)

Only 100 of 216 total aligned reads displayed.

GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCT  >  NC_000913/1537955‑1538134

gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                         >  1:2812726/1‑93 (MQ=255)

gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                         >  1:6010249/1‑93 (MQ=255)

gCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAga                                                                                         >  2:2032670/1‑93 (MQ=255)

  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                       >  2:594080/1‑93 (MQ=255)

  cacaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAgaga                                                                                       <  2:5074569/93‑1 (MQ=255)

   acaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCagagag                                                                                      <  1:6053795/93‑1 (MQ=255)

   acaATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCagagag                                                                                      <  2:1924989/93‑1 (MQ=255)

      atatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCa                                                                                   >  2:1173489/1‑93 (MQ=255)

      atatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCa                                                                                   >  2:1539968/1‑93 (MQ=255)

       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCAt                                                                                  <  2:840021/93‑1 (MQ=255)

       tatTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCAt                                                                                  >  2:6438061/1‑93 (MQ=255)

           cGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAg                                                                              <  2:6193111/93‑1 (MQ=255)

           cGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAg                                                                              <  1:2533051/93‑1 (MQ=255)

            ggCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGt                                                                             <  1:5518920/93‑1 (MQ=255)

             gcgcGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTa                                                                            <  2:5712404/93‑1 (MQ=255)

             gcgcGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTa                                                                            <  1:1627322/93‑1 (MQ=255)

              cgcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAg                                                                           >  2:5255113/1‑93 (MQ=255)

               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                          >  2:1081240/1‑93 (MQ=255)

               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                          >  1:331391/1‑93 (MQ=255)

               gcgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgg                                                                          <  1:1024478/93‑1 (MQ=255)

                cgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         >  1:1556767/1‑93 (MQ=255)

                cgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         <  1:829903/93‑1 (MQ=255)

                cgGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         <  1:2059940/93‑1 (MQ=255)

                cgGTGTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAggg                                                                         <  1:2182904/93‑1 (MQ=255)

                 ggTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAgggg                                                                        <  1:5121536/93‑1 (MQ=255)

                       gCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGATGTaa                                                                  >  1:2943567/1‑93 (MQ=255)

                        ccTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCTTGCAGTAGGGGTTGTAAt                                                                 <  1:4078111/93‑1 (MQ=255)

                        ccTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAAt                                                                 <  1:1776714/93‑1 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAAtc                                                                >  2:4777538/1‑92 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                >  2:4901202/1‑93 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                >  1:5026015/1‑93 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                <  1:2383567/93‑1 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                <  2:1949510/93‑1 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                >  1:35655/1‑93 (MQ=255)

                         cTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATa                                                                <  2:4587955/93‑1 (MQ=255)

                           ggAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATAcc                                                              >  2:2903183/1‑93 (MQ=255)

                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCa                                                             <  2:134072/93‑1 (MQ=255)

                            gAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCa                                                             >  1:6360116/1‑93 (MQ=255)

                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            <  2:4130612/93‑1 (MQ=255)

                             aaTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCAc                                                            >  2:4505186/1‑93 (MQ=255)

                              aTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                           >  2:1402346/1‑93 (MQ=255)

                              aTCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                           >  1:3139346/1‑93 (MQ=255)

                               tCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGt                                                          >  2:4283035/1‑93 (MQ=255)

                                cAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTc                                                         <  1:1175899/93‑1 (MQ=255)

                                  gATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCgg                                                       >  1:3629198/1‑93 (MQ=255)

                                  gATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCgg                                                       <  2:639410/93‑1 (MQ=255)

                                   aTCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGt                                                      <  1:2371245/93‑1 (MQ=255)

                                    tCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGtt                                                     <  1:411812/93‑1 (MQ=255)

                                     cGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGttt                                                    >  2:2478177/1‑93 (MQ=255)

                                      gCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTc                                                   <  2:6112094/93‑1 (MQ=255)

                                      gCATTCGCCTATGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTc                                                   >  1:3929209/1‑93 (MQ=255)

                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGGTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                >  2:3443940/1‑93 (MQ=255)

                                         ttCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCttt                                                >  2:4376668/1‑93 (MQ=255)

                                           cGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGc                                              >  1:4886805/1‑93 (MQ=255)

                                           cGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGc                                              <  2:4603565/93‑1 (MQ=255)

                                            gCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCCCTGCCAGAGAGTCATGCAGTAGGGGGTGTAATACCACGTCGGTTTCTTTGcc                                             >  1:1636411/1‑93 (MQ=255)

                                            cccTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGcc                                             <  1:6183097/92‑1 (MQ=255)

                                              cTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCaa                                           >  1:3743304/1‑93 (MQ=255)

                                               tttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAg                                          <  2:3756308/93‑1 (MQ=255)

                                               tttGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAg                                          >  1:4563622/1‑93 (MQ=255)

                                                 tGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGAt                                        >  1:4449704/1‑93 (MQ=255)

                                                 tGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGAt                                        <  2:5179229/93‑1 (MQ=255)

                                                 tGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGAt                                        >  1:5991630/1‑93 (MQ=255)

                                                  gcgcCAGTCGAGCACTTCACACGTCTGTGAGAGCTCCGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATg                                       <  1:4276637/93‑1 (MQ=255)

                                                    gcCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAc                                     >  1:569231/1‑93 (MQ=255)

                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    >  2:4934632/1‑93 (MQ=255)

                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    >  1:1423285/1‑93 (MQ=255)

                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    >  1:4422822/1‑93 (MQ=255)

                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    <  1:1683178/93‑1 (MQ=255)

                                                     ccAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAcc                                    >  1:5976147/1‑93 (MQ=255)

                                                        gTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCAc                                 >  1:2200608/1‑93 (MQ=255)

                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                <  1:1902095/93‑1 (MQ=255)

                                                         tCGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACg                                <  2:3565916/93‑1 (MQ=255)

                                                          cGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGc                               <  1:4400141/93‑1 (MQ=255)

                                                          cGAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATACAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGc                               >  1:2798458/1‑93 (MQ=255)

                                                           gAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGca                              >  1:1927343/1‑93 (MQ=255)

                                                           gAGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGca                              >  1:3068378/1‑93 (MQ=255)

                                                            aGCACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGcac                             >  2:2582643/1‑93 (MQ=255)

                                                              cACTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCacac                           <  1:3423462/93‑1 (MQ=255)

                                                               aCTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACt                          <  2:3412443/93‑1 (MQ=255)

                                                               aCTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACt                          <  2:3849105/93‑1 (MQ=255)

                                                                cTTCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACtt                         >  2:1237558/1‑93 (MQ=255)

                                                                 ttCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTg                        >  2:5385534/1‑93 (MQ=255)

                                                                  tCACACGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGc                       <  1:593771/93‑1 (MQ=255)

                                                                    acacGGCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGa                     >  2:2718401/1‑93 (MQ=255)

                                                                        ggCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAAt                 >  1:5162944/1‑93 (MQ=255)

                                                                        ggCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAAt                 >  2:3140610/1‑93 (MQ=255)

                                                                         gCTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATg                <  1:3382376/93‑1 (MQ=255)

                                                                          cTGTGAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGc               <  2:2531440/93‑1 (MQ=255)

                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            >  2:3260012/1‑93 (MQ=255)

                                                                             tgAGAGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCttt            <  1:533236/93‑1 (MQ=255)

                                                                              gagaGCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTg           <  1:2476602/93‑1 (MQ=255)

                                                                               agagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTgg          <  2:475996/93‑1 (MQ=255)

                                                                                gagCTCTGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGc         <  1:4717673/93‑1 (MQ=255)

                                                                                     ctGCCAGAGAGTCATGCAGTAGGGTTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc    <  1:4953083/93‑1 (MQ=255)

                                                                                     ctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc    <  2:2619663/93‑1 (MQ=255)

                                                                                     ctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc    <  2:1781867/93‑1 (MQ=255)

                                                                                     ctGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAc    >  2:3765372/1‑93 (MQ=255)

                                                                                      tGCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAcc   >  1:6268947/1‑93 (MQ=255)

                                                                                       gCCAGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCt  >  2:5153545/1‑93 (MQ=255)

GCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCT  >  NC_000913/1537955‑1538134

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