BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	541,786	+G	coding (675/1146 nt)	ybbZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	541,782	1	.	G	97.2%	330.5 / 2.8	108	coding (671/1146 nt)	ybbZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/3); new base G (49/56); total (49/59)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 124 total aligned reads displayed.

CGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑G‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTA  >  NC_000913/541690‑541868

cGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑a                                                                                          <  2:3846106/93‑1 (MQ=255)

cGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑a                                                                                          <  1:4558008/93‑1 (MQ=255)

 gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑g                                                                                        <  1:2819103/93‑1 (MQ=255)

 gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGGCGAAGAACTTGAAC‑A‑g                                                                                        >  1:301218/1‑93 (MQ=255)

  ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑g‑g                                                                                      >  1:3484160/1‑93 (MQ=255)

  ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑g‑g                                                                                      >  2:1607941/1‑93 (MQ=255)

    gATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑g‑ggg                                                                                    >  2:3199736/1‑93 (MQ=255)

       aaCCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑g‑ggggct                                                                                 >  2:451756/1‑90 (MQ=255)

          ccGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATTGTCGAAGAACTTGAAC‑AGG‑GGGCTGCa                                                                               >  2:1686054/1‑93 (MQ=255)

             cTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCaaca                                                                            >  2:1510036/1‑91 (MQ=255)

             cTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCaaaa                                                                            <  1:3902606/93‑1 (MQ=255)

              tGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAAt                                                                           >  1:3313696/1‑93 (MQ=255)

               ggTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAAtt                                                                          <  1:4635549/93‑1 (MQ=255)

                gTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTa                                                                         >  2:3018945/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGc                                                                      >  2:1057357/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGc                                                                      >  2:3536759/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGc                                                                      <  1:2036102/93‑1 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGc                                                                      >  2:167077/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGc                                                                      >  1:1101264/1‑93 (MQ=255)

                     ggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGccc                                                                    <  2:1619163/93‑1 (MQ=255)

                      gCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCg                                                                   <  1:2855040/93‑1 (MQ=255)

                        aCGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCgtg                                                                 <  2:2846402/93‑1 (MQ=255)

                         cgcggcgcggCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCgtgt                                                                <  2:2612027/93‑1 (MQ=255)

                          gcggcgcggcTGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGtt                                                               <  2:3836247/93‑1 (MQ=255)

                          gcggcgcggcGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGtt                                                               <  1:4671533/93‑1 (MQ=255)

                           ggacgcggcgGCGGTGTTTGGCCGGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGACCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGttt                                                              <  2:3954479/90‑1 (MQ=255)

                            ggcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGtttt                                                             >  2:3608832/1‑93 (MQ=255)

                             gcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTa                                                            <  1:286354/93‑1 (MQ=255)

                             gcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTa                                                            <  2:3777765/93‑1 (MQ=255)

                              cgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTATGCCCGTGTTTTac                                                           <  2:4501927/93‑1 (MQ=255)

                              cgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTac                                                           >  2:2572151/1‑93 (MQ=255)

                               gcggcggcgGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTaca                                                          <  1:3839535/93‑1 (MQ=255)

                                cggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaa                                                         <  2:450145/93‑1 (MQ=255)

                                cggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaa                                                         >  2:328551/1‑93 (MQ=255)

                                cggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaa                                                         <  1:14614/93‑1 (MQ=255)

                                 ggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaac                                                        >  1:178404/1‑93 (MQ=255)

                                  gcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaaca                                                       <  2:369362/93‑1 (MQ=255)

                                  gcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTacaaca                                                       <  2:1994687/93‑1 (MQ=255)

                                   gggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAg                                                      <  1:3420734/92‑1 (MQ=255)

                                     gcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCa                                                    >  2:4011809/1‑93 (MQ=255)

                                     gcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCa                                                    >  2:4135096/1‑93 (MQ=255)

                                      cggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCaa                                                   >  2:2358702/1‑93 (MQ=255)

                                       ggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCaaa                                                  <  2:543959/93‑1 (MQ=255)

                                       ggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCaaa                                                  <  2:4012644/93‑1 (MQ=255)

                                       ggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCaaa                                                  >  2:1494925/1‑93 (MQ=255)

                                        gtgtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                 <  1:2050812/93‑1 (MQ=255)

                                        gtgtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                 >  2:506174/1‑93 (MQ=255)

                                        gtgtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                 >  2:385734/1‑93 (MQ=255)

                                          gtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTg                                               <  2:1661056/93‑1 (MQ=255)

                                          gtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTg                                               >  2:1706086/1‑93 (MQ=255)

                                            ttGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGaa                                             >  2:327131/1‑93 (MQ=255)

                                             tGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGaaa                                            >  1:3736404/1‑93 (MQ=255)

                                              ggCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAAt                                           >  2:3185867/1‑93 (MQ=255)

                                              ggCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAAt                                           <  2:1162074/93‑1 (MQ=255)

                                                 ccGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTaa                                        <  1:3179514/93‑1 (MQ=255)

                                                  cGCAAAAAGGGGCAACGCCGGAGATTGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAAt                                       <  1:1686054/93‑1 (MQ=255)

                                                  cGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAAt                                       <  2:1773479/93‑1 (MQ=255)

                                                  cGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAAt                                       <  1:2689175/93‑1 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACGG‑G‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      <  1:4140119/93‑1 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      <  1:3527529/93‑1 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      >  2:3053099/1‑93 (MQ=255)

                                                    gaaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGt                                     <  1:4072562/92‑1 (MQ=255)

                                                    cAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGt                                     <  2:1602145/93‑1 (MQ=255)

                                                    cAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGt                                     <  2:4063832/93‑1 (MQ=255)

                                                     aaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTc                                    <  1:3259154/93‑1 (MQ=255)

                                                      aaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCt                                   >  2:1290875/1‑93 (MQ=255)

                                                      aaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCt                                   <  1:3921438/93‑1 (MQ=255)

                                                      aaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCt                                   >  1:4532729/1‑93 (MQ=255)

                                                      aaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGCGTTTTACAACAGCAAACTGAAATTAATGTCt                                   >  1:111776/1‑93 (MQ=255)

                                                       aaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTg                                  <  2:2889651/93‑1 (MQ=255)

                                                       aaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑GCGGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTg                                  <  1:798283/93‑1 (MQ=255)

                                                        aaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGc                                 >  1:4362134/1‑93 (MQ=255)

                                                         aGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGcc                                >  2:270957/1‑93 (MQ=255)

                                                          ggggCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCa                               <  1:355511/93‑1 (MQ=255)

                                                          ggggCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCa                               >  1:3521903/1‑93 (MQ=255)

                                                          ggggCAACACCGGAGATGATCGACGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAACCTGAAATTAATGTCTGCCa                               <  2:2378475/93‑1 (MQ=255)

                                                             gCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGAt                            <  1:1466744/93‑1 (MQ=255)

                                                             gCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGAt                            <  1:653935/93‑1 (MQ=255)

                                                             gCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGAt                            <  1:609772/93‑1 (MQ=255)

                                                             gCAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGAt                            <  1:1752955/93‑1 (MQ=255)

                                                              cAACGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATg                           <  2:1923563/93‑1 (MQ=255)

                                                                aCGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggc                         <  2:4629740/93‑1 (MQ=255)

                                                                aCGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggc                         >  1:3601608/1‑93 (MQ=255)

                                                                 cGCCGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcg                        <  1:2067409/93‑1 (MQ=255)

                                                                    cGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcgggc                     >  2:1812535/1‑93 (MQ=255)

                                                                    cGGAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcgggc                     >  2:1695791/1‑93 (MQ=255)

                                                                     ggAGATGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcgggcg                    >  1:4164346/1‑93 (MQ=255)

                                                                       agaTGGTCGAAGAACTGGAAC‑AGG‑GGGCTGGCAAAGTACCCCCGTGTTTTACAAAAGCAAACTTAAATTAATGTCTTCCAGATTGCGggcggc                  <  2:644847/93‑1 (MQ=255)

                                                                        gaTGGTCGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGggcggcg                 >  2:406732/1‑93 (MQ=255)

                                                                              cGAAGAACTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTg           >  2:1818425/1‑93 (MQ=255)

                                                                               gaagaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCCGGCGGCGGCGCTgc          <  2:294373/93‑1 (MQ=255)

                                                                                aagaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcg         >  2:3922334/1‑93 (MQ=255)

                                                                                aagaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcg         <  2:3595025/93‑1 (MQ=255)

                                                                                aagaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcg         >  1:125550/1‑93 (MQ=255)

                                                                                 agaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgg        >  2:3893272/1‑93 (MQ=255)

                                                                                  gaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcggg       <  1:3000298/93‑1 (MQ=255)

                                                                                  gaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcggg       >  1:2953926/1‑93 (MQ=255)

                                                                                  gaaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcggg       >  1:2478944/1‑93 (MQ=255)

                                                                                   aaCTTGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgggc      <  1:3645585/93‑1 (MQ=255)

                                                                                       tGAAC‑AGG‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  <  1:1989976/93‑1 (MQ=255)

CGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAAC‑A‑G‑GGGCTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTA  >  NC_000913/541690‑541868

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