BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	680,175	A→G	E247E (GAA→GAG)	ybeV →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	680,175	0	A	G	86.2%	718.1 / 97.2	275	E247E (GAA→GAG)	ybeV	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (25/13); new base G (105/132); total (130/145)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

Only 100 of 286 total aligned reads displayed.

AAGTTACATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGCGCGTACT  >  NC_000913/680083‑680266

aaGTTACATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGAAACTCATTCCTGTCGCTCAACAGTGGCAAACGgaagaa                                                                                             >  1:1040180/1‑93 (MQ=255)

  gTTACATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGcc                                                                                           >  2:2505692/1‑93 (MQ=255)

   ttACATCGCAACCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCa                                                                                          <  2:5913705/93‑1 (MQ=255)

   ttACATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGCACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCa                                                                                          >  2:4854922/1‑93 (MQ=255)

     aCATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGCAGAGCCAGa                                                                                        >  1:4707648/1‑93 (MQ=255)

      cATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGaa                                                                                       >  2:3754605/1‑93 (MQ=255)

       aTCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGaaa                                                                                      <  1:3258475/93‑1 (MQ=255)

       aTCACAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGaaa                                                                                      <  1:5114288/93‑1 (MQ=255)

        tCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAg                                                                                     >  2:2751505/1‑93 (MQ=255)

        tCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAg                                                                                     <  1:5164119/93‑1 (MQ=255)

        tCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAg                                                                                     >  2:6056701/1‑93 (MQ=255)

         cGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGc                                                                                    <  2:5446701/93‑1 (MQ=255)

           cAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGa                                                                                  <  1:3552582/93‑1 (MQ=255)

             aaCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAg                                                                                <  2:1132139/93‑1 (MQ=255)

             aaCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAg                                                                                <  1:6461087/93‑1 (MQ=255)

             aaCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAg                                                                                >  1:461178/1‑93 (MQ=255)

             aaCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAg                                                                                >  2:2660404/1‑93 (MQ=255)

             aaCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAg                                                                                >  2:1358891/1‑93 (MQ=255)

              aCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGa                                                                               <  2:5772796/93‑1 (MQ=255)

              aCTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGa                                                                               >  2:6361876/1‑93 (MQ=255)

                tGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATg                                                                             <  1:5410226/93‑1 (MQ=255)

                tGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATg                                                                             <  1:4778638/93‑1 (MQ=255)

                tGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATg                                                                             >  2:5956963/1‑93 (MQ=255)

                 gTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGc                                                                            <  1:6644460/93‑1 (MQ=255)

                   tGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGccc                                                                          <  2:6602388/93‑1 (MQ=255)

                    gcgcTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCg                                                                         >  1:3229477/1‑93 (MQ=255)

                    gcgcTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCg                                                                         <  2:693533/93‑1 (MQ=255)

                      gcTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCgga                                                                       >  1:3060212/1‑92 (MQ=255)

                       cTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTa                                                                      <  2:6004249/93‑1 (MQ=255)

                       cTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTa                                                                      <  1:2182250/93‑1 (MQ=255)

                           tACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTAc                                                                  >  2:4815578/1‑93 (MQ=255)

                             cAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACtt                                                                <  1:6191957/93‑1 (MQ=255)

                             cAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACtt                                                                <  1:1916552/93‑1 (MQ=255)

                               gCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTtgt                                                              <  2:2816203/93‑1 (MQ=255)

                               gCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTtgt                                                              >  2:655584/1‑93 (MQ=255)

                               gCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTtgt                                                              >  1:4461668/1‑93 (MQ=255)

                                 tCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTtgtgt                                                            <  2:77688/93‑1 (MQ=255)

                                   ggTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGc                                                          >  2:5284964/1‑93 (MQ=255)

                                      gCAAACAGGTATGGCGGGACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGcaca                                                       >  1:5872705/1‑93 (MQ=255)

                                       cAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACaa                                                      >  2:2557119/1‑93 (MQ=255)

                                          aCAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGt                                                   <  2:3248332/93‑1 (MQ=255)

                                          aCAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGt                                                   >  1:3985302/1‑93 (MQ=255)

                                            aGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGttt                                                 >  2:1296863/1‑93 (MQ=255)

                                            aGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGttt                                                 >  2:3897839/1‑93 (MQ=255)

                                             ggTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTg                                                <  1:5549923/93‑1 (MQ=255)

                                             ggTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTg                                                <  2:74363/93‑1 (MQ=255)

                                              gTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGt                                               <  1:1998054/93‑1 (MQ=255)

                                               tATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTa                                              >  1:3366628/1‑93 (MQ=255)

                                                aTGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTAc                                             <  1:4215904/93‑1 (MQ=255)

                                                aTGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTAc                                             >  2:121376/1‑93 (MQ=255)

                                                  ggCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTg                                           <  2:4759254/93‑1 (MQ=255)

                                                  ggCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTg                                           >  2:2607844/1‑93 (MQ=255)

                                                   gCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgc                                          <  2:936269/93‑1 (MQ=255)

                                                   gCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgc                                          <  1:2590747/93‑1 (MQ=255)

                                                   gCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgc                                          <  1:3719039/93‑1 (MQ=255)

                                                   gCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgc                                          <  2:1839272/93‑1 (MQ=255)

                                                    cGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgcg                                         >  2:3102335/1‑93 (MQ=255)

                                                     ggAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgcgg                                        <  2:34376/93‑1 (MQ=255)

                                                     ggAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgcgg                                        <  1:781974/93‑1 (MQ=255)

                                                       aaCTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTgcggcg                                      >  2:4018922/1‑93 (MQ=255)

                                                        aCTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGa                                     >  1:407123/1‑93 (MQ=255)

                                                        aCTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGa                                     <  2:2885274/93‑1 (MQ=255)

                                                          tCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAgg                                   <  1:6315628/93‑1 (MQ=255)

                                                          tCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAgg                                   >  1:5468600/1‑93 (MQ=255)

                                                           cATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggg                                  <  1:3688438/93‑1 (MQ=255)

                                                           cATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggg                                  <  1:4501403/93‑1 (MQ=255)

                                                           cATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggg                                  >  2:2366528/1‑93 (MQ=255)

                                                           cATTCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggg                                  >  2:4232791/1‑93 (MQ=255)

                                                            aTTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAgggg                                 >  1:5741748/1‑93 (MQ=255)

                                                             ttCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGgggggg                                >  1:2444554/1‑93 (MQ=255)

                                                             ttCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggggg                                >  2:2424480/1‑93 (MQ=255)

                                                             ttCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAggggg                                >  2:5632339/1‑93 (MQ=255)

                                                              tCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGa                               <  1:3467448/93‑1 (MQ=255)

                                                              tCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGa                               >  1:2405031/1‑93 (MQ=255)

                                                              tCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGa                               <  1:2827525/93‑1 (MQ=255)

                                                              tCCTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGa                               <  1:3780063/93‑1 (MQ=255)

                                                                cTGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGaaa                             <  2:2053277/93‑1 (MQ=255)

                                                                 tGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAg                            <  1:4390369/93‑1 (MQ=255)

                                                                 tGTCGCTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAg                            >  1:1861036/1‑93 (MQ=255)

                                                                      cTCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGc                       <  2:4603396/93‑1 (MQ=255)

                                                                       tCAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCt                      >  2:6271743/1‑93 (MQ=255)

                                                                       tCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCt                      >  1:1677042/1‑93 (MQ=255)

                                                                       tCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGGATTACTTGCGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCt                      >  2:5708882/1‑93 (MQ=255)

                                                                        cAACAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCtt                     <  2:3931927/93‑1 (MQ=255)

                                                                          aCAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGc                   <  1:732021/93‑1 (MQ=255)

                                                                           cAGTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGcc                  <  2:6157733/93‑1 (MQ=255)

                                                                             gTGGCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGa                <  1:3765970/93‑1 (MQ=255)

                                                                               ggCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGAtc              <  1:4306897/93‑1 (MQ=255)

                                                                               ggCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGAtc              <  2:5813911/93‑1 (MQ=255)

                                                                               ggCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGAtc              <  2:3630955/93‑1 (MQ=255)

                                                                                gCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGAtct             <  1:1123105/93‑1 (MQ=255)

                                                                                gCAAACGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGAtct             <  2:5145421/93‑1 (MQ=255)

                                                                                  aaaCGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATctct           <  1:5971877/93‑1 (MQ=255)

                                                                                   aaCGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTg          <  2:4461220/93‑1 (MQ=255)

                                                                                   aaCGGAAGAGCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTg          >  1:1352541/1‑93 (MQ=255)

                                                                                       gaagagCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGcgcg      >  2:3784289/1‑93 (MQ=255)

                                                                                        aagagCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGCGCGt     >  1:981148/1‑93 (MQ=255)

                                                                                         agagCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGCGCGTa    >  2:2230816/1‑93 (MQ=255)

                                                                                         agagCCAGAAAGCGAAGATGCCCGGAATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGCAAGCCTGCTTGCCGATCTCTGCGCGTa    >  1:1273116/1‑93 (MQ=255)

                                                                                           aaCCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGCGCGTACt  >  2:3339377/1‑93 (MQ=255)

AAGTTACATCGCAAACTGTTGCGCTGGTACAGCTCGGTGCAAACAGGTATGGCGGAACTCATTCCTGTCGCTCAACAGTGGCAAACGGAAGAACCAGAAAGCGAAGATGCCCGGTATTACTTGTGTGCACAACGTTTGTACTGCGGCGAGGGGGAAAGCCTGCTTGCCGATCTCTGCGCGTACT  >  NC_000913/680083‑680266

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