BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,382	A→G	H50R (CAC→CGC)	b0005 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,382	0	A	G	99.0%	692.8 / ‑3.8	198	H50R (CAC→CGC)	b0005	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (131/65); total (132/66)
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.71e-01

Only 100 of 216 total aligned reads displayed.

CAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACC  >  NC_000913/5292‑5471

cagcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCt                                                                                         <  1:256715/93‑1 (MQ=255)

cagcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCt                                                                                         <  2:2365750/93‑1 (MQ=255)

 agcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTg                                                                                        >  1:1419382/1‑93 (MQ=255)

 agcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTg                                                                                        >  2:1103530/1‑93 (MQ=255)

  gcaCAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTgg                                                                                       >  2:3120964/1‑93 (MQ=255)

   cacaGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGgc                                                                                      <  1:3295404/93‑1 (MQ=255)

    acaGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGgcg                                                                                     >  2:5084708/1‑93 (MQ=255)

     caGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGgcgc                                                                                    <  2:4919356/93‑1 (MQ=255)

      aGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGcgcg                                                                                   >  1:1830371/1‑93 (MQ=255)

       ggCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGa                                                                                  >  1:4668102/1‑93 (MQ=255)

        gCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGAc                                                                                 >  1:1622868/1‑93 (MQ=255)

          tGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCa                                                                               <  1:1145346/93‑1 (MQ=255)

           gctgAAATTCCGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCAc                                                                              <  2:5005542/90‑1 (MQ=255)

           gCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCAc                                                                              >  1:4446491/1‑93 (MQ=255)

           gCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGGCGCTGGCGCGACCAc                                                                              >  1:2734603/1‑93 (MQ=255)

            cGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACcgcg                                                                             <  2:1147819/93‑1 (MQ=255)

             ggAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACgg                                                                            >  1:3800925/1‑93 (MQ=255)

              gAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGc                                                                           >  2:490351/1‑93 (MQ=255)

                aaTTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCtg                                                                         >  1:5243529/1‑93 (MQ=255)

                 aTTACGTTAGCCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCtgg                                                                        <  1:5763702/93‑1 (MQ=255)

                  ttACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCtggg                                                                       >  2:2180406/1‑92 (MQ=255)

                      gTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGaa                                                                   >  1:2284370/1‑93 (MQ=255)

                      gTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGaa                                                                   >  1:526167/1‑93 (MQ=255)

                          gTCTCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAaacaa                                                               <  2:6087285/93‑1 (MQ=255)

                          gTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAaacaa                                                               >  2:5073136/1‑93 (MQ=255)

                           tCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAaacaac                                                              >  2:1335926/1‑93 (MQ=255)

                           tCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAaacaac                                                              <  2:5938139/93‑1 (MQ=255)

                           gcccGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAaacaac                                                              <  1:2031295/92‑1 (MQ=255)

                            cccGTCAGTAAAATTACAGATGGGCGATCGTGATAATCGTGGCTATTACTGGGAGGGAGGTCGCTGGCGCGCCCACGGCTGGTGGAAacaaca                                                             <  2:2675535/93‑1 (MQ=255)

                            cccGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAacaaca                                                             >  2:556818/1‑93 (MQ=255)

                               gTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTa                                                          >  2:681330/1‑93 (MQ=255)

                                tCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTAt                                                         <  2:1869769/93‑1 (MQ=255)

                                tCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTAt                                                         >  1:1919518/1‑93 (MQ=255)

                                 cAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATg                                                        <  2:2293985/93‑1 (MQ=255)

                                  aGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGGGGCAACAACATTATGa                                                       >  1:5423265/1‑93 (MQ=255)

                                   gTAAAATTACAGATAGGCGATCGTGATAATCGAGGCTATTACTGGGATCGAGGTCGCGCGCGAGACCCCGGCTGGCGGCAGAAACATTATGaa                                                      >  2:2547666/1‑93 (MQ=255)

                                    tAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAAt                                                     <  1:5112534/93‑1 (MQ=255)

                                    gaaaaTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAAt                                                     >  2:5466030/2‑93 (MQ=255)

                                    aaaaaATACAGAGAGGGGATGGTGATAGTCGTGGTTTTGAGGGGGATGGGGGCCGCGGGCGCGCCCACGGCTGGTGGAAACAACATTATGAAt                                                     <  2:4511618/92‑1 (MQ=255)

                                     aaaaTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATg                                                    <  1:1994762/93‑1 (MQ=255)

                                     aaaaTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGCAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATg                                                    <  2:2327604/93‑1 (MQ=255)

                                      aaaTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATgg                                                   <  2:3275624/93‑1 (MQ=255)

                                        aTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCg                                                 <  1:2154014/93‑1 (MQ=255)

                                        aTTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCg                                                 >  1:1049020/1‑93 (MQ=255)

                                         ttACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGa                                                >  2:3247890/1‑93 (MQ=255)

                                          tACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAg                                               >  2:4762966/1‑93 (MQ=255)

                                           aCAGATAGGCGATCGTGATAATCGTGGTTATTACTGGGATGGAGGTCGCTGGCGCGCCCACGGCTGGTGGAAACAACATTATGAATGGCGAgg                                              <  2:4966779/93‑1 (MQ=255)

                                           aCAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAgg                                              >  2:2793186/1‑93 (MQ=255)

                                            cAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGc                                             <  1:4077229/93‑1 (MQ=255)

                                              gATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCaa                                           >  1:4466921/1‑93 (MQ=255)

                                              gATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCaa                                           >  2:5293768/1‑93 (MQ=255)

                                               aTAGGCGATGGTGATAATCGTGGCTATTACTGGGAGGGAGGTCGTGGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAAt                                          <  2:3952910/93‑1 (MQ=255)

                                               aTAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAAt                                          >  1:2338737/1‑93 (MQ=255)

                                                tAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATc                                         >  2:922434/1‑93 (MQ=255)

                                                tAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATc                                         >  1:3569404/1‑93 (MQ=255)

                                                 aGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCg                                        >  1:6128823/1‑93 (MQ=255)

                                                  ggCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGc                                       <  2:3177497/93‑1 (MQ=255)

                                                    ggATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCGGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTg                                     <  2:2565079/92‑1 (MQ=255)

                                                     gATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTgg                                    >  1:4387661/1‑93 (MQ=255)

                                                      aTCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGc                                   >  1:4975021/1‑93 (MQ=255)

                                                      aTCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGc                                   >  2:5342593/1‑93 (MQ=255)

                                                       tCGTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGcc                                  >  1:3587070/1‑92 (MQ=255)

                                                       tCGTGATAATCGTGGCTATTACTGGGATGGAGGGCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCa                                  >  1:2864404/1‑93 (MQ=255)

                                                         gTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCAcc                                >  1:5894803/1‑93 (MQ=255)

                                                         gTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCAcc                                >  1:5444160/1‑93 (MQ=255)

                                                         gTGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCAcc                                >  2:4290613/1‑93 (MQ=255)

                                                          tGATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCt                               >  1:4293953/1‑93 (MQ=255)

                                                           gATAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTa                              >  2:600845/1‑93 (MQ=255)

                                                            aTAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTac                             >  2:2907320/1‑93 (MQ=255)

                                                             tAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTaca                            <  1:6011212/93‑1 (MQ=255)

                                                             tAATCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTaca                            >  1:5367857/1‑93 (MQ=255)

                                                             tAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTaca                            <  1:5137401/93‑1 (MQ=255)

                                                               aTCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACg                          >  2:5654644/1‑93 (MQ=255)

                                                                tCGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACgg                         <  1:3138936/93‑1 (MQ=255)

                                                                 cGTGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGa                        <  2:1007476/93‑1 (MQ=255)

                                                                 cGTGGCTATTACGGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGa                        <  1:1921281/93‑1 (MQ=255)

                                                                   tGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAcc                      >  1:1685019/1‑93 (MQ=255)

                                                                   tGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAcc                      <  2:2916212/93‑1 (MQ=255)

                                                                   tGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAcc                      >  2:3930419/1‑93 (MQ=255)

                                                                   tGGCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACATACACGGAcc                      >  1:1122408/1‑93 (MQ=255)

                                                                    ggCTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAccc                     >  2:3329525/1‑92 (MQ=255)

                                                                      cTATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAccgcc                   >  2:3508465/1‑93 (MQ=255)

                                                                       tATTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGAcccccc                  >  1:5322013/1‑92 (MQ=255)

                                                                        aTTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACcgccgc                 <  2:4194693/93‑1 (MQ=255)

                                                                        aTTACTGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACcgccgc                 >  1:3271336/1‑93 (MQ=255)

                                                                        aTTACTGGGAGGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACcgccgc                 <  2:1478394/93‑1 (MQ=255)

                                                                           tgggggATGGGGGCCGCGGGCGCGCCCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAc              <  2:4635036/90‑1 (MQ=255)

                                                                           acggggATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAc              <  1:756230/90‑1 (MQ=255)

                                                                             tGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccg            >  2:631674/1‑93 (MQ=255)

                                                                             tGGGATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccc            >  1:2768061/1‑92 (MQ=255)

                                                                              gggATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccgc           >  1:3280857/1‑93 (MQ=255)

                                                                                gATGGGGGTGGCGGGCGCTACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccgccg         <  1:5793819/93‑1 (MQ=255)

                                                                                gATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccgccg         >  2:2860213/1‑93 (MQ=255)

                                                                                gATGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCAccgccg         >  2:1599305/1‑93 (MQ=255)

                                                                                 aTGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACcgccgc        <  2:5421270/93‑1 (MQ=255)

                                                                                  tGGAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCcgcg       >  2:364565/1‑93 (MQ=255)

                                                                                  tGGAGGTCGCTGGCGCGACCACGGCTGGGGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCcgcg       >  1:1952765/1‑93 (MQ=255)

                                                                                    gAGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGcc     <  1:3988101/93‑1 (MQ=255)

                                                                                     aGGTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGcca    >  1:4252773/1‑93 (MQ=255)

                                                                                      ggTCGCTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGccac   >  1:4348923/1‑93 (MQ=255)

CAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACC  >  NC_000913/5292‑5471

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