BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,225,050	A→G	E1215G (GAG→GGG)	metH →	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,225,050	0	A	G	100.0%	848.4 / NA	257	E1215G (GAG→GGG)	metH	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (131/126); total (131/126)

Only 100 of 265 total aligned reads displayed.

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138

cGGACAGCAAGTACTACGCTGTAGCACAATTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgc                                                                                        <  2:138198/93‑1 (MQ=255)

cGGACAGAAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgc                                                                                        <  1:4369684/93‑1 (MQ=255)

  gACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcgc                                                                                      >  2:3753823/1‑93 (MQ=255)

   aaaGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCGGTTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                     <  1:1091243/91‑1 (MQ=255)

   aCAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                     >  1:2446911/1‑93 (MQ=255)

    cAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAACGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                    >  1:4273156/1‑93 (MQ=255)

      ggaaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGCGCgctggc                                                                                  >  1:6414555/3‑93 (MQ=255)

      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                  >  1:3324802/1‑93 (MQ=255)

      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                  >  2:946032/1‑93 (MQ=255)

        aaGTACGACGCTGTAGCACAAATTGAGCACGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggctg                                                                                <  2:1773987/93‑1 (MQ=255)

         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               >  1:473521/1‑93 (MQ=255)

          gTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              >  1:3124513/1‑93 (MQ=255)

            actacGCTGTAGCACAAATTCAGCGCAATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                            <  2:694460/93‑1 (MQ=255)

             ctacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                           <  2:4531290/93‑1 (MQ=255)

              tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                          >  2:4819315/1‑93 (MQ=255)

                cGCTGTAGCACAAATTAAGCGCAATCAGGTTGAATATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGaa                                                                        <  1:6323445/93‑1 (MQ=255)

                  cgggaGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                      <  2:2479146/89‑1 (MQ=255)

                  cTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                      >  1:1352495/1‑93 (MQ=255)

                  cTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                      >  1:3397863/1‑93 (MQ=255)

                   tGTAGCACAAATTCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     <  2:4763072/93‑1 (MQ=255)

                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                    >  2:2583692/1‑93 (MQ=255)

                     tAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgg                                                                   <  2:6394954/93‑1 (MQ=255)

                      aGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                  >  2:1643704/1‑93 (MQ=255)

                      aGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                  <  1:4141714/93‑1 (MQ=255)

                       acaaaaaTTCAGCGCGATAAGGTTAAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                 <  1:4636607/89‑1 (MQ=255)

                         acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                               >  2:1300295/1‑93 (MQ=255)

                         aaaaaTTCAGCGCTATGAGTTAGAATATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                               <  1:447536/91‑1 (MQ=255)

                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                             >  2:1489824/1‑93 (MQ=255)

                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                             >  1:2829320/1‑93 (MQ=255)

                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            >  1:4860584/1‑93 (MQ=255)

                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                         >  1:825916/1‑93 (MQ=255)

                                 aGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCCCCGAATCTGGGGCATGACGCgg                                                       >  2:4924409/1‑93 (MQ=255)

                                      gATCGGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  <  1:2262006/93‑1 (MQ=255)

                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  <  2:6107156/93‑1 (MQ=255)

                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 >  1:495607/1‑93 (MQ=255)

                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 <  1:5431237/93‑1 (MQ=255)

                                        tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                <  1:1400851/93‑1 (MQ=255)

                                        tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                >  2:795159/1‑93 (MQ=255)

                                        tCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                <  2:5644737/93‑1 (MQ=255)

                                         aaGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                               <  2:5284384/92‑1 (MQ=255)

                                          aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                              >  1:6032909/1‑93 (MQ=255)

                                          aGGTTGAAGATTATGCCCGACGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                              >  1:6550930/1‑93 (MQ=255)

                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAGTCTGGGGTATGAAGCGGACTGATGcac                                             >  2:3827054/1‑93 (MQ=255)

                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  2:2756255/1‑93 (MQ=255)

                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  1:627513/1‑93 (MQ=255)

                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  1:6729268/1‑93 (MQ=255)

                                            gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcaca                                            >  1:1514206/1‑93 (MQ=255)

                                             ttGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaa                                           >  2:2109581/1‑93 (MQ=255)

                                             ttGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaa                                           >  2:714392/1‑93 (MQ=255)

                                             ttGAAAATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaa                                           <  1:3812005/93‑1 (MQ=255)

                                              tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                          >  2:585627/1‑93 (MQ=255)

                                              tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                          >  2:4434105/1‑93 (MQ=255)

                                               gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAt                                         >  1:5151907/1‑93 (MQ=255)

                                                 aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGCAGTTGGGCGCTGGCTGGCACCGAATCTGAGGACTGCCGCGACCTGATTCACAActct                                       >  2:6408148/1‑93 (MQ=255)

                                                 aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCt                                       >  1:2877964/1‑93 (MQ=255)

                                                  gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTg                                      >  1:3235560/1‑93 (MQ=255)

                                                   aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                     >  1:5811508/1‑93 (MQ=255)

                                                   aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGt                                     >  2:4318982/1‑93 (MQ=255)

                                                    ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                    <  1:2968054/93‑1 (MQ=255)

                                                     tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCa                                   <  2:2071595/93‑1 (MQ=255)

                                                     tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCa                                   <  2:6030245/93‑1 (MQ=255)

                                                      ttGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                  <  1:23503/92‑1 (MQ=255)

                                                        gCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtt                                <  2:6458553/93‑1 (MQ=255)

                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              >  2:5215025/1‑93 (MQ=255)

                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  1:2028279/93‑1 (MQ=255)

                                                             ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                           <  1:1645339/93‑1 (MQ=255)

                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                         >  1:1161947/1‑93 (MQ=255)

                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                         <  1:6603351/93‑1 (MQ=255)

                                                                tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaa                        <  2:6239835/93‑2 (MQ=255)

                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAGTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                       >  2:1614673/1‑93 (MQ=255)

                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                       >  2:1556777/1‑93 (MQ=255)

                                                                  aaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaa                      >  2:1076339/1‑93 (MQ=255)

                                                                   aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                     <  2:2002388/93‑1 (MQ=255)

                                                                     gTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAcaacaa                   >  1:3464308/1‑93 (MQ=255)

                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                  >  2:5964181/1‑93 (MQ=255)

                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 <  2:316308/93‑1 (MQ=255)

                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 <  1:2195590/93‑1 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCGTTACaacaaaca                >  1:3852108/1‑93 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                <  1:3623004/93‑1 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                >  2:3909142/1‑93 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                <  2:2213557/93‑1 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                >  1:4504717/1‑93 (MQ=255)

                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                >  1:2617679/1‑93 (MQ=255)

                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               >  2:3652291/1‑93 (MQ=255)

                                                                           gCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAggg             <  1:874792/93‑1 (MQ=255)

                                                                            cGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgc            <  2:575600/93‑1 (MQ=255)

                                                                            cGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgc            <  1:977369/93‑1 (MQ=255)

                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc          <  2:5714483/93‑1 (MQ=255)

                                                                              ttACCGAAGTGGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc          >  1:4071096/1‑93 (MQ=255)

                                                                               tACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGCt         >  2:2245188/1‑93 (MQ=255)

                                                                                aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc        <  2:3355263/93‑1 (MQ=255)

                                                                                aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc        >  1:5331910/1‑93 (MQ=255)

                                                                                 ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa       <  1:6455846/93‑1 (MQ=255)

                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa      >  2:4121556/1‑93 (MQ=255)

                                                                                   caaGTTGGGCGCTGGCTGTCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAGCTGGCACTTTTCCTTACAACAAACAGGGCGCTCAAt     <  1:1822312/92‑1 (MQ=255)

                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg    <  2:4209498/93‑1 (MQ=255)

                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg    >  2:1834546/1‑93 (MQ=255)

                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  <  1:5765578/93‑1 (MQ=255)

                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  >  2:5834607/1‑93 (MQ=255)

                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  <  2:2977920/93‑1 (MQ=255)

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138

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