BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,051,170	+C	intergenic (+172/+279)	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,170	1	.	C	100.0%	417.6 / NA	151	intergenic (+172/+279)	hemN/glnG	O2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (90/61); total (90/61)

Only 100 of 163 total aligned reads displayed.

TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC  >  NC_000913/4051077‑4051255

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

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

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

       caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                                                 >  2:332717/1‑93 (MQ=25)

        aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                <  1:2419866/93‑1 (MQ=25)

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

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

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

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

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

             cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                           >  2:4644048/1‑93 (MQ=21)

              tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATTCCTGATGCGGCGCTGGTCGCGTCTTa                                                                          >  1:4526028/1‑93 (MQ=21)

              tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                                          >  1:1675235/1‑93 (MQ=21)

                ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                                        >  2:2173172/1‑93 (MQ=21)

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

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

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

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

                   aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                                                     >  1:590833/1‑93 (MQ=18)

                    aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                    >  1:2283751/1‑93 (MQ=18)

                     attattGCTTTGGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                   <  1:1303906/93‑1 (MQ=18)

                      tttttGTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                  <  2:3626102/90‑1 (MQ=17)

                      ttattTTTTTGGTAGGCCGAATAAGGCTTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                  <  2:3626103/93‑1 (MQ=17)

                      ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                  <  1:2725699/93‑1 (MQ=18)

                        attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                                >  2:3175637/1‑93 (MQ=18)

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

                        attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                                >  2:1163243/1‑93 (MQ=18)

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

                         ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                                               <  2:4637773/93‑1 (MQ=17)

                           gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                                             >  2:3729947/1‑93 (MQ=17)

                            gttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                            <  2:255828/92‑1 (MQ=17)

                            cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                            <  2:255829/93‑1 (MQ=17)

                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                           >  1:2643272/1‑93 (MQ=17)

                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                           <  1:4189293/93‑1 (MQ=17)

                               ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                                         >  2:2716982/1‑93 (MQ=18)

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

                                  tAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                      <  1:2927793/93‑1 (MQ=18)

                                   aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                                     <  1:1998263/93‑1 (MQ=18)

                                   aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGACTACGACTTTCTCCg                                                     >  2:1707862/1‑93 (MQ=18)

                                     gCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                                   <  1:2527844/93‑1 (MQ=18)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                   ttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                                     >  1:3333703/1‑93 (MQ=25)

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

                                                        gcgcATCAGG‑CATCTTAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                <  1:3933284/92‑1 (MQ=25)

                                                          gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                              >  2:3417207/1‑93 (MQ=37)

                                                          gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                              <  2:4501118/93‑1 (MQ=25)

                                                          gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                              >  2:3051145/1‑93 (MQ=37)

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

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

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

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

                                                            aTCAGG‑CATCTGAGCAGGAGGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                            <  1:4565877/93‑1 (MQ=25)

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

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

                                                             tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGACTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                           >  2:187820/1‑93 (MQ=37)

                                                               caggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                          >  2:3534963/5‑93 (MQ=38)

                                                                 ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                        >  2:663990/3‑93 (MQ=255)

                                                                  gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa                       >  2:2901779/2‑93 (MQ=255)

                                                                   cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa                      <  1:4592056/93‑1 (MQ=255)

                                                                    aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                     <  2:3395679/93‑1 (MQ=255)

                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  1:2755044/1‑93 (MQ=255)

                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  1:1403398/1‑93 (MQ=255)

                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                  <  1:3139770/93‑1 (MQ=255)

                                                                        gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                 >  2:221350/1‑93 (MQ=255)

                                                                        gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                 >  1:3246683/1‑93 (MQ=255)

                                                                         aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                >  1:513923/1‑93 (MQ=255)

                                                                          gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCa               >  1:4319235/1‑93 (MQ=255)

                                                                          gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCa               >  1:3112097/1‑93 (MQ=255)

                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc            <  2:3284584/93‑1 (MQ=255)

                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc            <  2:1483808/93‑1 (MQ=255)

                                                                              gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCa           >  2:2156861/1‑93 (MQ=255)

                                                                              gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCa           >  1:73267/1‑93 (MQ=255)

                                                                                tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt         <  2:2431461/93‑1 (MQ=255)

                                                                                agCCGGATGCGGCGCTGGTCGCGTCTTATCAAGCCTTCGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt         <  2:2111428/92‑1 (MQ=38)

                                                                                 gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTg        <  2:4175646/93‑1 (MQ=255)

                                                                                  ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc       >  2:2502772/1‑93 (MQ=255)

                                                                                    tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg     >  2:3878665/1‑93 (MQ=255)

                                                                                     gATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGc    <  1:1513609/93‑1 (MQ=255)

                                                                                      aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGca   <  1:1876690/93‑1 (MQ=255)

                                                                                       tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac  >  1:3510887/1‑93 (MQ=255)

                                                                                       tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac  >  2:2989296/1‑93 (MQ=255)

TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC  >  NC_000913/4051077‑4051255

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