BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,451	T→C	A279A (GCA→GCG)	uidB ←	glucuronide permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,451	0	T	C	89.7%	479.6 / 42.0	195	A279A (GCA→GCG)	uidB	glucuronide permease
Reads supporting (aligned to +/- strand): ref base T (7/13); new base C (97/78); total (104/91)
Fisher's exact test for biased strand distribution p-value = 9.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-02

Only 100 of 204 total aligned reads displayed.

GATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGG  >  NC_000913/1691363‑1691542

aaaaaCCGCAGGTCCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCACCCCCGGCACCAGCGGTGCCGACGCca                                                                                         <  1:7383808/90‑1 (MQ=255)

   aaCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAg                                                                                      <  1:6566806/93‑1 (MQ=255)

   aaCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAg                                                                                      <  2:3300788/93‑1 (MQ=255)

    aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGt                                                                                     <  2:3637996/93‑1 (MQ=255)

     ccGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTa                                                                                    <  1:6013933/93‑1 (MQ=255)

      cGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTac                                                                                   >  2:1259713/1‑93 (MQ=255)

      cGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACCGTac                                                                                   >  2:4801110/1‑93 (MQ=255)

       gCAGGTTCCCAGGAAAGGCCAAATCAGGAAGGTATTCTTTTTACCGATACTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTacc                                                                                  <  1:3286203/93‑1 (MQ=255)

       gCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTacc                                                                                  >  2:1698734/1‑93 (MQ=255)

        cAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTacca                                                                                 >  1:9443385/1‑93 (MQ=255)

         aGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTaccaa                                                                                <  2:3427057/93‑1 (MQ=255)

           gTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaacc                                                                              <  1:274250/93‑1 (MQ=255)

           gTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTaccaacc                                                                              >  2:6513286/1‑93 (MQ=255)

             tCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAg                                                                            >  2:8350096/1‑93 (MQ=255)

             tCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAg                                                                            >  1:9961982/1‑93 (MQ=255)

               ccAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGt                                                                          >  1:252411/1‑93 (MQ=255)

               ccAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGt                                                                          <  2:6834282/93‑1 (MQ=255)

                 aGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTAACAACCAGGttt                                                                        >  2:4175783/1‑93 (MQ=255)

                  gCAAAGCGCCAATCGGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCACGGGACCCGCGACGACGACGCCACAGTACCAAACAGGtttt                                                                       >  2:6897863/1‑93 (MQ=255)

                  gCAAAGCGCCAATCAGGAAGCTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGtttt                                                                       >  1:1163973/1‑93 (MQ=255)

                   cAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTg                                                                      >  1:314050/1‑93 (MQ=255)

                   cAAAGCGCCAAGCAGGAAGGTATTTTTTTTACTAATCGTCGCGACCACGCCCGGCCCGATCGTTGCCGACGCCACAGTACCAACCAGGTTTTg                                                                      <  2:5014247/93‑1 (MQ=255)

                    aaaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGGTTTGc                                                                     >  1:1971938/1‑93 (MQ=255)

                    aaaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGCTTTGc                                                                     >  1:7642481/1‑93 (MQ=255)

                     aaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCa                                                                    >  1:5020630/1‑93 (MQ=255)

                     aaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCa                                                                    >  2:9434629/1‑93 (MQ=255)

                       gccccAATCAGGAAGTTATTCTTTTTCCCGACCCCCGCGACCACCCCCGGCACCGGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCAcc                                                                  <  2:261061/90‑1 (MQ=255)

                          ccAATCAGTAAGTTATTTTTTTTCCCGATCCCCGCGACCATCCCCGGCCCCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGt                                                               <  2:7986616/93‑1 (MQ=255)

                          ccAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGt                                                               <  2:6107509/93‑1 (MQ=255)

                          ccAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCCGCCCCAGCGGGGCCGTCGCCACGGTACCAACCAGGTTTTCCACCAGt                                                               >  1:7577575/1‑93 (MQ=255)

                            aaTCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGGCGACGCCACAGTACCAACCAGGCTTTGCACCAGTAc                                                             >  1:958949/1‑93 (MQ=255)

                            aaTCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTAc                                                             >  2:5086826/1‑93 (MQ=255)

                              tCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCa                                                           >  1:1186331/1‑93 (MQ=255)

                              tCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCa                                                           >  1:8564036/1‑93 (MQ=255)

                              tCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCa                                                           <  1:641050/93‑1 (MQ=255)

                              tCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCa                                                           >  2:8876167/1‑93 (MQ=255)

                                aGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCaccagcaccagc                                                         >  2:4121214/1‑93 (MQ=255)

                                aGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGc                                                         >  2:3161316/1‑93 (MQ=255)

                                  gAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGcac                                                       >  2:671032/1‑93 (MQ=255)

                                    aGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCCCCACCCGTGCCGACGCCCCAGTACCGACCAGGTTTTGCAACACTACCACCACAg                                                     >  1:3878472/1‑93 (MQ=255)

                                    aGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAg                                                     <  1:6296270/93‑1 (MQ=255)

                                     ggTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGt                                                    >  2:7558972/1‑93 (MQ=255)

                                     ggTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTCGCACCAGTACCAGCACAGt                                                    >  1:4149238/1‑93 (MQ=255)

                                      gTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                   >  1:7168702/1‑93 (MQ=255)

                                      gTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                   >  2:9514324/1‑93 (MQ=255)

                                      gTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGGGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                   >  1:8209988/1‑93 (MQ=255)

                                       tATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGa                                                  <  2:2101391/93‑1 (MQ=255)

                                       tATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGa                                                  >  2:1951885/1‑93 (MQ=255)

                                        aTTCTTTTTACCGACCCTCGCGACCACCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGaa                                                 <  1:7863906/93‑1 (MQ=255)

                                         ttCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                >  1:4629905/1‑93 (MQ=255)

                                         ttCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                <  2:8521603/93‑1 (MQ=255)

                                         ttCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                <  1:1299843/93‑1 (MQ=255)

                                         ttCTTTTTACCGACCCTCGCGACCATCCCCGGCCCCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                <  2:8015897/93‑1 (MQ=255)

                                            tttttACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGGTTTGCACCAGTACCAGCACAGTGAACaac                                             >  2:8274180/1‑91 (MQ=255)

                                                tACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCgg                                         <  2:9358366/93‑1 (MQ=255)

                                                 aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGt                                        >  1:5687624/1‑93 (MQ=255)

                                                 aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGt                                        >  1:8041933/1‑93 (MQ=255)

                                                 aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGCGAACAGCCCgga                                        >  2:9061746/1‑92 (MQ=255)

                                                  ccGATCCTCGCGACCATCCCCGGCACCAGCGGGGCCGACGCCACAGTACCAACCAGGGTTTGCACCAGTACCAGCACAGTGAACAGCCCggga                                       >  1:9031905/1‑91 (MQ=255)

                                                   cGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGGACAATGAACAACCCGGTAt                                      >  2:8084235/1‑93 (MQ=255)

                                                   cGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTAt                                      <  2:4406262/93‑1 (MQ=255)

                                                   cGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTAt                                      >  2:9625119/1‑93 (MQ=255)

                                                   cGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTAt                                      >  2:2192872/1‑93 (MQ=255)

                                                      tCCTCGCGCCCATCCCCGGCCCCAGCGGTCCCGACGCCCCAGTACCAGCCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                   <  1:4955737/93‑1 (MQ=255)

                                                      tCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                   >  1:7773469/1‑93 (MQ=255)

                                                      tCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                   <  2:5364013/93‑1 (MQ=255)

                                                      tCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                   >  1:3778008/1‑93 (MQ=255)

                                                      tCCTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                   >  2:8100902/1‑93 (MQ=255)

                                                        cTCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAttt                                 >  2:3799908/1‑93 (MQ=255)

                                                        cTCGCGACCACCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACATCCCGGTATCAttt                                 <  2:3203247/93‑1 (MQ=255)

                                                         tCGCGACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTa                                >  1:8559402/1‑93 (MQ=255)

                                                         tCGCGACCATCCCCGGCACCAGCGGGGCCGACGCCACAGAACCAACCAGGTTTGGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATCTa                                >  2:8972427/1‑93 (MQ=255)

                                                          cgcgACCATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTaa                               >  2:7578744/1‑93 (MQ=255)

                                                               ccATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACAt                          >  2:8798597/1‑93 (MQ=255)

                                                                cccccccGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATa                         <  1:1084479/90‑1 (MQ=255)

                                                                cATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCTCAGTGAACAGCACGGTATCATTTAACACATa                         >  2:3365124/1‑93 (MQ=255)

                                                                cATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATa                         <  2:2817240/93‑1 (MQ=255)

                                                                cATCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATa                         >  1:6328913/1‑93 (MQ=255)

                                                                 aTCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGACCAGCCCGGTATCATTTAACACATAg                        >  1:413348/1‑93 (MQ=255)

                                                                 aTCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAg                        >  1:36208/1‑93 (MQ=255)

                                                                 aTCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAg                        >  2:7069753/1‑93 (MQ=255)

                                                                  tCCCCGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAgc                       >  1:8560753/1‑93 (MQ=255)

                                                                    cccGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAgcgc                     <  2:1276615/93‑1 (MQ=255)

                                                                     ccGGCACCAGTGGTGCCGACGCCACAGTACCAACCAGGGTTTGCACCAGTACCAGCACAGTGAACAGCCCGGGATCATTTAACACATAGCGCa                    >  2:9239861/1‑93 (MQ=255)

                                                                     ccGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCa                    <  2:9777371/93‑1 (MQ=255)

                                                                      cGGCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCAc                   >  2:5838015/1‑93 (MQ=255)

                                                                       ggCACCAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACg                  <  1:8768176/93‑1 (MQ=255)

                                                                           ccAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGa              <  1:8148504/93‑1 (MQ=255)

                                                                           ccAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGa              <  2:7880996/93‑1 (MQ=255)

                                                                           ccAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGa              >  2:6044522/1‑93 (MQ=255)

                                                                            cAGCGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaa             <  1:3251308/93‑1 (MQ=255)

                                                                            cAGCGGTGCCGACGCCACAGTACCAACCAGGTTTCGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaa             <  1:1013533/93‑1 (MQ=255)

                                                                             aGCGGGGCCGACGCCACGGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaac            <  1:7718889/93‑1 (MQ=255)

                                                                               cGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaacaa          >  1:2056323/1‑93 (MQ=255)

                                                                               cGGTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaacaa          <  1:1574628/93‑1 (MQ=255)

                                                                                ggTCCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaacaac         <  2:408287/93‑1 (MQ=255)

                                                                                 gTGCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacg        <  1:6165969/93‑1 (MQ=255)

                                                                                   gCCGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacgac      <  1:356611/93‑1 (MQ=255)

                                                                                     cGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAAcgacga    <  1:529100/93‑1 (MQ=255)

                                                                                     cGACGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAAcgacga    <  2:4041538/93‑1 (MQ=255)

GATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGG  >  NC_000913/1691363‑1691542

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