BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,026,252	2 bp→GC	coding (110‑111/1164 nt)	fadA ←	thiolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,026,252	0	C	G	100.0%	764.6 / NA	237	T37S (ACG→AGC)	fadA	thiolase I
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (125/112); total (125/112)
*	NC_000913	4,026,253	0	G	C	100.0%	774.2 / NA	242	T37S (ACG→AGC)	fadA	thiolase I
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (126/116); total (126/116)

Only 100 of 245 total aligned reads displayed.

AAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCA  >  NC_000913/4026166‑4026339

||

aaaCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGc                                                                                   >  2:645180/1‑93 (MQ=255)

 aaCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCa                                                                                  <  2:1043154/93‑1 (MQ=255)

  aCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                                 <  2:4252268/93‑1 (MQ=255)

  aCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                                 >  2:1397902/1‑93 (MQ=255)

  aCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                                 >  2:2848498/1‑93 (MQ=255)

  aCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAt                                                                                 <  1:332085/93‑1 (MQ=255)

   cccTGCTCCAGCGTCTGCTGCACACACCCCCAGTAAATATCGTCGGGGCCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAtt                                                                                <  2:3511198/93‑1 (MQ=255)

   cccTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCAtt                                                                                >  1:2290795/1‑93 (MQ=255)

    ccTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTa                                                                               >  2:2634860/1‑93 (MQ=255)

     cTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaa                                                                              >  1:84562/1‑93 (MQ=255)

      cgccccACCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                                             <  1:2014380/89‑1 (MQ=255)

       gCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAAt                                                                            >  1:473115/1‑93 (MQ=255)

        cTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATg                                                                           >  1:4830911/1‑93 (MQ=255)

         tCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          <  2:301242/93‑1 (MQ=255)

         tCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          <  1:5159355/93‑1 (MQ=255)

         tCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          >  1:2340371/1‑93 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         <  2:2172958/93‑1 (MQ=255)

           cAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGc                                                                        >  1:326797/1‑93 (MQ=255)

           cAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGc                                                                        >  2:1310632/1‑93 (MQ=255)

           cAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGc                                                                        <  2:2807204/93‑1 (MQ=255)

            cgCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                                       <  2:727401/92‑1 (MQ=255)

            aGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                                       >  2:4531984/1‑93 (MQ=255)

             gCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCgg                                                                      >  2:2712501/1‑93 (MQ=255)

               gTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgag                                                                    <  2:5095099/93‑1 (MQ=255)

                tCTGCTGCACACATCCCTAGTAAATAGGGTCGGGGGCCGCCGCTTCGAGCGCCGGGTTGCGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                                   <  1:9692/93‑1 (MQ=255)

                tCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                                   >  2:4049343/1‑93 (MQ=255)

                tCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                                   >  1:5306283/1‑93 (MQ=255)

                tCTGCTCCACAAACCCCCAGAAAATA‑CGGGGAGGCCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagag                                                                  <  1:3519715/93‑1 (MQ=255)

                 ctgctgCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagag                                                                  <  2:47535/93‑1 (MQ=255)

                  tgctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                                 <  1:3480193/93‑1 (MQ=255)

                  tgctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                                 <  2:3333662/93‑1 (MQ=255)

                  tgctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCACGCTACGCATTAAATGAACAgagaga                                                                 >  1:2904200/1‑93 (MQ=255)

                   gctgcACACACCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                <  1:1617043/93‑1 (MQ=255)

                   gctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                <  2:6581571/93‑1 (MQ=255)

                   gctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                >  1:6171493/1‑93 (MQ=255)

                    ctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtc                                                               >  2:6039510/1‑93 (MQ=255)

                    ctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtc                                                               <  1:6041181/93‑1 (MQ=255)

                    ctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGCGAtc                                                               >  1:3991397/1‑93 (MQ=255)

                     tgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtct                                                              <  1:5048903/93‑1 (MQ=255)

                      ccacacaACCCCAGAAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtctt                                                             <  1:27002/92‑1 (MQ=255)

                        ccacaCCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttct                                                           <  1:1766417/92‑1 (MQ=255)

                        acacaACCCCAGTAAATATCGTCGAGGGCTGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttct                                                           <  1:1051201/93‑1 (MQ=255)

                        acacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttct                                                           >  2:3980904/1‑93 (MQ=255)

                          acaacaCCAGTAAATATCGCCGAGGGCCGCCGCTTCCAGCGCCGAGTTACGCGCAAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgc                                                         <  1:6387569/93‑1 (MQ=255)

                          acaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgc                                                         <  1:364857/93‑1 (MQ=255)

                            aaCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcac                                                       >  2:178960/1‑93 (MQ=255)

                              ccccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                                     >  1:6112010/1‑93 (MQ=255)

                               cccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgca                                                    <  2:4427528/93‑1 (MQ=255)

                                  aGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGt                                                 >  2:5439582/1‑93 (MQ=255)

                                      aaTATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                                             >  1:1230205/1‑93 (MQ=255)

                                        tatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                                           <  1:1815797/93‑1 (MQ=255)

                                        tatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                                           <  2:4796759/93‑1 (MQ=255)

                                        tatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                                           >  2:5053285/1‑93 (MQ=255)

                                           cgtcgGGGGCGGCCTCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                                        <  1:5604494/93‑1 (MQ=255)

                                           cgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                                        >  2:5642980/1‑93 (MQ=255)

                                           cgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                                        >  1:5332275/1‑93 (MQ=255)

                                             ccgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCa                                      <  1:3015338/92‑1 (MQ=255)

                                              cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCCCGCACGTTACGAAAAGCAc                                     <  1:3185501/93‑1 (MQ=255)

                                              cgAGGGCCGCCGCTTCCAGCACCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                                     >  1:6456659/1‑93 (MQ=255)

                                                aGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                                   >  2:1582854/1‑93 (MQ=255)

                                                 gggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGc                                  >  2:222206/1‑93 (MQ=255)

                                                 gggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGc                                  >  1:4236968/1‑93 (MQ=255)

                                                 gggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGc                                  >  1:4027338/1‑93 (MQ=255)

                                                  ggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGcc                                 >  1:6308142/1‑93 (MQ=255)

                                                    ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt                               <  2:1230953/93‑1 (MQ=255)

                                                    ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt                               <  1:6588875/93‑1 (MQ=255)

                                                      gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc                             >  2:2408352/1‑93 (MQ=255)

                                                      gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc                             <  2:3085528/93‑1 (MQ=255)

                                                      gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc                             >  1:6508302/1‑93 (MQ=255)

                                                       ccgcTTCCAGCGCCGGTTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg                            <  1:2522398/93‑1 (MQ=255)

                                                        cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa                           >  1:314493/1‑93 (MQ=255)

                                                         gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa                          <  1:4642301/93‑1 (MQ=255)

                                                         gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa                          <  1:2119069/93‑1 (MQ=255)

                                                          cTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAAc                         >  2:4342618/1‑93 (MQ=255)

                                                           ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg                        >  1:4394653/1‑93 (MQ=255)

                                                            tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACgg                       >  2:4190263/1‑93 (MQ=255)

                                                            tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACgg                       >  1:3837600/1‑93 (MQ=255)

                                                              cAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGGACCGcc                     >  1:4147917/1‑93 (MQ=255)

                                                              cAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGcc                     >  2:2514720/1‑93 (MQ=255)

                                                               aGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGccc                    <  2:1952514/93‑1 (MQ=255)

                                                                gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa                   <  1:1545667/93‑1 (MQ=255)

                                                                gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa                   >  1:909443/1‑93 (MQ=255)

                                                                 cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCAt                  <  1:6720628/93‑1 (MQ=255)

                                                                  gcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATc                 <  2:3837604/93‑1 (MQ=255)

                                                                   ccGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCg                >  2:1675346/1‑93 (MQ=255)

                                                                   ccGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCg                >  1:5848068/1‑93 (MQ=255)

                                                                    cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgg               >  2:3645057/1‑93 (MQ=255)

                                                                    cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgg               >  2:2957750/1‑93 (MQ=255)

                                                                     gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCggg              >  2:259235/1‑93 (MQ=255)

                                                                      ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgggg             >  2:2530968/1‑93 (MQ=255)

                                                                      ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgggg             <  1:1000684/93‑1 (MQ=255)

                                                                       gTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGt            >  2:6289226/1‑93 (MQ=255)

                                                                        ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTg           <  1:662976/93‑1 (MQ=255)

                                                                           cgcgCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGa        >  1:2786845/1‑93 (MQ=255)

                                                                             cgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAAt      >  2:5180695/1‑93 (MQ=255)

                                                                             cgcAAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAAt      <  1:3448658/93‑1 (MQ=255)

                                                                              gcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAAtt     <  1:4511795/93‑1 (MQ=255)

                                                                              gcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAAtt     <  2:5374952/93‑1 (MQ=255)

                                                                               ccAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTg    >  2:3914518/1‑93 (MQ=255)

                                                                                 agcaggCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCa  <  1:797934/93‑1 (MQ=255)

||

AAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCA  >  NC_000913/4026166‑4026339

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