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

Only 100 of 436 total aligned reads displayed.

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  1:8995000/93‑1 (MQ=255)

gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  1:6523804/93‑1 (MQ=255)

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

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

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

  tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcat                                                                                        >  1:9899360/1‑91 (MQ=255)

   gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                                       <  1:8216932/93‑4 (MQ=255)

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

     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGTCGCATCAGG‑CATCTg                                                                                    >  2:8688770/1‑93 (MQ=255)

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

      aaCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGa                                                                                   <  2:1805777/93‑1 (MQ=255)

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

        cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGc                                                                                 >  1:2155438/1‑93 (MQ=255)

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

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

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

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

           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATCGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatg                                                                              <  1:3826224/93‑1 (MQ=255)

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

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

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

             ttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatgat                                                                            >  2:9053956/1‑93 (MQ=255)

               gCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGc                                                                          >  2:2684711/1‑93 (MQ=255)

                  tGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCatga                                                                             >  2:2176077/1‑87 (MQ=255)

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

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

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

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

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

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

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

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

                         cTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggc                                                                >  2:6973864/1‑93 (MQ=255)

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

                          tGGGAAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCATTTACGGCGCAGCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                               >  1:7073/1‑93 (MQ=37)

                           gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                              <  2:3184049/93‑2 (MQ=255)

                           gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                              <  1:2916914/93‑2 (MQ=255)

                           gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                              <  1:3333490/93‑2 (MQ=255)

                             gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgctg                                                            >  2:858703/1‑90 (MQ=255)

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

                                acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                         >  1:6256167/1‑93 (MQ=25)

                                 caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                        <  1:659312/93‑1 (MQ=25)

                                 caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                        <  2:7874885/93‑1 (MQ=25)

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

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

                                    aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                     >  1:3537377/1‑93 (MQ=25)

                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    >  2:4304886/1‑93 (MQ=21)

                                       cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGTCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                  >  1:9431477/1‑93 (MQ=21)

                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                <  1:9488079/93‑1 (MQ=21)

                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                <  1:1825819/93‑1 (MQ=21)

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

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

                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATcg                                              >  2:7764541/1‑92 (MQ=18)

                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              <  1:9701206/93‑1 (MQ=21)

                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              <  2:3686645/93‑1 (MQ=21)

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

                                                  attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       >  1:6927576/1‑93 (MQ=17)

                                                  attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       >  1:852204/1‑93 (MQ=17)

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

                                                    tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGTTGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCTACGa                                     >  2:8111594/1‑93 (MQ=17)

                                                     gCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCCGGCCTACGAc                                    >  1:2498755/1‑93 (MQ=17)

                                                      tttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                   <  2:1871221/92‑1 (MQ=17)

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

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

                                                         ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                <  1:7875585/93‑1 (MQ=18)

                                                         ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGACGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                <  1:9339560/93‑1 (MQ=17)

                                                          tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                               >  1:5009243/1‑93 (MQ=18)

                                                          tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                               >  2:295324/1‑93 (MQ=18)

                                                          tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCACGCCTACGACTTtct                               >  2:4817594/1‑93 (MQ=18)

                                                           gTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                              >  2:635910/1‑93 (MQ=18)

                                                            tAGGCTGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                             <  2:3755332/93‑1 (MQ=17)

                                                             aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                            >  2:5153008/1‑93 (MQ=18)

                                                                ccGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                         >  2:9698449/1‑93 (MQ=18)

                                                                 cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTTTCCGCATc                        >  1:428612/1‑93 (MQ=17)

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

                                                                 cGGATAAGGCGTTTACGCCGCATCA‑GCCATCTGAGCATGATGCATGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                        <  1:9805949/93‑1 (MQ=18)

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

                                                                  ggATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       >  1:4514964/1‑93 (MQ=18)

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

                                                                    aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     >  2:9517052/1‑93 (MQ=21)

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

                                                                     tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  2:8415279/1‑93 (MQ=21)

                                                                     tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  1:8796611/1‑93 (MQ=21)

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

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

                                                                        ggCGTTTACGCCTCATCAGG‑CATCTGAGCATGATGCCTGATGCGTCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                 <  2:653792/93‑1 (MQ=18)

                                                                        ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                 <  2:3469917/93‑1 (MQ=21)

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

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

                                                                           gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac              >  1:4724903/1‑93 (MQ=25)

                                                                           gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac              <  1:4115005/93‑1 (MQ=21)

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

                                                                            tgtACGCCGCATCAGG‑CAACTGAGCATGAGTCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             <  1:2323143/91‑1 (MQ=21)

                                                                               aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg          >  2:4014049/1‑93 (MQ=25)

                                                                               aCGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg          <  2:4118477/93‑1 (MQ=25)

                                                                                 gccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg        <  1:1569107/93‑1 (MQ=25)

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

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

                                                                                      gtCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa   <  1:5103577/92‑1 (MQ=37)

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

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

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