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	99.4%	533.2 / ‑4.5	166	Q224R (CAG→CGG)	ybbZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); major base G (71/93); minor base T (1/0); total (73/93)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 197 total aligned reads displayed.

CGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑GGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTA  >  NC_000913/541690‑541868

cGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGCAGAACTTGAACa                                                                                          <  2:5043444/93‑1 (MQ=255)

cGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACa                                                                                          <  1:1225500/93‑1 (MQ=255)

cGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACa                                                                                          <  1:2260071/93‑1 (MQ=255)

 gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑g                                                                                        <  1:1827388/93‑1 (MQ=255)

 gTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑g                                                                                        >  1:2906245/1‑93 (MQ=255)

  ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGCACTTGAACA‑gg                                                                                       >  2:284060/1‑93 (MQ=255)

  ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gg                                                                                       <  2:1153112/93‑1 (MQ=255)

  ttGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gg                                                                                       <  1:3623719/93‑1 (MQ=255)

    gATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gggg                                                                                     >  1:4831406/1‑93 (MQ=255)

      tAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gggggc                                                                                   <  1:4921062/93‑3 (MQ=255)

       aaCCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gggggct                                                                                  <  1:5024325/93‑4 (MQ=255)

        aCCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑gggggctg                                                                                 >  2:5817912/1‑89 (MQ=255)

          ccGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCa                                                                               >  1:1930402/1‑93 (MQ=255)

          ccGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATAGTCGAAGAACTTGAACAGGGGG‑CTGCa                                                                               >  2:1061978/1‑93 (MQ=255)

            gCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAATAACTTGAACAGGGGG‑CTGCaaa                                                                             >  1:6112529/1‑93 (MQ=255)

             cTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGCAGAACTTGAACAGGGGG‑CTGCaaaa                                                                            >  2:2709306/1‑93 (MQ=255)

              tGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAAt                                                                           >  1:1033004/1‑93 (MQ=255)

                gTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTa                                                                         <  1:1569789/93‑1 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGc                                                                      >  2:4694392/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGc                                                                      >  1:5910110/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGc                                                                      >  2:1224726/1‑93 (MQ=255)

                   ggggCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGc                                                                      >  1:2389385/1‑93 (MQ=255)

                      gCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCg                                                                   <  1:2581037/93‑1 (MQ=255)

                       cACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCgt                                                                  >  1:2046976/1‑93 (MQ=255)

                        aCGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCgtg                                                                 <  1:3458166/93‑1 (MQ=255)

                         cgcggcgcggCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCgtgt                                                                >  1:3756754/1‑93 (MQ=255)

                         cgcggcgcggCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCgtgt                                                                >  1:6282444/1‑93 (MQ=255)

                           cggcgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGttt                                                              <  2:5371963/93‑1 (MQ=255)

                              cgcggcgGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTac                                                           <  2:2862058/93‑1 (MQ=255)

                               gcggcggcgGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTaca                                                          >  1:3183416/1‑93 (MQ=255)

                                cggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTacaa                                                         <  2:518334/93‑1 (MQ=255)

                                cggcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTacaa                                                         <  1:1866102/93‑1 (MQ=255)

                                  gcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTacaaca                                                       <  1:1956144/93‑1 (MQ=255)

                                  gcggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTacaaca                                                       <  2:6075483/93‑1 (MQ=255)

                                    ggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGc                                                     <  1:3997220/93‑1 (MQ=255)

                                    ggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGc                                                     >  2:3778056/1‑93 (MQ=255)

                                    ggcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGc                                                     >  1:1299177/1‑93 (MQ=255)

                                     gcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCa                                                    <  1:3250406/93‑1 (MQ=255)

                                     gcggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCa                                                    >  1:4898153/1‑93 (MQ=255)

                                      cggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCaa                                                   >  2:862054/1‑93 (MQ=255)

                                       ggTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCaaa                                                  >  1:4715441/1‑93 (MQ=255)

                                        gtgtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAAc                                                 <  2:4538213/93‑1 (MQ=255)

                                         tgtTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACt                                                >  2:1923489/1‑93 (MQ=255)

                                             tGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGaaa                                            <  1:4554796/93‑1 (MQ=255)

                                              ggCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAAt                                           <  1:655176/93‑1 (MQ=255)

                                               gCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAAtt                                          >  2:242083/1‑93 (MQ=255)

                                               gCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAAtt                                          >  1:4887927/1‑93 (MQ=255)

                                                 ccGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTaa                                        <  2:1751611/93‑1 (MQ=255)

                                                 ccGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTaa                                        <  2:2735519/93‑1 (MQ=255)

                                                  cGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAAt                                       <  1:5316567/93‑1 (MQ=255)

                                                  cGCAAAAAGGCGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑GGGGACATCAAAATAAAGCACGCGTTTTACGACAGCAAAATGAAATTAAt                                       >  2:4176550/1‑93 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      <  1:2511786/93‑1 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      <  2:84225/93‑1 (MQ=255)

                                                   gCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATg                                      >  1:4088767/1‑93 (MQ=255)

                                                     aaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTc                                    >  1:2606351/1‑93 (MQ=255)

                                                     aaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTc                                    <  1:3894098/93‑1 (MQ=255)

                                                     aaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTc                                    >  1:3008856/1‑93 (MQ=255)

                                                     aaaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTc                                    <  2:4049126/93‑1 (MQ=255)

                                                      aaaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCt                                   >  1:1333478/1‑93 (MQ=255)

                                                       aaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTg                                  <  2:5211715/93‑1 (MQ=255)

                                                       aaaGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTg                                  <  2:5118382/93‑1 (MQ=255)

                                                          ggggCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCa                               >  1:2111698/1‑93 (MQ=255)

                                                            ggCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGa                             >  2:1729377/1‑93 (MQ=255)

                                                             gCAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGAt                            >  2:335052/1‑93 (MQ=255)

                                                              cAACGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATg                           <  2:625143/93‑1 (MQ=255)

                                                                aCGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggc                         >  2:4312084/1‑93 (MQ=255)

                                                                 cGCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcg                        >  1:5643375/1‑93 (MQ=255)

                                                                  gCCGGAGATGGTCGAAGGACTTGAACAGGGGG‑CCGCAAAATTATGCCCGTGTTTTTCAACAGCAAACTGAACATAATGTCTTCCAGATggcgg                       >  2:5562863/1‑93 (MQ=255)

                                                                  gCCGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcgg                       <  1:6308735/93‑1 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAACCTGAAATTAATGACTGCCAGATggcggg                      >  1:5215006/1‑93 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                      >  2:2122945/1‑93 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                      <  1:5531683/93‑1 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                      <  1:6141749/93‑1 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                      <  2:3595917/93‑1 (MQ=255)

                                                                   ccGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcggg                      <  2:1589626/93‑1 (MQ=255)

                                                                    cGGAGATGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATggcgggg                     >  1:4244726/1‑92 (MQ=255)

                                                                      gagaTGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGgcgggcgg                   >  2:372523/1‑93 (MQ=255)

                                                                      gagaTGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGgcgggcgg                   >  1:306169/1‑93 (MQ=255)

                                                                      gagaTGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGgcgggcgg                   >  1:753424/1‑93 (MQ=255)

                                                                          tGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGggcggcggc               >  2:5252126/1‑93 (MQ=255)

                                                                          tGGTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGggcggcggc               >  1:1484374/1‑93 (MQ=255)

                                                                           ggTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGgcggcggcg              <  2:2781079/93‑1 (MQ=255)

                                                                            gTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGgcgc             <  2:5344565/93‑1 (MQ=255)

                                                                            gTCGAAGAACTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGgcgc             <  1:4350768/93‑1 (MQ=255)

                                                                               gaagaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgc          <  1:1008832/93‑1 (MQ=255)

                                                                               gaagaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACGGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgc          <  2:5087483/93‑1 (MQ=255)

                                                                                aagaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcg         <  2:4390026/93‑1 (MQ=255)

                                                                                 agaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgg        <  2:4854807/93‑1 (MQ=255)

                                                                                 agaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgg        >  2:3158413/1‑93 (MQ=255)

                                                                                  gaaCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcggg       >  2:560800/1‑93 (MQ=255)

                                                                                    aCTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgggcg     >  2:1676214/1‑93 (MQ=255)

                                                                                     cTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgggcgg    >  2:3558024/1‑93 (MQ=255)

                                                                                     cTTGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTgcgggcgg    >  2:5637740/1‑93 (MQ=255)

                                                                                      ttGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGt   >  1:1092312/1‑93 (MQ=255)

                                                                                      ttGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGt   >  2:1980782/1‑93 (MQ=255)

                                                                                       tGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  <  1:2688549/93‑1 (MQ=255)

                                                                                       tGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  >  2:5422251/1‑93 (MQ=255)

                                                                                       tGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  <  1:3359207/93‑1 (MQ=255)

                                                                                       tGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  <  2:2207373/93‑1 (MQ=255)

                                                                                       tGAACAGGGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTa  <  1:3657702/93‑1 (MQ=255)

CGTTGATAACCCGCTGGTAGGGGCACGCGGCGCGGCGGCGGTGTTTGGCCCGCAAAAAGGGGCAACGCCGGAGATGGTCGAAGAACTTGAACA‑GGGG‑CTGCAAAATTACGCCCGTGTTTTACAACAGCAAACTGAAATTAATGTCTGCCAGATGGCGGGCGGCGGCGCTGCGGGCGGTA  >  NC_000913/541690‑541868

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: `‑`