BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,021,684	Δ1 bp	coding (542/609 nt)\|coding (1/243 nt)	b4556 → tatD →	hypothetical protein Mg‑dependent DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,021,683	0	C	.	100.0%	1107.7 / NA	239	coding (541/609 nt)	b4556	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (113/126); total (113/126)

Only 100 of 261 total aligned reads displayed.

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          <  2:178212/93‑1 (MQ=255)

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          <  2:520536/93‑1 (MQ=255)

ttcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                                          <  1:1755747/93‑1 (MQ=255)

  cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                        >  2:5043472/1‑91 (MQ=255)

     cATTGCTTTACCGGCACACGCGAAGAGGTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgc                                                                                    <  2:1376384/93‑1 (MQ=255)

     cATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgc                                                                                    <  2:542642/93‑1 (MQ=255)

      aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                                   >  1:4582068/1‑93 (MQ=255)

      aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                                   >  1:3863498/1‑93 (MQ=255)

       ttGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCgg                                                                                  <  2:2681812/93‑1 (MQ=255)

         tcTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                <  1:6132741/92‑1 (MQ=255)

         gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                <  1:6639659/93‑1 (MQ=255)

         gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                >  1:4022863/1‑93 (MQ=255)

          cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACt                                                                               >  2:3187423/1‑93 (MQ=255)

           tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                                              >  2:5087387/1‑93 (MQ=255)

           tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                                              >  2:3492254/1‑93 (MQ=255)

            ttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTGCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTgg                                                                             <  2:3963016/93‑1 (MQ=255)

            ttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTgg                                                                             <  1:2087497/93‑1 (MQ=255)

               ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGc                                                                          <  2:5801509/93‑1 (MQ=255)

                cGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCt                                                                         <  1:5672546/93‑1 (MQ=255)

                cGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCt                                                                         <  1:1366425/93‑1 (MQ=255)

                 ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACGGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTg                                                                        <  2:6285947/93‑1 (MQ=255)

                 ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTg                                                                        >  1:2045585/1‑93 (MQ=255)

                 ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTg                                                                        <  2:2541182/93‑1 (MQ=255)

                  gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGATTGCGATGAA‑CGACGCGGACTGGAGCTga                                                                       >  2:6435327/1‑92 (MQ=255)

                   cacaCGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                                      >  2:5023574/1‑93 (MQ=255)

                       cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGaa                                                                  <  1:3196734/93‑1 (MQ=255)

                        gcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAAc                                                                 <  2:2812918/93‑1 (MQ=255)

                         cgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACt                                                                <  1:4253974/93‑1 (MQ=255)

                          gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                               <  1:1569580/93‑1 (MQ=255)

                           aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGCCGCGGACTGGAGCTGCGGGCACttt                                                              >  1:6075135/1‑93 (MQ=255)

                            agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                                             >  1:2741223/1‑93 (MQ=255)

                            agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                                             >  2:4615991/1‑93 (MQ=255)

                             gagaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTg                                                            >  2:6116534/1‑93 (MQ=255)

                             gagaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTg                                                            <  2:5122089/93‑1 (MQ=255)

                              agaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                                           <  2:6240221/93‑1 (MQ=255)

                              agaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                                           >  1:6298666/1‑93 (MQ=255)

                                 tGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                                        <  2:1093817/93‑1 (MQ=255)

                                    aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     <  2:1868316/93‑1 (MQ=255)

                                    aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     >  1:1628358/1‑93 (MQ=255)

                                     ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                    <  1:3133441/93‑1 (MQ=255)

                                     ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                    <  1:503282/93‑1 (MQ=255)

                                        gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                                 <  1:5369275/93‑1 (MQ=255)

                                          gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcgg                                               >  2:4871339/1‑93 (MQ=255)

                                           cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                              <  1:4744944/93‑1 (MQ=255)

                                           cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                                              <  2:4280635/93‑1 (MQ=255)

                                                 gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                        <  1:750483/93‑1 (MQ=255)

                                                  cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       >  1:2939007/1‑93 (MQ=255)

                                                  cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       >  1:2506865/1‑93 (MQ=255)

                                                   aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                      >  2:1599834/1‑93 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     >  2:3632506/1‑93 (MQ=255)

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     >  1:2587096/1‑93 (MQ=255)

                                                      gcaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTat                                   >  1:1491288/3‑92 (MQ=255)

                                                       aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                                  <  1:5163481/93‑1 (MQ=255)

                                                       aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                                  >  2:6605286/1‑93 (MQ=255)

                                                        aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                 <  1:3325022/93‑1 (MQ=255)

                                                        aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                 >  2:685789/1‑93 (MQ=255)

                                                        aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                 >  2:3553000/1‑93 (MQ=255)

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                <  2:1663888/93‑1 (MQ=255)

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                <  2:533145/93‑1 (MQ=255)

                                                          ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                               <  1:5801746/93‑1 (MQ=255)

                                                           tataTCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc                              <  2:4122777/93‑1 (MQ=255)

                                                           tataTCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc                              >  1:1960380/1‑93 (MQ=255)

                                                           tataTCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc                              >  2:2644012/1‑93 (MQ=255)

                                                            atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg                             >  1:1002964/1‑93 (MQ=255)

                                                            atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg                             <  1:3460629/93‑1 (MQ=255)

                                                             tatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGa                            >  1:4093707/1‑93 (MQ=255)

                                                             tatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGa                            >  2:4363084/1‑93 (MQ=255)

                                                               tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTATTGATCGaaa                          >  1:3947343/1‑93 (MQ=255)

                                                               tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa                          >  1:4869985/1‑93 (MQ=255)

                                                               tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa                          >  1:1123952/1‑93 (MQ=255)

                                                                cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                         >  1:2914977/1‑93 (MQ=255)

                                                                cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                         <  1:4715833/93‑1 (MQ=255)

                                                                cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                         <  1:6171281/93‑1 (MQ=255)

                                                                  gCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCTGCGGAAAAAGTACTGATTGAAACTg                       >  2:3869039/1‑93 (MQ=255)

                                                                  gCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTg                       >  1:2983019/1‑93 (MQ=255)

                                                                  gCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTg                       >  2:2057953/1‑93 (MQ=255)

                                                                   cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                      <  1:3596762/93‑1 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCTGGCGGAAAAACTACTGATCCAAACTGATg                    >  2:5496378/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    <  2:1301190/93‑1 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  1:972964/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  2:5227276/1‑93 (MQ=255)

                                                                     ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg                    >  1:1237057/1‑93 (MQ=255)

                                                                       tccGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcg                  <  2:1538936/92‑1 (MQ=255)

                                                                           gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt              >  2:3936439/1‑93 (MQ=255)

                                                                              gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc           >  2:2187430/1‑93 (MQ=255)

                                                                              gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc           >  2:5852537/1‑93 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt          <  2:5261465/93‑1 (MQ=255)

                                                                               ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt          <  1:3264717/93‑1 (MQ=255)

                                                                                gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTg         >  2:5003128/1‑93 (MQ=255)

                                                                                  ttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCt       <  2:4232511/93‑1 (MQ=255)

                                                                                   tGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc      >  1:1081194/1‑93 (MQ=255)

                                                                                   tGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc      <  1:4478771/93‑1 (MQ=255)

                                                                                     cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc    <  2:3785686/93‑1 (MQ=255)

                                                                                     cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc    <  2:2434954/93‑1 (MQ=255)

                                                                                     cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc    <  1:6498850/93‑1 (MQ=255)

                                                                                     cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc    <  1:4969455/93‑1 (MQ=255)

                                                                                     cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc    <  2:827539/93‑1 (MQ=255)

                                                                                      gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCt   <  1:5005870/93‑1 (MQ=255)

                                                                                      gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCt   >  1:4215607/1‑93 (MQ=255)

                                                                                       aTGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc  >  2:5684357/1‑93 (MQ=255)

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

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