BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,654,259	T→C	T38A (ACC→GCC)	hdeA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,654,259	0	T	C	83.1%	652.7 / 130.0	301	T38A (ACC→GCC)	hdeA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (20/31); new base C (121/129); total (141/160)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 303 total aligned reads displayed.

TCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATA  >  NC_000913/3654168‑3654350

  ttctGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCa                                                                                          <  1:7044352/93‑1 (MQ=255)

  ttctGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCa                                                                                          <  2:3911290/93‑1 (MQ=255)

   tctGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAg                                                                                         <  2:1898690/93‑1 (MQ=255)

   tctGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAg                                                                                         <  2:6832348/93‑1 (MQ=255)

     tGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAggc                                                                                       <  2:7405617/93‑2 (MQ=255)

     tGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGa                                                                                       <  1:3005995/93‑1 (MQ=255)

     tGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGa                                                                                       >  1:7163488/1‑93 (MQ=255)

      ggTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAg                                                                                      <  1:1856471/93‑1 (MQ=255)

      ggTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAg                                                                                      <  2:6589537/93‑1 (MQ=255)

       gTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGt                                                                                     <  2:7117362/93‑1 (MQ=255)

       gTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGt                                                                                     <  1:3340217/93‑1 (MQ=255)

         ttATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTg                                                                                   <  2:1447217/93‑1 (MQ=255)

          tATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGa                                                                                  <  2:6464208/93‑1 (MQ=255)

           aTCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGAc                                                                                 <  1:1038705/93‑1 (MQ=255)

           aTCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGAc                                                                                 <  1:1589214/93‑1 (MQ=255)

           aTCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGAc                                                                                 <  1:1441787/93‑1 (MQ=255)

            tCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGAcc                                                                                <  1:2972293/93‑1 (MQ=255)

             cTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCg                                                                               <  2:5008589/93‑1 (MQ=255)

               ttgttgTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGt                                                                             <  2:3626737/93‑1 (MQ=255)

               ttgttgTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGt                                                                             >  1:1258524/1‑93 (MQ=255)

                 gttgttCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGttt                                                                           <  1:3974358/93‑1 (MQ=255)

                  ttgttCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGtttt                                                                          >  1:3381449/1‑93 (MQ=255)

                   tgttCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGttttt                                                                         <  1:6610267/93‑1 (MQ=255)

                   tgttCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGttttt                                                                         <  1:940014/93‑1 (MQ=255)

                      tCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGt                                                                      >  1:2855435/1‑93 (MQ=255)

                       cAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGtt                                                                     <  1:3068498/93‑1 (MQ=255)

                           gcTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCa                                                                 <  2:4586042/93‑1 (MQ=255)

                           gcTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCa                                                                 <  2:5844270/93‑1 (MQ=255)

                            cTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAg                                                                >  2:3739537/1‑93 (MQ=255)

                            cTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAg                                                                >  1:5077975/1‑93 (MQ=255)

                             ttCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAgc                                                               <  2:2834020/93‑1 (MQ=255)

                             ttCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAgc                                                               >  1:1198788/1‑93 (MQ=255)

                             ttCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAgc                                                               <  1:1420978/93‑1 (MQ=255)

                              tCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAgct                                                              >  1:24098/1‑93 (MQ=255)

                              tCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAgct                                                              <  2:5613581/93‑1 (MQ=255)

                              tCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAgct                                                              <  1:1823294/93‑1 (MQ=255)

                                aGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAgctgc                                                            <  1:3469883/93‑1 (MQ=255)

                                aGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAgctgc                                                            <  1:6848855/93‑1 (MQ=255)

                                 gCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAgctgct                                                           >  2:5459284/1‑93 (MQ=255)

                                  cAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCtt                                                          >  2:4109822/1‑93 (MQ=255)

                                   aaaaCCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCttt                                                         <  1:5235445/93‑1 (MQ=255)

                                       ccAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTgc                                                     <  2:2294283/93‑1 (MQ=255)

                                        cAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTgcg                                                    >  1:2111959/1‑93 (MQ=255)

                                        cAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTgcg                                                    >  1:1648606/1‑93 (MQ=255)

                                          aCTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCa                                                  >  1:735483/1‑93 (MQ=255)

                                           cTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCNt                                                 >  2:7172436/1‑93 (MQ=255)

                                           cTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCAt                                                 >  1:4456171/1‑93 (MQ=255)

                                            tGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATc                                                <  2:2657574/93‑1 (MQ=255)

                                             gCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATcc                                               <  1:3866551/93‑1 (MQ=255)

                                               aGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGc                                             <  1:4571742/93‑1 (MQ=255)

                                                gTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCt                                            >  2:136114/1‑93 (MQ=255)

                                                gTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCt                                            >  1:538411/1‑93 (MQ=255)

                                                gTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCt                                            <  2:1998173/93‑1 (MQ=255)

                                                gTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCt                                            >  1:6302958/1‑93 (MQ=255)

                                                 ttGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTg                                           >  1:5803889/1‑93 (MQ=255)

                                                  tGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGc                                          <  2:2228498/93‑1 (MQ=255)

                                                  tGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGc                                          >  2:6186927/1‑93 (MQ=255)

                                                   ggCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCa                                         <  1:4290344/93‑1 (MQ=255)

                                                     cTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAtt                                       <  2:7024145/93‑1 (MQ=255)

                                                     cTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAtt                                       <  2:4766973/93‑1 (MQ=255)

                                                     cTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCAtt                                       <  2:5812451/93‑1 (MQ=255)

                                                       ggAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGc                                     >  1:2635388/1‑93 (MQ=255)

                                                        gAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCt                                    >  2:3824771/1‑93 (MQ=255)

                                                          aGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTca                                  <  2:1619183/93‑1 (MQ=255)

                                                           ggATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTcac                                 >  2:4490924/1‑93 (MQ=255)

                                                            gATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTcaca                                <  2:4961099/93‑1 (MQ=255)

                                                            gATTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTcaca                                >  1:5995902/1‑93 (MQ=255)

                                                             aTTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACaa                               <  1:3566661/93‑1 (MQ=255)

                                                             aTTCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACaa                               >  2:1167626/1‑93 (MQ=255)

                                                              ttCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAAc                              <  2:4159275/93‑1 (MQ=255)

                                                               tCGTCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACt                             <  2:1767557/93‑1 (MQ=255)

                                                                  tCCACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGc                          >  2:506151/1‑93 (MQ=255)

                                                                   ccACAGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCa                         >  1:2581430/1‑93 (MQ=255)

                                                                    cacaGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCag                        >  1:2284991/1‑93 (MQ=255)

                                                                    cacaGCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCag                        >  1:3146691/1‑93 (MQ=255)

                                                                       aGCCAGGAAATCTTCACAGGTCCAGGAGTTGACAGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaag                     <  1:461850/93‑1 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    <  2:2559022/93‑1 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    >  1:6180734/1‑93 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    >  1:5676462/1‑93 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    >  1:5266273/1‑93 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    <  2:5350966/93‑1 (MQ=255)

                                                                        gCCAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCagaaga                    <  2:5579196/93‑1 (MQ=255)

                                                                          cAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGaagaag                  <  2:6676628/93‑1 (MQ=255)

                                                                          cAGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGaagaag                  >  2:2461210/1‑93 (MQ=255)

                                                                           aGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGc                 <  2:3520385/93‑1 (MQ=255)

                                                                           aGGAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGc                 >  2:3241476/1‑93 (MQ=255)

                                                                             gAAATCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAg               >  2:789046/1‑93 (MQ=255)

                                                                              aaaTCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGa              <  2:3115732/93‑1 (MQ=255)

                                                                              aaaTCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGa              <  1:2488988/93‑1 (MQ=255)

                                                                              aaaTCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGa              <  1:4509468/93‑1 (MQ=255)

                                                                               aaTCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGac             <  2:3995840/93‑1 (MQ=255)

                                                                                aTCTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGacc            >  1:4826137/1‑93 (MQ=255)

                                                                                  cTTCACAGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGaccac          <  2:4230012/93‑1 (MQ=255)

                                                                                      acaGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAg      >  1:4173665/1‑93 (MQ=255)

                                                                                       caGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGa     >  1:7378878/1‑93 (MQ=255)

                                                                                        aGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGaa    <  2:2393634/93‑1 (MQ=255)

                                                                                        aGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGaa    <  2:2059504/93‑1 (MQ=255)

                                                                                        aGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGaa    >  1:3699958/1‑93 (MQ=255)

                                                                                        aGGCCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGaa    <  1:6475940/93‑1 (MQ=255)

                                                                                          gTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATa  <  2:1399788/93‑1 (MQ=255)

TCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATA  >  NC_000913/3654168‑3654350

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