BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	628,816	A→G	V98V (GTA→GTG)	ybdB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,816	0	A	G	97.0%	497.4 / 0.9	164	V98V (GTA→GTG)	ybdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (90/69); total (93/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 176 total aligned reads displayed.

GCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628727‑628905

  cGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcgg                                                                                      >  2:2835721/1‑93 (MQ=255)

   ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     <  1:1364585/93‑1 (MQ=255)

   ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     <  1:3923135/93‑1 (MQ=255)

      tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  >  1:3007391/1‑93 (MQ=255)

       ttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCt                                                                                 >  1:2381842/1‑93 (MQ=255)

         atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               >  1:1797926/1‑93 (MQ=255)

          tgatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              <  2:585753/93‑1 (MQ=255)

           gatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGGGCGCGGCGTCTGCCa                                                                             >  1:3402716/1‑93 (MQ=255)

            atgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAg                                                                            >  2:2042308/1‑93 (MQ=255)

             tgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGc                                                                           <  1:2900660/93‑1 (MQ=255)

              gaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGcc                                                                          <  2:271132/93‑1 (MQ=255)

               aCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCg                                                                         >  1:2883714/1‑93 (MQ=255)

                cccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGc                                                                        >  1:2516844/1‑93 (MQ=255)

                 ccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCt                                                                       >  2:1667480/1‑93 (MQ=255)

                   gcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGc                                                                     >  2:935658/1‑93 (MQ=255)

                   gcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGc                                                                     >  1:2296373/1‑93 (MQ=255)

                     gACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCAt                                                                   >  2:4544446/1‑93 (MQ=255)

                      aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATc                                                                  >  1:593316/1‑93 (MQ=255)

                      aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGGACGCGGCGTCTGCCAGCCGCTGCATc                                                                  >  1:2996746/1‑93 (MQ=255)

                       cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  2:1114267/93‑1 (MQ=255)

                       cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  2:3046329/93‑1 (MQ=255)

                       cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  1:1216443/93‑1 (MQ=255)

                        ggACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCAAGCCGCTGCATCtt                                                                >  2:1326650/1‑93 (MQ=255)

                         gACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTg                                                               >  2:816548/1‑93 (MQ=255)

                          aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              <  1:1086390/93‑1 (MQ=255)

                          aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              <  2:4371171/93‑1 (MQ=255)

                          aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              >  1:2227329/1‑93 (MQ=255)

                          aCAGTGTGTGGTAGGCACAGAACTTAAGGCAACAATCCAACGCCGGGTGTCTGAGGGAAGGGTGCCCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              <  2:215374/93‑1 (MQ=255)

                                    gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    <  1:3283560/93‑1 (MQ=255)

                                    gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  2:3453874/1‑93 (MQ=255)

                                      aGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                  >  2:3785238/1‑93 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  1:1383504/93‑1 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  2:1100686/93‑1 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  1:3125094/93‑1 (MQ=255)

                                        gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  2:3456529/93‑1 (MQ=255)

                                         cacaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               <  1:1886715/93‑1 (MQ=255)

                                         cacaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               >  2:3505482/1‑93 (MQ=255)

                                         cacaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               <  1:2149202/93‑1 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagagt                                              >  1:3812345/1‑92 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  2:3363368/1‑93 (MQ=255)

                                          acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  1:3085217/1‑93 (MQ=255)

                                           caGAACTTAATGCAACCCACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCt                                             <  1:962783/93‑1 (MQ=255)

                                           caGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCt                                             <  2:4466352/93‑1 (MQ=255)

                                           caGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGAGCGCGGCGTCTGCCAGCCGCTGCATCTGGGTCGGCAAAATCAGAGCt                                             >  2:4225277/1‑93 (MQ=255)

                                            aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            <  1:2671038/93‑1 (MQ=255)

                                            aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            <  2:4414843/93‑1 (MQ=255)

                                            aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            >  2:1929105/1‑93 (MQ=255)

                                            aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            <  1:4679198/93‑1 (MQ=255)

                                             gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           >  1:461890/1‑93 (MQ=255)

                                             gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           >  1:4275999/1‑93 (MQ=255)

                                              aaCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTggg                                          >  2:3992083/1‑93 (MQ=255)

                                               aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         >  1:1573678/1‑93 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  1:1955172/1‑93 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       <  1:1677667/93‑1 (MQ=255)

                                                 ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  1:4611901/1‑93 (MQ=255)

                                                   aaTGCAACACACCATCGCCTGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtc                                     <  2:3700244/93‑1 (MQ=255)

                                                   aaTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtc                                     <  2:4568470/93‑1 (MQ=255)

                                                    aTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcg                                    <  2:1552799/93‑1 (MQ=255)

                                                      tcAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  <  2:294585/92‑1 (MQ=255)

                                                      gCAACACACCAGCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  <  2:2449585/93‑1 (MQ=255)

                                                        aaCACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                >  2:4290082/1‑93 (MQ=255)

                                                        aaCACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                <  1:921000/93‑1 (MQ=255)

                                                          cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  2:360690/93‑1 (MQ=255)

                                                          cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              >  1:994207/1‑93 (MQ=255)

                                                           ccacCATCGCGCGGTGTCTAAGGGAAAGGTGCGGGGCGTCTTCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtttt                             <  1:155283/92‑1 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  1:395488/1‑93 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  2:3717141/1‑93 (MQ=255)

                                                             acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  2:766917/1‑93 (MQ=255)

                                                               cATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         <  2:4535772/93‑1 (MQ=255)

                                                               cATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         <  1:929092/93‑1 (MQ=255)

                                                                aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        <  2:4341042/93‑1 (MQ=255)

                                                                aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        <  1:1832534/93‑1 (MQ=255)

                                                                 tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGGGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       >  2:3936136/1‑93 (MQ=255)

                                                                   gCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAAc                     >  2:486362/1‑93 (MQ=255)

                                                                    cccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACa                    >  1:2392698/1‑93 (MQ=255)

                                                                     ccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAg                   <  1:4083496/93‑1 (MQ=255)

                                                                     ccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAg                   >  1:2252575/1‑93 (MQ=255)

                                                                     ccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTAGTTTTCGATGAACAg                   >  1:2906685/1‑93 (MQ=255)

                                                                      cGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  >  2:1064704/1‑93 (MQ=255)

                                                                      cGGTGTCTGAGAGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  <  2:1496404/93‑1 (MQ=255)

                                                                      cGGCGTCTGAGGGAAAGGTGCGCGGCGTCTTCCAGCAGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  >  2:3302808/1‑93 (MQ=255)

                                                                       ggggtCTGAGGGAGGGGGGCGCGCGGTCTGCCTGCCGCTGCATCTTGGGCGACAAATTCGGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 <  2:2605256/90‑1 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 >  2:3372453/1‑93 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 <  1:2955681/93‑1 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 >  1:1728619/1‑93 (MQ=255)

                                                                       ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 >  2:3804296/1‑93 (MQ=255)

                                                                        gtgtCTGGGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:3466018/93‑1 (MQ=255)

                                                                        gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                >  1:2524260/1‑93 (MQ=255)

                                                                        gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:50373/93‑1 (MQ=255)

                                                                         tgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggc               <  1:738653/93‑1 (MQ=255)

                                                                          gtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg              >  2:2606317/1‑93 (MQ=255)

                                                                           tCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcgg             >  2:2222617/1‑93 (MQ=255)

                                                                           tCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcgg             >  2:476484/1‑93 (MQ=255)

                                                                            cTGAGGGAAAGGTGCGCGGCGTCGGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcggc            <  1:1412307/93‑1 (MQ=255)

                                                                                gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtg        >  1:4245175/1‑93 (MQ=255)

                                                                                gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtg        >  1:3438497/1‑93 (MQ=255)

                                                                                gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtg        >  1:4528471/1‑93 (MQ=255)

                                                                                 ggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgc       <  1:4387939/93‑1 (MQ=255)

                                                                                  gAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgct      >  1:1488402/1‑93 (MQ=255)

                                                                                   aaaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctg     >  1:2864993/1‑93 (MQ=255)

GCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628727‑628905

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