BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,934,461	A→G	P186P (CCA→CCG)	rbsB →	D‑ribose periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,934,461	0	A	G	100.0%	806.6 / NA	246	P186P (CCA→CCG)	rbsB	D‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (110/136); total (110/136)

Only 100 of 283 total aligned reads displayed.

CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550

 aTTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  2:6395128/93‑1 (MQ=255)

  ttGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggcag                                                                                       >  2:1397602/1‑93 (MQ=255)

   tGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGa                                                                                      <  1:335977/93‑1 (MQ=255)

    gCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAt                                                                                     >  1:4204556/1‑93 (MQ=255)

      tGGTATATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAttt                                                                                   <  2:4871214/93‑1 (MQ=255)

      tGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAttt                                                                                   <  1:6299509/93‑1 (MQ=255)

      tGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAttt                                                                                   >  2:4085373/1‑93 (MQ=255)

       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:719651/93‑1 (MQ=255)

       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  1:3493224/93‑1 (MQ=255)

       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:3076637/93‑1 (MQ=255)

       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:4418070/93‑1 (MQ=255)

        gTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTg                                                                                 >  1:2911193/1‑93 (MQ=255)

          aCATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAt                                                                               >  2:4390415/1‑93 (MQ=255)

           cATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATc                                                                              <  2:1751355/93‑1 (MQ=255)

           cATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATc                                                                              >  1:1794306/1‑93 (MQ=255)

               cGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAt                                                                          <  1:848211/93‑1 (MQ=255)

                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         <  1:2000378/93‑1 (MQ=255)

                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         >  1:6451134/1‑93 (MQ=255)

                  aGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaa                                                                       >  1:2443036/1‑93 (MQ=255)

                  aGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaa                                                                       <  2:4643893/93‑1 (MQ=255)

                  aGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaa                                                                       >  1:6697052/1‑93 (MQ=255)

                   gCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaaa                                                                      <  2:2254113/93‑1 (MQ=255)

                     ccGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAgg                                                                    >  2:1213828/1‑93 (MQ=255)

                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  2:3700519/93‑1 (MQ=255)

                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  2:5579662/93‑1 (MQ=255)

                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  1:5969687/93‑1 (MQ=255)

                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  <  1:3554601/93‑1 (MQ=255)

                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  <  1:1277318/93‑1 (MQ=255)

                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 >  1:3102174/1‑93 (MQ=255)

                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 <  1:2309810/93‑1 (MQ=255)

                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 <  1:1184483/93‑1 (MQ=255)

                         gAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTg                                                                <  1:4219556/93‑1 (MQ=255)

                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              >  1:5919007/1‑93 (MQ=255)

                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              >  2:3921623/1‑93 (MQ=255)

                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              >  2:4144272/1‑93 (MQ=255)

                           aCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              >  1:1062838/1‑93 (MQ=255)

                           aCGGGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGaa                                                              <  1:1494527/93‑1 (MQ=255)

                            cGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAAc                                                             <  2:3236587/93‑1 (MQ=255)

                            cGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAAc                                                             <  2:976279/93‑1 (MQ=255)

                            cGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAAc                                                             <  1:2706428/93‑1 (MQ=255)

                            cGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAAc                                                             <  1:497383/93‑1 (MQ=255)

                              tGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGt                                                           <  2:6155754/93‑1 (MQ=255)

                              tGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGt                                                           <  1:4693940/93‑1 (MQ=255)

                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  1:2164158/93‑1 (MQ=255)

                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  1:2108726/93‑1 (MQ=255)

                                 cGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAAt                                                        <  1:1903130/93‑1 (MQ=255)

                                 cGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAAt                                                        <  2:630919/93‑1 (MQ=255)

                                  gAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATg                                                       >  2:5829903/1‑93 (MQ=255)

                                  gAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATg                                                       >  2:1672029/1‑93 (MQ=255)

                                    aGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTCATGCa                                                     >  1:3205286/1‑93 (MQ=255)

                                    aGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCa                                                     >  2:1623317/1‑93 (MQ=255)

                                     ggCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGACAGCCAGCCGGCAGATTTTGATCGCATTAAAAGTTTGCACGTAATGCAg                                                    >  1:1016902/1‑93 (MQ=255)

                                      gCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGa                                                   <  2:2544786/93‑1 (MQ=255)

                                         tCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAcc                                                <  2:3991792/93‑1 (MQ=255)

                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCt                                               <  1:1885394/93‑1 (MQ=255)

                                           cagcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTg                                              >  1:6080485/1‑93 (MQ=255)

                                            agcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGt                                             <  1:154982/93‑1 (MQ=255)

                                            agcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGt                                             <  2:3566956/93‑1 (MQ=255)

                                            agcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGt                                             >  1:5768498/1‑93 (MQ=255)

                                            agcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGt                                             >  1:383295/1‑93 (MQ=255)

                                             gcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGtt                                            >  1:6194637/1‑93 (MQ=255)

                                             gcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGCTTGAACGTAATGCAGAACCTGtt                                            >  1:1439927/1‑93 (MQ=255)

                                              cagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTACAGGTTTGAACGTAATGCAGAACCTGTTg                                           >  2:6486281/1‑93 (MQ=255)

                                               agGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          >  2:806124/1‑93 (MQ=255)

                                                 gCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAcc                                        <  2:3319234/93‑1 (MQ=255)

                                                  ccGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCg                                       >  2:4281366/1‑93 (MQ=255)

                                                   cGTTGCTGCTCACAAGTTTCATGTTCTTGCCTGCCAACCGACAGATTTTGCTCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGc                                      <  2:6722992/93‑1 (MQ=255)

                                                   cGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGc                                      <  2:1096952/93‑1 (MQ=255)

                                                   cGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGc                                      <  1:895666/93‑1 (MQ=255)

                                                    gTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCt                                     >  1:5265348/1‑93 (MQ=255)

                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  2:2774609/1‑93 (MQ=255)

                                                      tgctgcTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCa                                   <  2:5346485/93‑1 (MQ=255)

                                                        ctgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  2:5962355/93‑1 (MQ=255)

                                                          tctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAATGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCg                               <  1:3206739/92‑1 (MQ=255)

                                                             cacaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAt                            <  1:2840691/93‑1 (MQ=255)

                                                             cacaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAt                            <  2:3419858/93‑1 (MQ=255)

                                                             cacaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAt                            >  1:3958845/1‑93 (MQ=255)

                                                               caAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATgc                          <  2:2039926/93‑2 (MQ=255)

                                                               caAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATgc                          <  2:2396507/93‑2 (MQ=255)

                                                               caAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATgc                          >  1:1615767/1‑92 (MQ=255)

                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       >  2:5762189/1‑93 (MQ=255)

                                                                   tttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAg                      <  1:5471092/93‑1 (MQ=255)

                                                                    ttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAgg                     >  2:2947652/1‑93 (MQ=255)

                                                                     tattGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  1:5062031/90‑1 (MQ=255)

                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  1:3894052/93‑1 (MQ=255)

                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  2:2109355/93‑1 (MQ=255)

                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  2:5127088/93‑1 (MQ=255)

                                                                      aaTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCt                   <  2:3482189/93‑1 (MQ=255)

                                                                       aTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTg                  <  2:4844537/93‑1 (MQ=255)

                                                                       aTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTg                  >  2:3274995/1‑93 (MQ=255)

                                                                          ttCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTAt               <  1:950328/93‑1 (MQ=255)

                                                                             ttGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcg            >  2:1373907/1‑93 (MQ=255)

                                                                             ttGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcg            <  2:382400/93‑1 (MQ=255)

                                                                                ccagccagccggCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCgcgc         >  1:4081986/1‑93 (MQ=255)

                                                                                  agccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAg       >  2:5025969/1‑93 (MQ=255)

                                                                                   gccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGa      >  1:4757675/1‑93 (MQ=255)

                                                                                      agccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatg   >  1:5008711/1‑93 (MQ=255)

                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  >  2:5633559/1‑93 (MQ=255)

                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  >  2:2482671/1‑93 (MQ=255)

                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  <  1:4287874/93‑1 (MQ=255)

CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550

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