BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,364	A→G	H339R (CAC→CGC)	gsk →	inosine‑guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,364	0	A	G	83.0%	595.7 / 117.2	271	H339R (CAC→CGC)	gsk	inosine‑guanosine kinase
Reads supporting (aligned to +/- strand): ref base A (19/27); new base G (104/121); total (123/148)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 281 total aligned reads displayed.

CTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCG  >  NC_000913/500273‑500456

 tgcCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGcaca                                                                                            <  2:4994227/93‑1 (MQ=255)

   ccGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCAtt                                                                                          <  1:1094698/93‑1 (MQ=255)

   ccGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACAtt                                                                                          <  2:6977256/93‑1 (MQ=255)

     gggCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTgc                                                                                        <  1:2178756/93‑1 (MQ=255)

        cgcTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGcc                                                                                     <  1:834790/93‑1 (MQ=255)

         gcTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACACGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCg                                                                                    >  2:3156698/1‑93 (MQ=255)

         gcTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCg                                                                                    >  1:6472742/1‑93 (MQ=255)

         gcTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCg                                                                                    >  1:663235/1‑93 (MQ=255)

          cTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGt                                                                                   >  1:3639164/1‑93 (MQ=255)

             taGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACa                                                                                >  1:1729445/1‑93 (MQ=255)

             taGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACa                                                                                >  2:2726687/1‑93 (MQ=255)

              aGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACAt                                                                               <  2:5652951/93‑1 (MQ=255)

               gCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATg                                                                              >  2:4084100/1‑93 (MQ=255)

               gCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATg                                                                              <  1:471418/93‑1 (MQ=255)

               gCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATg                                                                              <  2:6046555/93‑1 (MQ=255)

               gCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATg                                                                              >  1:1579070/1‑93 (MQ=255)

                 ggAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATggg                                                                            >  1:4168665/1‑93 (MQ=255)

                 ggAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATggg                                                                            >  1:4792371/1‑93 (MQ=255)

                   aaTTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATgggcg                                                                          <  2:3637088/93‑1 (MQ=255)

                    aTTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATgggcgg                                                                         <  1:6893457/93‑1 (MQ=255)

                    aTTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgggcgg                                                                         <  2:5491375/93‑1 (MQ=255)

                      tCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATgggcgggc                                                                       <  2:850853/93‑1 (MQ=255)

                       cAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGcc                                                                      >  2:7513028/1‑93 (MQ=255)

                        aaCCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCg                                                                     >  2:5561640/1‑93 (MQ=255)

                        aaCCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCg                                                                     >  1:7290443/1‑93 (MQ=255)

                         aCCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCgg                                                                    <  2:65767/93‑1 (MQ=255)

                         aCCAGTATGAGTTTAGCCGCGCCATGCGCCAC‑AGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGa                                                                   >  1:5906905/1‑93 (MQ=255)

                          ccAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGa                                                                   >  2:9718/1‑93 (MQ=255)

                            aGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGaaa                                                                 >  1:5030489/1‑93 (MQ=255)

                               aTGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGaaaaaa                                                              <  2:5865825/93‑1 (MQ=255)

                               aTGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGaaaaaa                                                              >  1:5176299/1‑93 (MQ=255)

                                tGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAAt                                                             >  2:659415/1‑93 (MQ=255)

                                tGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAAt                                                             <  2:5409643/93‑1 (MQ=255)

                                tGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAAt                                                             <  2:4469098/93‑1 (MQ=255)

                                tGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAAt                                                             >  1:3025310/1‑93 (MQ=255)

                                  aGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCa                                                           <  1:1086110/93‑1 (MQ=255)

                                   gTTCAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCAt                                                          >  1:6978587/1‑93 (MQ=255)

                                    tttAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATg                                                         <  1:3276515/93‑1 (MQ=255)

                                    tttAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATg                                                         >  2:2801532/1‑93 (MQ=255)

                                        gCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAaca                                                     >  1:699584/1‑93 (MQ=255)

                                            cgcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTaa                                                 <  2:1155769/93‑1 (MQ=255)

                                            cgcCATGCGCCACAAGGATTGCCAGAACCCGCTGCGTGTATATTAGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTaa                                                 <  2:3525501/93‑1 (MQ=255)

                                             gcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAAt                                                <  1:1907379/93‑1 (MQ=255)

                                             gcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAAt                                                >  2:1324918/1‑93 (MQ=255)

                                              ccATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATg                                               <  1:3567478/93‑1 (MQ=255)

                                              ccATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATg                                               <  2:7218894/93‑1 (MQ=255)

                                               cATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATgg                                              <  1:6212337/93‑1 (MQ=255)

                                               cATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATgg                                              <  1:2221529/93‑1 (MQ=255)

                                                aTGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGa                                             <  2:3341716/93‑1 (MQ=255)

                                                aTGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGa                                             <  2:1709924/93‑1 (MQ=255)

                                                 tGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAg                                            <  1:2039372/93‑1 (MQ=255)

                                                   cgcCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCg                                          <  1:2917764/93‑1 (MQ=255)

                                                    gcCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCgg                                         >  2:2465925/1‑93 (MQ=255)

                                                     ccACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCggg                                        <  2:7022526/93‑1 (MQ=255)

                                                      cacaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCgggg                                       >  2:1531504/1‑93 (MQ=255)

                                                      cacaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCgggg                                       >  2:6848385/1‑93 (MQ=255)

                                                      cacaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCgggg                                       >  2:113918/1‑93 (MQ=255)

                                                      cacaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCgggg                                       >  2:7434595/1‑93 (MQ=255)

                                                       acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCggggg                                      >  1:6030154/1‑93 (MQ=255)

                                                       acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCggggg                                      >  1:6666697/1‑93 (MQ=255)

                                                       acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCggggg                                      >  2:6199382/1‑93 (MQ=255)

                                                        caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGa                                     >  2:197569/1‑93 (MQ=255)

                                                        caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGa                                     <  1:5896168/93‑1 (MQ=255)

                                                        caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGa                                     <  1:4080464/93‑1 (MQ=255)

                                                        caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGa                                     <  1:5640730/93‑1 (MQ=255)

                                                         aaGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGAt                                    <  1:2668808/93‑1 (MQ=255)

                                                         aaGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGAt                                    <  1:1312340/93‑1 (MQ=255)

                                                          aGGATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATg                                   >  2:2886474/1‑93 (MQ=255)

                                                           ggATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATgg                                  <  2:6761893/93‑1 (MQ=255)

                                                           ggATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATgg                                  >  2:1916638/1‑93 (MQ=255)

                                                            gATTGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgc                                 <  1:5559393/93‑1 (MQ=255)

                                                              ttGCCAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgcgc                               >  2:6780488/1‑93 (MQ=255)

                                                                gCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCAt                             <  1:984672/93‑1 (MQ=255)

                                                                 ccAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCAtt                            <  1:6306498/93‑1 (MQ=255)

                                                                 ccAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCAtt                            <  1:6780837/93‑1 (MQ=255)

                                                                  cAGAATCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTg                           >  1:4542126/1‑93 (MQ=255)

                                                                   aGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTgg                          <  2:212102/93‑1 (MQ=255)

                                                                      aTCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAg                       >  2:6268992/1‑93 (MQ=255)

                                                                      aTCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAg                       <  1:3202370/93‑1 (MQ=255)

                                                                      aTCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAg                       <  2:2256534/93‑1 (MQ=255)

                                                                      aTCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAg                       >  2:2635039/1‑93 (MQ=255)

                                                                       tCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGc                      <  2:4765233/93‑1 (MQ=255)

                                                                       tCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGc                      <  2:4392064/93‑1 (MQ=255)

                                                                       tCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGc                      >  2:3772992/1‑93 (MQ=255)

                                                                       tCCGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGc                      <  2:252091/93‑1 (MQ=255)

                                                                         cGCTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGt                    <  2:3429877/93‑1 (MQ=255)

                                                                           cTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtg                  >  2:3615604/1‑93 (MQ=255)

                                                                           cTGCGTGTATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtg                  <  2:6707911/93‑1 (MQ=255)

                                                                               gtgtATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctgc              <  1:7442744/93‑1 (MQ=255)

                                                                               gtgtATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctgc              >  1:4336277/1‑93 (MQ=255)

                                                                                tgtATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCa             >  2:4168022/1‑93 (MQ=255)

                                                                                 gtATATTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCAt            >  2:7154151/1‑93 (MQ=255)

                                                                                    tatTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGTATGAc         >  1:2462903/1‑93 (MQ=255)

                                                                                    tatTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGAc         <  1:4913568/93‑1 (MQ=255)

                                                                                    tatTCGCGCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGAc         <  2:5122823/93‑1 (MQ=255)

                                                                                      ttCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACAt       >  2:3118446/1‑93 (MQ=255)

                                                                                      ttCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACAt       <  1:1849634/93‑1 (MQ=255)

                                                                                        cgcgCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTa     >  1:3449490/1‑93 (MQ=255)

                                                                                         gCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTAc    <  1:2548820/93‑1 (MQ=255)

                                                                                           acaTTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCg  >  2:2364066/1‑93 (MQ=255)

CTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCG  >  NC_000913/500273‑500456

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