BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	908,737	T→C	Q512Q (CAA→CAG)	poxB ←	pyruvate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	908,737	0	T	C	99.1%	688.7 / ‑5.6	215	Q512Q (CAA→CAG)	poxB	pyruvate oxidase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (92/121); total (94/121)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

CTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGT  >  NC_000913/908647‑908827

cTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTg                                                                                          >  1:9291079/1‑93 (MQ=255)

 tGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGAGCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGc                                                                                         >  2:9867585/1‑93 (MQ=255)

 tGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGc                                                                                         <  2:4958311/93‑1 (MQ=255)

  gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                        >  1:8212359/1‑93 (MQ=255)

  gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                        <  2:9689761/93‑1 (MQ=255)

   ttCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAg                                                                                       >  1:6241075/1‑93 (MQ=255)

      gAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGc                                                                                    >  1:5838312/1‑93 (MQ=255)

      gAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGc                                                                                    <  2:8296052/93‑1 (MQ=255)

      gAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGc                                                                                    <  1:998213/93‑1 (MQ=255)

       aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                   >  1:3828343/1‑93 (MQ=255)

        gTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCtt                                                                                  <  1:4040261/93‑1 (MQ=255)

        gTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCtt                                                                                  <  1:3384969/93‑1 (MQ=255)

          ttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtca                                                                                <  2:7331185/93‑1 (MQ=255)

             aTCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatca                                                                             >  2:3680155/1‑93 (MQ=255)

               cTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                           <  1:8758185/93‑1 (MQ=255)

               cTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                           <  1:4331984/93‑1 (MQ=255)

                  cGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAGTTc                                                                        >  2:761912/1‑93 (MQ=255)

                   ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                       >  2:1602926/1‑93 (MQ=255)

                      ggAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATTAACTTCAGAc                                                                    <  2:6554412/93‑1 (MQ=255)

                      ggAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGAc                                                                    >  1:6445033/1‑93 (MQ=255)

                        aaTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGc                                                                  >  2:5719737/1‑93 (MQ=255)

                         aTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCt                                                                 <  1:5434571/93‑1 (MQ=255)

                          tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                                <  1:3960645/93‑1 (MQ=255)

                          tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                                >  1:2909414/1‑93 (MQ=255)

                          tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                                >  1:9536517/1‑93 (MQ=255)

                          tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                                <  1:5746190/93‑1 (MQ=255)

                           ggCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                               <  2:6048939/93‑1 (MQ=255)

                           ggCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                               <  2:7771283/93‑1 (MQ=255)

                           ggCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                               >  1:2475527/1‑93 (MQ=255)

                             cTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttttt                                                             >  2:8742832/1‑93 (MQ=255)

                             cTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttttt                                                             <  2:5576586/93‑1 (MQ=255)

                              tAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTc                                                            >  1:2078225/1‑93 (MQ=255)

                               aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                           >  1:1829804/1‑93 (MQ=255)

                               aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                           >  2:3532846/1‑93 (MQ=255)

                               aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                           >  2:6606841/1‑93 (MQ=255)

                               aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                           >  1:3089350/1‑93 (MQ=255)

                                aCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                          <  2:3234249/93‑1 (MQ=255)

                                aCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                          >  1:2359439/1‑93 (MQ=255)

                                 ctctTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                                                         >  2:596134/1‑93 (MQ=255)

                                 ctctTCTCTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACTCTTTTTCTac                                                         <  1:1298991/93‑1 (MQ=255)

                                  tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacc                                                        >  1:8381926/1‑92 (MQ=255)

                                  tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                        <  1:8319837/93‑1 (MQ=255)

                                  tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                        >  2:8859944/1‑93 (MQ=255)

                                  tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                        <  2:9534162/93‑1 (MQ=255)

                                   cttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                                                       >  2:3980611/1‑93 (MQ=255)

                                   cttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                                                       <  1:7331052/93‑1 (MQ=255)

                                       tttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGAt                                                   <  1:4126099/93‑1 (MQ=255)

                                       tttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGAt                                                   <  1:9036638/93‑1 (MQ=255)

                                        ttGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATa                                                  >  2:1889491/1‑93 (MQ=255)

                                         tGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAATTCAGACGCTTTTTCTACACGGATAc                                                 <  1:531243/93‑1 (MQ=255)

                                           gCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACTCTTTTTCTACACGGATAccc                                               <  1:4753672/93‑1 (MQ=255)

                                               ccaccaCATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTaa                                           <  1:5919032/93‑1 (MQ=255)

                                               ccaccaCATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTaa                                           <  2:9421227/93‑1 (MQ=255)

                                               ccaccaCATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTaa                                           >  2:6202000/1‑93 (MQ=255)

                                                caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                          <  1:9135974/93‑1 (MQ=255)

                                                caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                          >  2:3063787/1‑93 (MQ=255)

                                                caccacATCCACCAACACCGGACCGGCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTaag                                          >  2:9602820/1‑92 (MQ=255)

                                                 cacacaTCCCCCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATg                                         <  2:8238720/91‑1 (MQ=255)

                                                 accacATCCACCAACGCCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATg                                         <  1:7082961/93‑1 (MQ=255)

                                                      aaccaccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCa                                    <  2:8953458/91‑1 (MQ=255)

                                                       tCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCAc                                   >  2:2473363/1‑93 (MQ=255)

                                                       tCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCAc                                   >  1:1563147/1‑93 (MQ=255)

                                                        ccaccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCcgcgcg                                  >  1:4923827/1‑93 (MQ=255)

                                                        ccaccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACg                                  <  2:7888680/93‑1 (MQ=255)

                                                         caccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGc                                 <  2:9655271/93‑1 (MQ=255)

                                                         caccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGc                                 >  2:4306704/1‑93 (MQ=255)

                                                          cccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCt                                <  2:9329747/92‑1 (MQ=255)

                                                            caACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTc                              >  2:9887192/1‑93 (MQ=255)

                                                            caACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTc                              >  2:4746560/1‑93 (MQ=255)

                                                            caACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTc                              >  2:4746557/1‑93 (MQ=255)

                                                             aaCACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCg                             <  1:7375596/93‑1 (MQ=255)

                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                         <  1:2506694/93‑1 (MQ=255)

                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                         <  1:1695156/93‑1 (MQ=255)

                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                         >  2:9232839/1‑93 (MQ=255)

                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                         <  2:1735828/93‑1 (MQ=255)

                                                                   ggACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCGACACGGATACCCGTAATGCCGCACGCTTCGGCAATg                       >  2:8460675/1‑93 (MQ=255)

                                                                     gccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCg                     <  1:3212903/92‑1 (MQ=255)

                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATGCCCGTAATGCCGCACGCTTCGGCAATGCgg                    >  2:5050135/1‑93 (MQ=255)

                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCgg                    <  1:2605079/93‑1 (MQ=255)

                                                                       cGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCggg                   >  1:4691418/1‑93 (MQ=255)

                                                                         tCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCa                 >  2:528816/1‑93 (MQ=255)

                                                                          cGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCaa                <  1:6123550/93‑1 (MQ=255)

                                                                           gATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCaaa               <  2:4525922/93‑1 (MQ=255)

                                                                           gATGGACAAGGCGCGCTGCAGCGCTTCATCAACTCCAGACGCTTTTTCTACACGGATCCCCGTAATGCCGCCCGCTTCGGCAATGCGGGCaaa               <  2:1436680/93‑1 (MQ=255)

                                                                             tGGAGAAGGCGCGCTGCAGGGCTTCATCCACTTCAGACGCTTTTTCTACCCGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGt             >  1:7008888/1‑93 (MQ=255)

                                                                              ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGtt            <  2:3849121/93‑1 (MQ=255)

                                                                              ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCATAATGCCGCACGCTTCGGCAATGCGGGCAAAGtt            >  1:4613357/1‑93 (MQ=255)

                                                                               gagaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGttt           <  1:6949944/93‑1 (MQ=255)

                                                                                agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtg          <  1:7698694/93‑1 (MQ=255)

                                                                                agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtg          >  1:5365487/1‑93 (MQ=255)

                                                                                 gaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtgt         <  1:6729384/93‑1 (MQ=255)

                                                                                 gaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtgt         <  1:1054803/93‑1 (MQ=255)

                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt       >  2:5914926/1‑93 (MQ=255)

                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt       <  1:9393916/93‑1 (MQ=255)

                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt       <  1:3774744/93‑1 (MQ=255)

                                                                                    ggCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTc      <  2:6467033/93‑1 (MQ=255)

                                                                                      cgcgCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCgt    <  1:8008550/93‑1 (MQ=255)

                                                                                      cgcgCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCgt    <  2:7159922/93‑1 (MQ=255)

                                                                                       gcgcTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCgtg   <  1:9061254/93‑1 (MQ=255)

                                                                                        cgTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCgtgt  >  2:5112085/1‑93 (MQ=255)

CTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGT  >  NC_000913/908647‑908827

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