BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,632	T→C	Q84R (CAG→CGG)	fecB ←	FecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,632	0	T	C	99.5%	613.1 / ‑5.6	200	Q84R (CAG→CGG)	fecB	FecB
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (81/118); total (82/118)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 216 total aligned reads displayed.

CCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCG  >  NC_000913/4511542‑4511721

cccGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGc                                                                                         >  1:421936/1‑93 (MQ=255)

cccGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGc                                                                                         >  2:2208895/1‑93 (MQ=255)

cccGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACCGc                                                                                         >  1:2431811/1‑93 (MQ=255)

  cGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCa                                                                                       <  1:2803124/93‑1 (MQ=255)

   gCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCAc                                                                                      >  1:2904065/1‑93 (MQ=255)

     aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCATGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCGCgg                                                                                    >  2:1890602/1‑93 (MQ=255)

     aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACgg                                                                                    >  2:1164588/1‑93 (MQ=255)

     aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACgg                                                                                    >  2:1707469/1‑93 (MQ=255)

      tGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGt                                                                                   >  1:3968254/1‑93 (MQ=255)

       ggCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGtt                                                                                  <  1:3820524/93‑1 (MQ=255)

          gACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCa                                                                               <  2:2326910/93‑1 (MQ=255)

          gACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCa                                                                               <  2:821864/93‑1 (MQ=255)

          gACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCa                                                                               <  1:2177619/93‑1 (MQ=255)

          gACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCa                                                                               <  1:2479719/93‑1 (MQ=255)

           aCTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAg                                                                              <  2:458224/93‑1 (MQ=255)

            ctgctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAgg                                                                             <  1:2581152/93‑1 (MQ=255)

            ctgctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGCTTCAgg                                                                             >  2:2396968/1‑93 (MQ=255)

             tgctgTCGGCAATGATCGGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGt                                                                            >  2:3039891/1‑93 (MQ=255)

               ctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTgc                                                                          <  2:1858747/93‑1 (MQ=255)

                 gTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGGTTCAGGTgggc                                                                        >  2:3955824/1‑93 (MQ=255)

                   cGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGGCCGCCACGGTTTCTGGTGcgcac                                                                      >  1:1589521/1‑93 (MQ=255)

                     gCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGcgcacgc                                                                    >  2:1357405/1‑93 (MQ=255)

                      cAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGcgcacgca                                                                   <  1:4431636/93‑1 (MQ=255)

                       aaTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCgcacgcac                                                                  >  2:4003495/1‑93 (MQ=255)

                       aaTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCgcacgcac                                                                  >  1:1068847/1‑93 (MQ=255)

                        aTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACt                                                                 >  1:3140231/1‑93 (MQ=255)

                        aTGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACt                                                                 >  2:2928412/1‑93 (MQ=255)

                         tGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACtt                                                                >  2:3286279/1‑93 (MQ=255)

                          gATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTc                                                               <  2:2088542/93‑1 (MQ=255)

                            tCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCgg                                                             >  2:1140454/1‑93 (MQ=255)

                            tCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCgg                                                             >  1:1265412/1‑93 (MQ=255)

                            tCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCgg                                                             <  2:3622581/93‑1 (MQ=255)

                            tCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCCCACGCACTTCgg                                                             >  2:570296/1‑93 (MQ=255)

                               ggTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCa                                                          <  1:1983468/93‑1 (MQ=255)

                               ggTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGGTTCAGGTGCGCACGCACTTCGGGCa                                                          >  1:102349/1‑93 (MQ=255)

                                gTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAg                                                         <  1:3842337/93‑1 (MQ=255)

                                 tCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAgg                                                        >  2:4171361/1‑93 (MQ=255)

                                   tGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCCGGAt                                                      >  1:4126624/1‑93 (MQ=255)

                                   tGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGAt                                                      <  2:1407920/93‑1 (MQ=255)

                                   tGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGAt                                                      >  2:2625475/1‑93 (MQ=255)

                                    ggTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATg                                                     <  2:4087983/93‑1 (MQ=255)

                                     gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGc                                                    <  2:3129570/93‑1 (MQ=255)

                                        tCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtt                                                 <  1:2497049/93‑1 (MQ=255)

                                         cAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCTGGCAGGATGGGttt                                                >  1:2901018/1‑93 (MQ=255)

                                          agagCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACTGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtttt                                               <  1:1245564/93‑1 (MQ=255)

                                             gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                            <  1:34807/93‑1 (MQ=255)

                                             gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                            <  2:2948765/93‑1 (MQ=255)

                                              cGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCAt                                           <  2:4479207/93‑1 (MQ=255)

                                               ggCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATc                                          <  2:515995/93‑1 (MQ=255)

                                                gCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCg                                         >  2:3207665/1‑93 (MQ=255)

                                                 cAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGt                                        <  1:3567161/93‑1 (MQ=255)

                                                  aaTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGtt                                       <  1:3383003/93‑1 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                      >  1:2130406/1‑93 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                      >  1:2875695/1‑93 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                      <  2:4377069/93‑1 (MQ=255)

                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTa                                      >  1:2303753/1‑93 (MQ=255)

                                                    tGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAt                                     <  2:3731883/93‑1 (MQ=255)

                                                    tGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAt                                     <  1:3599903/93‑1 (MQ=255)

                                                     ggCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                    <  1:148075/93‑1 (MQ=255)

                                                     ggCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                    <  1:2718444/93‑1 (MQ=255)

                                                      gCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGGTTTAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcg                                   >  2:1855593/1‑93 (MQ=255)

                                                       cTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgt                                  <  2:1235824/93‑1 (MQ=255)

                                                       cTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgt                                  <  2:3883786/93‑1 (MQ=255)

                                                       cTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgt                                  <  1:4296807/93‑1 (MQ=255)

                                                        ttCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtc                                 <  2:3946862/93‑1 (MQ=255)

                                                        ttCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtc                                 >  2:2427240/1‑93 (MQ=255)

                                                         tCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgtcg                                >  1:693620/1‑93 (MQ=255)

                                                           cAGGCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGc                              >  1:930630/1‑93 (MQ=255)

                                                             ggCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCaa                            >  2:3810844/1‑93 (MQ=255)

                                                             ggCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCaa                            <  1:1766951/93‑1 (MQ=255)

                                                              gCTCGGCTGCGCGCGCGTTCCGACGGACCTCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAAt                           <  1:402544/93‑1 (MQ=255)

                                                              gCTCGGCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAAt                           <  1:2646430/93‑1 (MQ=255)

                                                                   gCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGa                      <  1:4578166/93‑1 (MQ=255)

                                                                   gCTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGa                      <  1:3785814/93‑1 (MQ=255)

                                                                    cTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGAt                     <  2:1975568/93‑1 (MQ=255)

                                                                    cTGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGAt                     <  2:1639121/93‑1 (MQ=255)

                                                                    cTGCGCGCGCGGTCCGATGGACCGCCACGGTTTCAGGTGCTCACGCACTTCGGGCAGGATGCGTTGTGTATCGTTATCGTCGGCAATACCGAt                     <  2:3481276/93‑1 (MQ=255)

                                                                     tGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATAGCGAtt                    >  2:3243662/1‑92 (MQ=255)

                                                                     tGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATc                    >  1:875922/1‑93 (MQ=255)

                                                                     tGCGCGCGCGTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATc                    >  1:3102648/1‑93 (MQ=255)

                                                                     cgcgcgcgcGTTCCGACGCACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATc                    <  2:2447140/92‑1 (MQ=255)

                                                                       cgcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCgg                  >  1:4502395/1‑93 (MQ=255)

                                                                       cgcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCgg                  <  1:3197851/93‑1 (MQ=255)

                                                                        gcgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCggg                 <  1:1376268/93‑1 (MQ=255)

                                                                         cgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGa                <  2:883302/93‑1 (MQ=255)

                                                                         cgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGa                <  1:72921/93‑1 (MQ=255)

                                                                         cgcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGa                >  1:275989/1‑93 (MQ=255)

                                                                          gcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGAt               <  2:3421069/93‑1 (MQ=255)

                                                                          gcgcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGAt               <  2:500785/93‑1 (MQ=255)

                                                                            gcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGa             <  2:2721128/93‑1 (MQ=255)

                                                                            gcgTTCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGa             <  2:723504/93‑1 (MQ=255)

                                                                               ttCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGt          <  1:1029743/93‑1 (MQ=255)

                                                                                tCCGACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTc         <  1:625008/93‑1 (MQ=255)

                                                                                 ccGACGGACCGCCACGGTTTCAGGTTCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTcc        >  2:2540453/1‑93 (MQ=255)

                                                                                 ccGACGGACCGCCACGGTTTCAGGTGCGCACGCACTCCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTcc        >  2:3526145/1‑93 (MQ=255)

                                                                                   gACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAc      <  1:555050/93‑1 (MQ=255)

                                                                                   gACGGACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAc      <  1:184145/93‑1 (MQ=255)

                                                                                      ggACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcggc   >  1:496365/1‑93 (MQ=255)

                                                                                      ggACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcggc   <  2:1454500/93‑1 (MQ=255)

                                                                                       gACCGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAcggcg  >  2:136442/1‑93 (MQ=255)

CCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCG  >  NC_000913/4511542‑4511721

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