BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	804,639	C→T	G646G (GGC→GGT)	ybhJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	804,639	0	C	T	99.6%	782.1 / ‑5.8	223	G646G (GGC→GGT)	ybhJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (121/100); total (123/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 230 total aligned reads displayed.

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

gtTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGa                                                                                         <  1:1149241/93‑1 (MQ=255)

gtTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGa                                                                                         >  1:6209475/1‑93 (MQ=255)

gtTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGa                                                                                         <  2:4597823/93‑1 (MQ=255)

 tttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGAt                                                                                        <  2:1453104/93‑1 (MQ=255)

 tttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGAt                                                                                        <  1:3023117/93‑1 (MQ=255)

  ttGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATg                                                                                       <  2:389637/93‑1 (MQ=255)

   tGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATgg                                                                                      <  2:2003574/93‑1 (MQ=255)

   tGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATgg                                                                                      <  2:2659768/93‑1 (MQ=255)

    gcgcgcATTAAGCAGATTGCTGGTCAGGAGGATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGt                                                                                     <  2:2706732/93‑1 (MQ=255)

    gcgcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGt                                                                                     <  2:4660619/93‑1 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  2:2561167/93‑1 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  2:6052026/93‑1 (MQ=255)

      gcgcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTc                                                                                   <  1:1041751/93‑1 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  1:2371658/1‑93 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  2:5524997/1‑93 (MQ=255)

        gcATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTg                                                                                 >  2:6221418/1‑93 (MQ=255)

          aTTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCCGCATGGTCTATGCGGTGAAACCAGGTGATGGCTCTgcg                                                                               >  2:3725834/1‑93 (MQ=255)

          aTTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgcg                                                                               >  1:3712241/1‑93 (MQ=255)

           ttAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTgcgc                                                                              >  1:4506348/1‑93 (MQ=255)

              aGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTg                                                                           >  2:6568665/1‑93 (MQ=255)

              aGCAGATTGCTGGTCAGGAGCATATTGATACGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTg                                                                           >  1:6405999/1‑93 (MQ=255)

               gCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGa                                                                          >  2:1252302/1‑93 (MQ=255)

               gCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGa                                                                          >  1:5221777/1‑93 (MQ=255)

               gCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGa                                                                          >  1:2048936/1‑93 (MQ=255)

               gCAGATTGCTGGTCAGGAGAATATTGATCCGCTGCAAACTGAAATTTGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                                                          >  1:4147991/1‑93 (MQ=255)

                 aGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAAc                                                                        >  1:5382960/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  1:2872637/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  1:3834514/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  2:5070296/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  2:2757323/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       >  1:3449965/1‑93 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       <  1:5104621/93‑1 (MQ=255)

                  gATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACa                                                                       <  1:4066537/93‑1 (MQ=255)

                    ttGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAgg                                                                     >  2:6119143/1‑93 (MQ=255)

                        tGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAggcggc                                                                 <  1:350355/93‑1 (MQ=255)

                             aGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                                                            >  2:3595836/1‑93 (MQ=255)

                             aGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                                                            >  1:6496499/1‑93 (MQ=255)

                              ggAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                                                           >  2:816491/1‑93 (MQ=255)

                               gAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc                                                          >  1:591924/1‑93 (MQ=255)

                                aGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGcc                                                         <  2:2935243/93‑1 (MQ=255)

                                 gCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                                                        >  2:3446906/1‑93 (MQ=255)

                                 gCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                                                        >  2:6065959/1‑93 (MQ=255)

                                    tatTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                                                     <  2:5107617/93‑1 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    >  2:5435918/1‑93 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    <  2:3248842/93‑1 (MQ=255)

                                     atTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt                                                    >  1:4030336/1‑93 (MQ=255)

                                      ttGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg                                                   >  1:5619007/1‑93 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  <  2:4731278/93‑1 (MQ=255)

                                       tGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt                                                  <  1:5823129/93‑1 (MQ=255)

                                        gATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGtgtg                                                 <  2:5264676/93‑1 (MQ=255)

                                        gATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGtgtg                                                 <  1:4428339/93‑1 (MQ=255)

                                         aTCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa                                                >  1:3049163/1‑93 (MQ=255)

                                          tCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt                                               <  1:61258/93‑1 (MQ=255)

                                           ccGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt                                              <  1:2634310/93‑1 (MQ=255)

                                           ccGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt                                              <  2:5492667/93‑1 (MQ=255)

                                            cGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg                                             <  1:3985086/93‑1 (MQ=255)

                                             gCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg                                            >  1:2484024/1‑93 (MQ=255)

                                               tGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCg                                          >  2:871819/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         >  2:6362668/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         >  2:5919882/1‑93 (MQ=255)

                                                gCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg                                         <  1:6582278/93‑1 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTcg                                      >  1:4804036/1‑92 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt                                      >  1:3189040/1‑93 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt                                      >  2:2466122/1‑93 (MQ=255)

                                                   aaCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt                                      <  1:885677/93‑1 (MQ=255)

                                                    aCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTg                                     >  1:5899066/1‑93 (MQ=255)

                                                     cTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTgg                                    >  2:2447073/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   >  1:4643496/1‑93 (MQ=255)

                                                      tGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGc                                   >  2:6094033/1‑93 (MQ=255)

                                                       gAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCg                                  >  1:4395423/1‑93 (MQ=255)

                                                       gAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCg                                  >  2:2080839/1‑93 (MQ=255)

                                                         aaTTGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGaa                                >  2:5859391/1‑93 (MQ=255)

                                                           ttGGCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAata                              >  2:1097989/1‑93 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            <  2:4379050/93‑1 (MQ=255)

                                                             ggCAGCATGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATAtt                            <  2:3190537/93‑1 (MQ=255)

                                                              gcagcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTg                           <  2:4272933/93‑1 (MQ=255)

                                                                agcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGCGATTGGCGGTCTGGCGAATATTGcc                         >  1:3881379/1‑93 (MQ=255)

                                                                 gcaTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCg                        >  2:746566/1‑93 (MQ=255)

                                                                  caTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCga                       <  1:1700661/93‑1 (MQ=255)

                                                                   aTGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgag                      >  2:3906137/1‑93 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  1:6677159/1‑93 (MQ=255)

                                                                    tGGTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagg                     >  2:3991480/1‑93 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgaggg                    >  2:1634965/1‑92 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  2:197354/93‑1 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  1:221382/93‑1 (MQ=255)

                                                                     ggTCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCgagga                    <  2:1772371/93‑1 (MQ=255)

                                                                       tCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGt                  >  2:1000364/1‑93 (MQ=255)

                                                                       tCTATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGt                  <  1:1917464/93‑1 (MQ=255)

                                                                         tATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAc                >  2:5149549/1‑93 (MQ=255)

                                                                         tATGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAc                >  1:689813/1‑93 (MQ=255)

                                                                          aTGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTAcg               <  1:46499/93‑1 (MQ=255)

                                                                          aTGCGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGGAcg               >  1:1410510/1‑93 (MQ=255)

                                                                             cGGTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGa            >  2:5781555/1‑93 (MQ=255)

                                                                               gTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACt          <  1:6024962/93‑1 (MQ=255)

                                                                               gTGAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGGGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACt          >  1:6176562/1‑93 (MQ=255)

                                                                                 gAAACCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTaa        >  2:3826123/1‑93 (MQ=255)

                                                                                  aaaCCAGGTGATGGTTCTGCGCGTGACCAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTaaa       <  1:2795387/93‑1 (MQ=255)

                                                                                    aCCAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGTCTAAACg     <  1:6607223/93‑1 (MQ=255)

                                                                                     ccAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGc    >  2:464126/1‑93 (MQ=255)

                                                                                      cAGGTGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCt   <  2:4159499/93‑1 (MQ=255)

GTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGGCGAATATTGCCGAGGAGTACGCGACTAAACGCTA  >  NC_000913/804549‑804728

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