BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,142,620	T→C	Q179R (CAG→CGG)	hybA ←	hydrogenase‑2 small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,142,620	0	T	C	100.0%	1243.7 / NA	355	Q179R (CAG→CGG)	hybA	hydrogenase‑2 small subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (187/168); total (187/168)

Only 100 of 370 total aligned reads displayed.

TCTTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCTGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGTAA  >  NC_000913/3142530‑3142709

tctTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCgg                                                                                         >  1:3835009/1‑93 (MQ=255)

tctTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCgg                                                                                         >  1:3480381/1‑93 (MQ=255)

tctTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCgg                                                                                         >  2:7745776/1‑93 (MQ=255)

tctTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCgg                                                                                         >  1:571198/1‑93 (MQ=255)

tctTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGCTCCGCACCTTTCCgg                                                                                         >  2:4276437/1‑93 (MQ=255)

  ttCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGtt                                                                                       >  1:4442069/1‑93 (MQ=255)

  ttCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGtt                                                                                       <  1:5927514/93‑1 (MQ=255)

   tCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTg                                                                                      >  2:3144080/1‑93 (MQ=255)

   tCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTg                                                                                      >  1:507534/1‑93 (MQ=255)

    cACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGc                                                                                     <  2:5986919/93‑1 (MQ=255)

     aCGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGca                                                                                    >  2:4522338/1‑93 (MQ=255)

      cgcgTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGcac                                                                                   <  2:7332067/93‑1 (MQ=255)

       gcgTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGcaca                                                                                  >  1:5833939/1‑93 (MQ=255)

        cgTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAg                                                                                 <  1:6486328/93‑1 (MQ=255)

        cgTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAg                                                                                 <  1:7547008/93‑1 (MQ=255)

         gTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGc                                                                                <  2:4176449/93‑1 (MQ=255)

         gTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGc                                                                                <  2:3704694/93‑1 (MQ=255)

         gTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGc                                                                                >  2:3330937/1‑93 (MQ=255)

         gTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGc                                                                                >  1:8673418/1‑93 (MQ=255)

                aaaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACt                                                                         >  1:3244407/1‑93 (MQ=255)

                aaaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACt                                                                         >  2:8345581/1‑93 (MQ=255)

                aaaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACt                                                                         >  2:7440230/1‑93 (MQ=255)

                aaaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACt                                                                         <  2:3579550/93‑1 (MQ=255)

                 aaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACtt                                                                        >  1:2950776/1‑93 (MQ=255)

                 aaTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACtt                                                                        >  2:757112/1‑93 (MQ=255)

                  aTCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTa                                                                       <  2:490029/93‑1 (MQ=255)

                   tcgccGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTAt                                                                      <  1:1004954/90‑1 (MQ=255)

                   tCACCGCGCCCGCCGGGCAAACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTAt                                                                      >  2:3879947/1‑93 (MQ=255)

                   cccccGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTAt                                                                      <  2:1343029/90‑1 (MQ=255)

                     aCCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGc                                                                    <  2:3461546/93‑1 (MQ=255)

                        gcgcCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAgc                                                                 <  1:2447322/93‑1 (MQ=255)

                        gcgcCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAgc                                                                 <  2:7925549/93‑1 (MQ=255)

                         cgcccgccGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAgcg                                                                >  2:4243279/1‑93 (MQ=255)

                            ccgccgGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCAc                                                             <  2:230623/93‑1 (MQ=255)

                            ccgccgGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTCCCGGTTGCACAGCTCGCACTTATGCAGCGCAc                                                             <  2:2533982/93‑1 (MQ=255)

                                cgGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCaaa                                                         <  2:5007314/93‑1 (MQ=255)

                                 gggCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAAc                                                        <  2:8060635/93‑1 (MQ=255)

                                  ggCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACg                                                       >  1:232037/1‑93 (MQ=255)

                                     ccacTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGgt                                                    <  1:932082/92‑1 (MQ=255)

                                     acacTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGgt                                                    <  2:877661/93‑1 (MQ=255)

                                     acacTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGgt                                                    <  1:1090093/93‑1 (MQ=255)

                                      cacTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGgtt                                                   <  2:6414563/93‑1 (MQ=255)

                                         ttCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGgttgtt                                                <  1:4340555/93‑1 (MQ=255)

                                           cTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTtgttgt                                              <  1:1824093/93‑1 (MQ=255)

                                            tACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTa                                             <  2:1418485/93‑1 (MQ=255)

                                             aCGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAg                                            >  2:2782954/1‑93 (MQ=255)

                                              cGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGt                                           <  1:7949410/93‑1 (MQ=255)

                                              cGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGt                                           <  2:529919/93‑1 (MQ=255)

                                              cGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGt                                           <  2:5840823/93‑1 (MQ=255)

                                                 aGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGt                                        >  2:6762865/1‑93 (MQ=255)

                                                   ccAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTAc                                      <  2:5188833/93‑1 (MQ=255)

                                                    cAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACt                                     <  1:318494/93‑1 (MQ=255)

                                                    cAGGTAGACCGCCTTTATCGAGACGTTCCACACATTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACt                                     <  1:8504275/93‑1 (MQ=255)

                                                     aGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACtt                                    >  2:5623275/1‑93 (MQ=255)

                                                     aGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACtt                                    >  1:715362/1‑93 (MQ=255)

                                                      ggTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTc                                   <  1:3931579/93‑1 (MQ=255)

                                                      ggTAGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTc                                   >  1:2546624/1‑93 (MQ=255)

                                                         aGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGc                                <  2:2652892/93‑1 (MQ=255)

                                                         aGACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGc                                >  1:4506844/1‑93 (MQ=255)

                                                          gACCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCa                               <  1:1490909/93‑1 (MQ=255)

                                                           aCCGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCAc                              >  1:5359276/1‑93 (MQ=255)

                                                             cGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGt                            <  2:7039388/93‑1 (MQ=255)

                                                             cGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGt                            <  2:4376552/93‑1 (MQ=255)

                                                             cGCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGt                            <  1:6568466/93‑1 (MQ=255)

                                                              gCCTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGtt                           <  2:4049219/93‑1 (MQ=255)

                                                                cTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGt                         >  2:5684399/1‑93 (MQ=255)

                                                                cTTTATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGt                         >  2:2242383/1‑93 (MQ=255)

                                                                 tttATCGAGACGTTCTACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTa                        >  1:2782554/1‑93 (MQ=255)

                                                                 tttATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTa                        >  2:6156525/1‑93 (MQ=255)

                                                                 tttATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTa                        >  2:5519244/1‑93 (MQ=255)

                                                                 tttATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGATGTTGTAGTCGTACTTCGGCACGTTGTa                        >  1:603521/1‑93 (MQ=255)

                                                                  ttATCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTAc                       <  1:6157355/93‑1 (MQ=255)

                                                                     tCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGa                    <  2:2147820/93‑1 (MQ=255)

                                                                     tCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGa                    >  1:8790674/1‑93 (MQ=255)

                                                                     tCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGa                    <  1:6178794/93‑1 (MQ=255)

                                                                     tCGAGACGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGa                    <  2:3020621/93‑1 (MQ=255)

                                                                      ccagaCGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCCGCGCCCCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGAc                   <  2:4473743/91‑1 (MQ=255)

                                                                        agaCGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAg                 <  1:6753389/93‑1 (MQ=255)

                                                                         gaCGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGTACGTTGTACGGACAgg                >  1:4034700/1‑93 (MQ=255)

                                                                         gaCGTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAgg                <  2:8734754/93‑1 (MQ=255)

                                                                            gTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGa             >  1:4948940/1‑93 (MQ=255)

                                                                            gTTCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGa             >  1:5586571/1‑93 (MQ=255)

                                                                             ttCCACACCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGAc            >  2:1784149/1‑93 (MQ=255)

                                                                                cacaCCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCAt         <  2:4653963/93‑1 (MQ=255)

                                                                                cacaCCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCAt         <  2:6588368/93‑1 (MQ=255)

                                                                                cacaCCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCAt         >  2:1576189/1‑93 (MQ=255)

                                                                                cacaCCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCAt         >  1:7715286/1‑93 (MQ=255)

                                                                                 ccacCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATg        <  1:1982320/92‑1 (MQ=255)

                                                                                 acacCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATg        >  1:7301293/1‑93 (MQ=255)

                                                                                  cacCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGc       <  1:8397336/93‑1 (MQ=255)

                                                                                  cacCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGc       <  2:7685327/93‑1 (MQ=255)

                                                                                   acCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCa      <  1:8058239/93‑1 (MQ=255)

                                                                                   acCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCa      <  1:2092832/93‑1 (MQ=255)

                                                                                   acCTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCa      >  2:4818485/1‑93 (MQ=255)

                                                                                    ccTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAg     <  2:6045054/93‑1 (MQ=255)

                                                                                    ccTTTCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAg     >  2:5452607/1‑93 (MQ=255)

                                                                                     cTTTCCGGTTGCACAGCTCGCAATTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGt    <  1:6882468/93‑1 (MQ=255)

                                                                                      tttCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGTa   <  1:6997312/93‑1 (MQ=255)

                                                                                       ttCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGTaa  >  1:2171273/1‑93 (MQ=255)

                                                                                       ttCCGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGTaa  >  1:1373559/1‑93 (MQ=255)

TCTTCACGCGTACCGAAAATCACCGCGCCCGCCGGGCACACTTCTACGCAGCCAGGTAGACCGCCTTTATCGAGACGTTCCACACCTTTCTGGTTGCACAGCTCGCACTTATGCAGCGCACCAAACGGGTTGTTGTAGTCGTACTTCGGCACGTTGTACGGACAGGCGACCATGCAGTAA  >  NC_000913/3142530‑3142709

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