BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,051,144	Δ1 bp	intergenic (+146/+305)	hemN → / ← glnG	O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,051,143	0	C	.	100.0%	635.9 / NA	191	intergenic (+145/+306)	hemN/glnG	O2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (99/92); total (99/92)

Only 100 of 207 total aligned reads displayed.

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         >  1:1028197/1‑92 (MQ=255)

 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         <  1:3683626/93‑2 (MQ=255)

 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         <  2:2572712/93‑2 (MQ=255)

     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTg                                                                                    >  2:1416667/1‑93 (MQ=255)

       aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAg                                                                                  <  2:1833232/93‑1 (MQ=255)

          aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCat                                                                               <  2:1977725/93‑1 (MQ=255)

            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatga                                                                             >  2:4433525/1‑93 (MQ=255)

            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatga                                                                             <  1:2716299/93‑1 (MQ=255)

            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatga                                                                             <  1:2171615/93‑1 (MQ=255)

             ttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatgat                                                                            <  1:1072642/93‑1 (MQ=255)

             ttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatgat                                                                            <  2:2719950/93‑1 (MQ=255)

                ccTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGcc                                                                         <  2:1121583/93‑1 (MQ=255)

                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCt                                                                        <  2:1621910/93‑1 (MQ=255)

                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCt                                                                        <  1:3493228/93‑1 (MQ=255)

                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCt                                                                        >  1:4653790/1‑93 (MQ=255)

                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCt                                                                        <  2:596567/93‑1 (MQ=255)

                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGa                                                                      >  2:169561/1‑93 (MQ=255)

                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGa                                                                      <  2:3309166/93‑1 (MQ=255)

                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGAt                                                                     <  2:2773324/93‑1 (MQ=255)

                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGAt                                                                     <  1:4069775/93‑1 (MQ=255)

                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGAt                                                                     <  1:3729722/93‑1 (MQ=255)

                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGAt                                                                     <  1:3862115/93‑1 (MQ=255)

                     tacgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATg                                                                    >  2:4617994/3‑93 (MQ=255)

                     tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATg                                                                    <  1:4494521/93‑1 (MQ=255)

                     tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATg                                                                    >  2:3204987/1‑93 (MQ=255)

                      gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgc                                                                   <  1:2337849/93‑1 (MQ=255)

                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcgg                                                                 <  1:3227246/93‑1 (MQ=255)

                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcgg                                                                 >  1:4631718/1‑93 (MQ=255)

                         cTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc                                                                <  1:4574444/93‑1 (MQ=255)

                          tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                               >  1:4298133/1‑93 (MQ=255)

                          tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                               <  2:1500528/93‑1 (MQ=255)

                             gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgctg                                                            <  1:3083123/93‑4 (MQ=255)

                              caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgctgg                                                           <  2:353699/93‑5 (MQ=255)

                              caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgctgg                                                           >  2:1324592/1‑89 (MQ=255)

                                acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                         <  1:4172080/93‑1 (MQ=25)

                                acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                         <  2:4377667/93‑1 (MQ=25)

                                 caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                        <  1:4099495/93‑1 (MQ=21)

                                  aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                       >  1:4274364/1‑93 (MQ=25)

                                   aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                      >  1:1380139/1‑93 (MQ=21)

                                    aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                     >  2:920241/1‑93 (MQ=21)

                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcggc                                                    >  1:1805492/1‑91 (MQ=21)

                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    <  2:754875/93‑1 (MQ=25)

                                       cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGCTGCGGCGCTCGTCGCGCCtt                                                  >  2:2379501/1‑93 (MQ=21)

                                       cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                  >  1:2828445/1‑93 (MQ=21)

                                        tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                 <  1:3309407/93‑1 (MQ=21)

                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                >  1:580885/1‑93 (MQ=21)

                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                >  1:1262400/1‑93 (MQ=21)

                                          ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                               <  1:1023688/93‑1 (MQ=21)

                                            gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                             <  2:3148945/93‑1 (MQ=21)

                                            gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                             >  2:2903942/1‑93 (MQ=18)

                                              aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                           <  1:2975417/93‑1 (MQ=18)

                                               attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                          >  2:505422/1‑93 (MQ=18)

                                                 tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                        <  1:2021160/93‑1 (MQ=18)

                                                  attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       <  1:1478071/93‑1 (MQ=17)

                                                  attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       <  1:2221574/93‑1 (MQ=18)

                                                    tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                     <  1:4460937/93‑1 (MQ=18)

                                                    tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                     <  1:2624497/93‑1 (MQ=18)

                                                      cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                   >  2:2540494/1‑93 (MQ=17)

                                                        tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                 <  2:724380/93‑1 (MQ=17)

                                                         ttGTAGGCCGAATAAGACTTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGGGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                <  2:4107031/93‑1 (MQ=17)

                                                          tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                               <  2:3464426/93‑1 (MQ=18)

                                                              ggCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                           >  2:78473/1‑93 (MQ=18)

                                                                 cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                        >  2:2053078/1‑93 (MQ=18)

                                                                 cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                        <  2:4556511/93‑1 (MQ=18)

                                                                  ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       <  1:1860882/93‑1 (MQ=18)

                                                                  ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       >  1:4107601/1‑93 (MQ=18)

                                                                  ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       <  1:1249823/93‑1 (MQ=18)

                                                                   gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      <  1:1066792/93‑1 (MQ=18)

                                                                   gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  2:3661323/1‑93 (MQ=18)

                                                                    aTAAGGCGTTTTGGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     <  2:1573182/93‑1 (MQ=18)

                                                                    aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     >  2:2836241/1‑93 (MQ=18)

                                                                    aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     <  2:1532933/93‑1 (MQ=18)

                                                                     tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  2:2641630/1‑93 (MQ=18)

                                                                      aaGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCa                   <  2:1234197/93‑1 (MQ=18)

                                                                      aaGGCGTTTACGCCCCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCa                   <  1:1350429/93‑1 (MQ=18)

                                                                       aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  2:2622347/1‑93 (MQ=21)

                                                                       aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  2:824330/1‑93 (MQ=21)

                                                                        ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                 >  2:119442/1‑93 (MQ=21)

                                                                         gCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                >  1:4367036/1‑93 (MQ=21)

                                                                         gCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                >  2:2789750/1‑93 (MQ=21)

                                                                          cGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa               <  1:1121679/93‑1 (MQ=21)

                                                                           gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac              >  2:3726761/1‑93 (MQ=21)

                                                                            tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             <  2:3129742/93‑1 (MQ=21)

                                                                            tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:2244885/1‑93 (MQ=25)

                                                                            tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  1:1173247/1‑93 (MQ=25)

                                                                            tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:4319252/1‑93 (MQ=21)

                                                                            tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCAGCATAaaca             >  2:3763695/1‑93 (MQ=21)

                                                                            tttACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:4404444/1‑93 (MQ=21)

                                                                             ttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            <  2:3347570/93‑1 (MQ=21)

                                                                              tACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  2:3990942/93‑1 (MQ=25)

                                                                              tACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  1:768906/93‑1 (MQ=25)

                                                                              tACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  1:4545802/93‑1 (MQ=25)

                                                                                   cgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt      <  1:4670449/93‑1 (MQ=25)

                                                                                     cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    >  1:2800949/1‑93 (MQ=37)

                                                                                     cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    >  1:478079/1‑93 (MQ=37)

                                                                                     cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    >  2:804455/1‑93 (MQ=37)

                                                                                      aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa   <  2:2996275/93‑1 (MQ=37)

                                                                                      aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa   >  2:3878399/1‑93 (MQ=37)

                                                                                       tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  >  2:2421598/1‑93 (MQ=37)

                                                                                       tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  >  1:3416583/1‑93 (MQ=37)

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

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: `‑`