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%	939.8 / NA	282	intergenic (+145/+306)	hemN/glnG	O2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (155/127); total (155/127)

Only 100 of 304 total aligned reads displayed.

CTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGC‑G‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051052‑4051230

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

  gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                                        <  2:3210958/93‑4 (MQ=255)

   agaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatct                                                                                       >  1:3063515/1‑89 (MQ=255)

    gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTg                                                                                     <  1:8526341/93‑1 (MQ=255)

     aaCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTgg                                                                                    >  2:5761651/1‑92 (MQ=255)

      aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAg                                                                                   >  2:3057322/1‑93 (MQ=255)

       cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGc                                                                                  <  1:7397364/93‑1 (MQ=255)

        aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCa                                                                                 >  1:5707189/1‑93 (MQ=255)

        aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCa                                                                                 <  1:8710514/93‑1 (MQ=255)

         aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGCTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCat                                                                                >  1:1166376/1‑93 (MQ=255)

         aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCat                                                                                >  1:4590959/1‑93 (MQ=255)

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

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

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

                cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAGGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCt                                                                         <  2:6110128/93‑1 (MQ=255)

                 tGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTg                                                                        >  1:2452154/1‑93 (MQ=255)

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

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

                   aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGAt                                                                      >  1:8638208/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

                        cTGTGCAACAAAGCCTGTTGAAATTATTGCTTTTGTAGTCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc                                                                 <  1:5054408/93‑1 (MQ=25)

                          gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc‑gc                                                              >  1:3391057/1‑92 (MQ=255)

                          gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc‑gc                                                              <  1:5375618/93‑2 (MQ=255)

                          gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc‑gc                                                              <  2:6622215/93‑2 (MQ=255)

                            gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc‑gctg                                                            >  2:7554699/1‑90 (MQ=255)

                            gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc‑gctg                                                            >  1:723732/1‑90 (MQ=255)

                               acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTc                                                         >  2:1171493/1‑93 (MQ=21)

                               acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTc                                                         >  1:8395528/1‑93 (MQ=21)

                                caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTcg                                                        >  2:4813793/1‑93 (MQ=25)

                                 aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTcgc                                                       >  1:3554677/1‑93 (MQ=25)

                                 aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTcgc                                                       >  1:7442560/1‑93 (MQ=25)

                                 aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTcgc                                                       >  2:7757538/1‑93 (MQ=25)

                                  aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTcgcg                                                      >  1:4959022/1‑93 (MQ=21)

                                   aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGt                                                     >  1:294391/1‑93 (MQ=25)

                                    aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTc                                                    >  2:2462086/1‑93 (MQ=21)

                                     ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCt                                                   <  1:4668367/93‑1 (MQ=21)

                                     ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCt                                                   <  2:229610/93‑1 (MQ=21)

                                      cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCtt                                                  >  2:3493408/1‑93 (MQ=21)

                                      cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCtt                                                  >  1:164985/1‑93 (MQ=21)

                                      cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCtt                                                  <  1:6679912/93‑1 (MQ=21)

                                      cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCtt                                                  >  2:3857893/1‑93 (MQ=21)

                                       tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTa                                                 <  2:1627686/93‑1 (MQ=21)

                                        attGAAATTAATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCC‑TCTGAGCATGATTCCTGATTCGGC‑GCTGGGCGGGTTTTAt                                                <  2:1343151/92‑1 (MQ=18)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               <  2:8426770/93‑1 (MQ=21)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               <  1:2349543/93‑1 (MQ=21)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               >  1:8073153/1‑93 (MQ=21)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               >  1:4205749/1‑93 (MQ=21)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               >  1:5559293/1‑93 (MQ=21)

                                         ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTCAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATc                                               >  2:6790824/1‑93 (MQ=21)

                                          tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCa                                              <  2:1848932/93‑1 (MQ=21)

                                              attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGcc                                          >  1:986686/1‑93 (MQ=18)

                                                 attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTAc                                       >  2:1955083/1‑93 (MQ=18)

                                                  ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTCCg                                      >  2:2580445/1‑93 (MQ=17)

                                                  ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACg                                      <  2:4894660/93‑1 (MQ=18)

                                                    gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGACTACGAc                                    >  2:5296076/1‑93 (MQ=17)

                                                       tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACttt                                 >  1:249623/1‑93 (MQ=17)

                                                          gTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                              <  1:692328/93‑1 (MQ=18)

                                                              gCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                          <  1:2007214/93‑1 (MQ=18)

                                                                 ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       >  2:8206804/1‑93 (MQ=18)

                                                                 ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       >  2:1455240/1‑93 (MQ=18)

                                                                  gATAAGGCGTTTACGTCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      <  1:5858604/93‑1 (MQ=18)

                                                                  gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      <  2:8365941/93‑1 (MQ=18)

                                                                  gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  1:3881458/1‑93 (MQ=18)

                                                                  gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  1:7938447/1‑93 (MQ=18)

                                                                  gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  1:1525293/1‑93 (MQ=18)

                                                                   aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     <  1:646461/93‑1 (MQ=18)

                                                                   aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     <  2:8232671/93‑1 (MQ=18)

                                                                   aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     >  1:4895122/1‑93 (MQ=18)

                                                                    tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  2:3383766/1‑93 (MQ=21)

                                                                    tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  2:8094725/1‑93 (MQ=21)

                                                                     aaGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCa                   <  2:2346727/93‑1 (MQ=18)

                                                                      aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  1:7064707/1‑93 (MQ=21)

                                                                      aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  2:7391734/1‑93 (MQ=21)

                                                                       ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                 >  1:196522/1‑93 (MQ=21)

                                                                         cGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa               <  2:7963341/93‑1 (MQ=21)

                                                                         cGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa               <  1:4014980/93‑1 (MQ=21)

                                                                           tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:865234/1‑93 (MQ=25)

                                                                           tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:3910655/1‑93 (MQ=25)

                                                                            ttccgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            <  2:6279336/90‑1 (MQ=21)

                                                                            ttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            <  2:435767/93‑1 (MQ=21)

                                                                             tACGCCGCATCAGG‑CATCTGAGCATGATGCGTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  2:714705/93‑1 (MQ=25)

                                                                               cgccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa         >  2:4302549/1‑93 (MQ=25)

                                                                                 ccgcATCAGG‑CATCTGAGCCGGATGCCTGATGCGGCTC‑TGGTCGCGTCTTATCAGGCCTACGACTTCCTCCGCATCCGGCATAAACAACGAGc       <  1:8306210/93‑1 (MQ=21)

                                                                                 ccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc       >  2:5834595/1‑93 (MQ=37)

                                                                                   gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     >  2:2397219/1‑93 (MQ=37)

                                                                                   gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     >  1:4279029/1‑93 (MQ=37)

                                                                                   gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     >  2:4562207/1‑93 (MQ=37)

                                                                                   gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     <  2:4434448/93‑1 (MQ=25)

                                                                                    cATCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    <  2:638547/93‑1 (MQ=37)

                                                                                      tCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGGGCTTCAg  >  1:2599184/1‑93 (MQ=37)

                                                                                      tCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  >  2:7877603/1‑93 (MQ=37)

                                                                                      tCAGG‑CATCTGAGCATGATGCCTGATGCGGC‑GCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  <  2:2737660/93‑1 (MQ=37)

CTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGC‑G‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051052‑4051230

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