BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,412	T→C	E82E (GAA→GAG)	relE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,412	0	T	C	98.5%	690.0 / ‑1.3	197	E82E (GAA→GAG)	relE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (84/110); total (85/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 206 total aligned reads displayed.

AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAT  >  NC_000913/1643322‑1643501

aaGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCtc                                                                                         <  1:876567/93‑1 (MQ=255)

aaGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCtc                                                                                         <  1:704600/93‑1 (MQ=255)

  gCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCtctc                                                                                       >  2:332514/1‑93 (MQ=255)

  gCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCtctc                                                                                       <  1:3432072/93‑1 (MQ=255)

  gCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTGTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCtctc                                                                                       <  1:31490/93‑1 (MQ=255)

   cAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGctctct                                                                                      <  2:5847960/93‑1 (MQ=255)

    aGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCtt                                                                                     <  2:6021201/93‑1 (MQ=255)

    aGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCtt                                                                                     <  2:1615965/93‑1 (MQ=255)

    aGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCtt                                                                                     <  1:4204968/93‑1 (MQ=255)

     gAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCttt                                                                                    >  1:5352900/1‑93 (MQ=255)

       aGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTc                                                                                  <  1:5169594/93‑1 (MQ=255)

        gtgtCGCCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTcc                                                                                 <  1:3578279/93‑1 (MQ=255)

         tgtCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTccc                                                                                >  2:2234604/1‑93 (MQ=255)

         tgtCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTccc                                                                                <  2:2407447/93‑1 (MQ=255)

         tgtCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTccc                                                                                >  2:1643242/1‑93 (MQ=255)

           tCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCaa                                                                              >  2:246995/1‑93 (MQ=255)

            cACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAAc                                                                             <  2:3793729/93‑1 (MQ=255)

            cACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAAc                                                                             <  2:227455/93‑1 (MQ=255)

             cccTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACa                                                                            <  1:2448991/92‑1 (MQ=255)

              ccTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAg                                                                           <  1:3471432/93‑1 (MQ=255)

               cTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGa                                                                          <  2:1706226/93‑1 (MQ=255)

               cTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGa                                                                          >  1:3616455/1‑93 (MQ=255)

                 tCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGaaa                                                                        <  1:1429598/93‑1 (MQ=255)

                 tCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGaaa                                                                        >  2:4759516/1‑93 (MQ=255)

                  cGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAAt                                                                       <  1:2261967/93‑1 (MQ=255)

                   ggTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATc                                                                      <  1:1335566/93‑1 (MQ=255)

                    gTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCa                                                                     >  1:5780455/1‑93 (MQ=255)

                        gAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCTAACAGAAATCACGaa                                                                 <  2:5995966/93‑1 (MQ=255)

                           aCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAAc                                                              >  1:5912079/1‑93 (MQ=255)

                           aCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAAc                                                              <  2:3662305/93‑1 (MQ=255)

                           aCAGAGATGTCATGCTTTGGCTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAAc                                                              >  1:5620476/1‑93 (MQ=255)

                            cAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACg                                                             <  1:2634113/93‑1 (MQ=255)

                             agagATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGa                                                            >  1:4465906/1‑93 (MQ=255)

                             agagATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGa                                                            >  2:6067857/1‑93 (MQ=255)

                              gagaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGAc                                                           >  2:2832568/1‑93 (MQ=255)

                                gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                         <  1:1052566/93‑1 (MQ=255)

                                gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                         >  2:2564131/1‑93 (MQ=255)

                                gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                         <  2:3951443/93‑1 (MQ=255)

                                gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                         >  1:1412667/1‑93 (MQ=255)

                                 aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                        <  2:4404561/93‑1 (MQ=255)

                                   gTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACtt                                                      >  2:1215496/1‑93 (MQ=255)

                                    tCATGCTTTGGTTCAGAGAATGCGTTTGCCCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                     <  2:309432/93‑1 (MQ=255)

                                    tCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                     >  1:4040540/1‑93 (MQ=255)

                                     cATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTtc                                                    <  2:5511422/93‑1 (MQ=255)

                                      aTGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTtct                                                   <  1:3553611/93‑1 (MQ=255)

                                         cTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGt                                                <  1:37573/93‑1 (MQ=255)

                                         cTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGt                                                <  1:1696219/93‑1 (MQ=255)

                                          tttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTc                                               >  2:1488945/1‑93 (MQ=255)

                                           ttGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCt                                              <  2:4127818/93‑1 (MQ=255)

                                             ggTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtat                                            <  2:4977933/93‑1 (MQ=255)

                                              gTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtata                                           <  2:6113388/93‑1 (MQ=255)

                                              gTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtata                                           <  2:720881/93‑1 (MQ=255)

                                              gTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCtata                                           <  1:2089309/93‑1 (MQ=255)

                                               ttCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATaa                                          <  2:3441443/93‑1 (MQ=255)

                                               ttCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATaa                                          <  2:458326/93‑1 (MQ=255)

                                                 cAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAAcc                                        <  1:1269127/93‑1 (MQ=255)

                                                   gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTg                                      >  2:3960711/1‑93 (MQ=255)

                                                   gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTg                                      <  1:556705/93‑1 (MQ=255)

                                                   gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTg                                      >  1:899174/1‑93 (MQ=255)

                                                     gaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGt                                    <  2:3549950/93‑1 (MQ=255)

                                                      aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGta                                   <  1:278304/93‑1 (MQ=255)

                                                      aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGta                                   >  2:3033235/1‑93 (MQ=255)

                                                      aaTGCGTTTGACCGCCTCGCTATATACTCCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGta                                   <  1:3129498/93‑1 (MQ=255)

                                                       aTGCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtat                                  <  1:4760133/93‑1 (MQ=255)

                                                         gCGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATAc                                <  1:4894416/93‑1 (MQ=255)

                                                          cGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACa                               <  2:5273065/93‑1 (MQ=255)

                                                          cGTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACa                               <  1:3515284/93‑1 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              >  2:5140531/1‑93 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              >  1:3194979/1‑93 (MQ=255)

                                                           gTTTGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACaa                              <  2:424039/93‑1 (MQ=255)

                                                            tttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAg                             >  1:2919871/1‑93 (MQ=255)

                                                            gttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAg                             <  1:3223162/92‑1 (MQ=255)

                                                             ttGACCGCCTCGCTATATTCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAgg                            <  2:1021511/93‑1 (MQ=255)

                                                             ttGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAgg                            <  1:6069881/93‑1 (MQ=255)

                                                              tGACCGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGc                           <  1:3859419/93‑1 (MQ=255)

                                                                 ccGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGTCAGCTTTCTCGTCTATAACCTGGTATACAAGGCGAt                        <  1:4164096/93‑1 (MQ=255)

                                                                 ccGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGAt                        <  1:946489/93‑1 (MQ=255)

                                                                 ccGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGAt                        <  2:2098403/93‑1 (MQ=255)

                                                                 ccGCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGAt                        <  2:2584426/93‑1 (MQ=255)

                                                                   gCCTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAg                      >  1:2015600/1‑93 (MQ=255)

                                                                     cTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGcc                    <  2:5773305/93‑1 (MQ=255)

                                                                     cTCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGcc                    <  2:5609915/93‑1 (MQ=255)

                                                                      tCGCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCt                   >  1:5536663/1‑93 (MQ=255)

                                                                        gCTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGa                 >  1:3376039/1‑93 (MQ=255)

                                                                         cTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGaa                >  1:5191464/1‑93 (MQ=255)

                                                                         cTATATACTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGaa                <  2:2948531/93‑1 (MQ=255)

                                                                          tatataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAATGCGATAGCCTGAAg               >  1:5109895/1‑93 (MQ=255)

                                                                           atataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGa              >  2:3804631/1‑93 (MQ=255)

                                                                            tataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGAc             <  2:5187534/93‑1 (MQ=255)

                                                                            tataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGAc             >  1:86164/1‑93 (MQ=255)

                                                                            tataCTTCCGAGCGCTCTCTTTTCACAACAGACATCACGACAACGTCAACTTTCTCGTCTATAACCTGGTATACAAGGCGATACACTGAAGAc             >  1:2316883/1‑93 (MQ=255)

                                                                             ataCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGAcc            <  2:4850819/93‑1 (MQ=255)

                                                                             ataCTCCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACATCTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGAcc            <  1:4612430/93‑1 (MQ=255)

                                                                              taCTTCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATATCCTGGTATACAAGGCGATAGCCTGAAGACCg           <  2:4008198/93‑1 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        <  2:3843674/93‑1 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        <  2:573354/93‑1 (MQ=255)

                                                                                 ttCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg        <  2:4328220/93‑1 (MQ=255)

                                                                                  tCCGAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGc       <  2:4468136/93‑1 (MQ=255)

                                                                                     gAGCGCTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTa    >  2:5355030/1‑93 (MQ=255)

                                                                                       gcgcTCTCTTTTCCCAACAGAAATCAGGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGTGATAGCCTGAAGACCGGAGCTTAAt  >  2:1445895/1‑93 (MQ=255)

AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAAT  >  NC_000913/1643322‑1643501

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