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%	1012.2 / NA	219	coding (541/609 nt)	b4556	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (90/129); total (90/129)

Only 100 of 248 total aligned reads displayed.

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

  cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                        <  1:100367/93‑3 (MQ=255)

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

    tcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                      <  2:2276536/93‑5 (MQ=255)

    tcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                      >  2:5554176/1‑89 (MQ=255)

    tcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                      >  2:5349805/1‑89 (MQ=255)

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

       ttGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCgg                                                                                  >  1:4941860/1‑93 (MQ=255)

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

        tGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGa                                                                                 <  1:5409285/93‑1 (MQ=255)

         gCTTGACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                                                <  1:6247134/93‑1 (MQ=255)

          cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAcg                                                                               >  1:4995296/1‑92 (MQ=255)

           tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                                              <  1:393796/93‑1 (MQ=255)

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

           tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                                              <  2:2721157/93‑1 (MQ=255)

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

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

             tACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTggg                                                                            >  2:2669390/1‑92 (MQ=255)

              aCCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGac                                                                           >  1:3201078/1‑92 (MQ=255)

              aCCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAg                                                                           <  1:3645373/93‑1 (MQ=255)

               ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGc                                                                          >  1:1101482/1‑93 (MQ=255)

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

               ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGAATTACAGGTTGGGTTTGCGATGAA‑CGACGCGGAATgtcgc                                                                          >  1:1886911/1‑89 (MQ=255)

                cGGCGCACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCt                                                                         <  1:830222/93‑1 (MQ=255)

                    acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCgg                                                                     >  2:2010044/1‑93 (MQ=255)

                      acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGa                                                                   >  2:3955344/1‑93 (MQ=255)

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

                          gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                               >  2:1360828/1‑93 (MQ=255)

                          gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                               >  1:2472740/1‑93 (MQ=255)

                          gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                                               >  1:4104621/1‑93 (MQ=255)

                           aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACttt                                                              <  1:4570491/93‑1 (MQ=255)

                            agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                                             <  2:4688816/93‑1 (MQ=255)

                            agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                                             <  2:110284/93‑1 (MQ=255)

                             gagaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTGCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTg                                                            <  1:5353236/93‑1 (MQ=255)

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

                               gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGcc                                                          <  1:2606168/93‑1 (MQ=255)

                                aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                                         <  2:6148294/93‑1 (MQ=255)

                                aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                                         >  1:2779740/1‑93 (MQ=255)

                                aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                                         >  2:3373181/1‑93 (MQ=255)

                                aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                                         >  1:576512/1‑93 (MQ=255)

                                  gCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                                       >  1:3666211/1‑93 (MQ=255)

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

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

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

                                    aGGCGGGGGGGGCGCATGGAATTTATATCGGCGTTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                     <  1:2247967/93‑1 (MQ=255)

                                      gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                                   >  2:5520927/1‑93 (MQ=255)

                                      gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                                   <  1:297084/93‑1 (MQ=255)

                                        gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                                 >  1:1985011/1‑93 (MQ=255)

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

                                         tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                >  1:603891/1‑93 (MQ=255)

                                         tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                <  1:6195934/93‑1 (MQ=255)

                                         tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                >  1:9157/1‑93 (MQ=255)

                                         tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                                                >  2:1846663/1‑93 (MQ=255)

                                            gtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                                             >  2:1750828/1‑93 (MQ=255)

                                            gtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                                             >  1:3061075/1‑93 (MQ=255)

                                            gtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                                             <  1:1482255/93‑1 (MQ=255)

                                             tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                                            >  2:4898109/1‑93 (MQ=255)

                                             tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                                            >  2:3756576/1‑93 (MQ=255)

                                              ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                                           >  2:3598825/1‑93 (MQ=255)

                                               gcgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaa                                          <  2:1282412/93‑1 (MQ=255)

                                                cgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaa                                         >  2:423946/1‑93 (MQ=255)

                                                cgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGAACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaa                                         <  1:1234331/93‑1 (MQ=255)

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

                                                  cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                       <  1:4669696/93‑1 (MQ=255)

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

                                                    tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                     <  2:2550294/93‑1 (MQ=255)

                                                     ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                    >  1:3115258/1‑93 (MQ=255)

                                                     ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                    >  2:3515785/1‑93 (MQ=255)

                                                      gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                   >  1:779515/1‑93 (MQ=255)

                                                      gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                   >  2:1467836/1‑93 (MQ=255)

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

                                                       aaTTTATATCGGCATTACCGGTTGGGTTGGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                                  <  1:1546185/93‑1 (MQ=255)

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

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                >  1:6092493/1‑93 (MQ=255)

                                                         tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                                >  2:159573/1‑93 (MQ=255)

                                                          ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                               >  2:2935540/1‑93 (MQ=255)

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

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

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

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

                                                              atCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaa                           >  2:3266399/1‑93 (MQ=255)

                                                               tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa                          <  1:6084095/93‑1 (MQ=255)

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

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

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

                                                                 ggCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt                        >  2:819180/1‑93 (MQ=255)

                                                                 ggCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt                        <  2:3678920/93‑1 (MQ=255)

                                                                   cATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                      >  2:2882100/1‑93 (MQ=255)

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

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

                                                                       aCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcg                  >  2:4733179/1‑93 (MQ=255)

                                                                       aCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgcg                  >  1:474245/1‑93 (MQ=255)

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

                                                                            ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTGTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTa             >  1:3378812/1‑93 (MQ=255)

                                                                              gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc           <  2:594887/93‑1 (MQ=255)

                                                                              gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc           <  1:4246554/93‑1 (MQ=255)

                                                                                 tttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATGCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGc        <  1:2925968/93‑1 (MQ=255)

                                                                                  ttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCACCGTATCTGCt       <  2:446381/93‑1 (MQ=255)

                                                                                   tGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCTGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc      <  1:4440750/93‑1 (MQ=255)

                                                                                    gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTcc     <  2:3017410/93‑1 (MQ=255)

                                                                                       aTGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc  <  2:4379687/93‑1 (MQ=255)

TTCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTC  >  NC_000913/4021591‑4021771

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