BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,056,882	T→C	Y99C (TAC→TGC)	torR ←	TorR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,056,882	0	T	C	99.5%	579.0 / ‑5.0	183	Y99C (TAC→TGC)	torR	TorR
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (87/95); total (87/96)
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.57e-01

Only 100 of 186 total aligned reads displayed.

TGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTGAGCGTTCT  >  NC_000913/1056792‑1056971

tgCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCaa                                                                                          <  1:6120105/93‑1 (MQ=255)

   ggTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCg                                                                                       <  2:5482259/93‑1 (MQ=255)

     ttGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTc                                                                                     <  1:3539730/93‑1 (MQ=255)

      tGAGCTTTTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCt                                                                                    <  1:1797450/93‑1 (MQ=255)

      tGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCt                                                                                    <  1:723464/93‑1 (MQ=255)

      tGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCt                                                                                    <  1:878676/93‑1 (MQ=255)

        aGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgc                                                                                  <  2:1334858/93‑1 (MQ=255)

        aGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgc                                                                                  >  1:2454092/1‑93 (MQ=255)

         gCTTGTCGCGCGAGGTCGACTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTgcg                                                                                 >  1:2456879/1‑93 (MQ=255)

            tGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGccc                                                                              >  2:3238903/1‑93 (MQ=255)

            tGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGccc                                                                              <  2:3984900/93‑1 (MQ=255)

             gTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCa                                                                             <  1:1757909/93‑1 (MQ=255)

             gTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCa                                                                             >  2:2683933/1‑93 (MQ=255)

               cgcgcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCAtt                                                                           >  1:980297/1‑93 (MQ=255)

                gcgcgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCAttt                                                                          >  1:2902094/1‑93 (MQ=255)

                   cgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCa                                                                       <  1:5126254/93‑1 (MQ=255)

                   cgAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCa                                                                       >  2:3928768/1‑93 (MQ=255)

                    gAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAg                                                                      >  1:337436/1‑93 (MQ=255)

                       gTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGccc                                                                   >  1:221061/1‑93 (MQ=255)

                       gTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGccc                                                                   >  1:5335201/1‑93 (MQ=255)

                       gTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGccc                                                                   >  2:3561392/1‑93 (MQ=255)

                        tCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTTCAGCCca                                                                  >  1:3384274/1‑93 (MQ=255)

                         cGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaa                                                                 <  1:3975532/93‑1 (MQ=255)

                         cGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaa                                                                 <  1:3410038/93‑1 (MQ=255)

                          gATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaac                                                                <  2:3677146/93‑1 (MQ=255)

                          gATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaac                                                                >  1:1722839/1‑93 (MQ=255)

                            ttCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCcaacaa                                                              >  1:2028143/1‑93 (MQ=255)

                             tCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAAt                                                             <  1:2095809/93‑1 (MQ=255)

                             tCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAAt                                                             <  2:3350608/93‑1 (MQ=255)

                              cGCCAGAGCAGATTTTTCCCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATa                                                            <  1:5956484/93‑1 (MQ=255)

                               gCCAGAGCATATTTTTCACCCGTACTACCAGTGCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCCCCCATTTCCAGCCCAACAATAc                                                           <  2:552616/93‑1 (MQ=255)

                               gCCAGAGCAGATTTTTCACCCGTACTCCCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATAc                                                           <  1:5926160/93‑1 (MQ=255)

                               gCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATAc                                                           >  2:3352856/1‑93 (MQ=255)

                                ccAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACg                                                          <  1:3595444/93‑1 (MQ=255)

                                 cAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCCACAATACgg                                                         >  1:4685253/1‑93 (MQ=255)

                                 cAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACgg                                                         <  1:4050768/93‑1 (MQ=255)

                                 cAGAGCAGATTTTTCACCCGAACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACgg                                                         >  2:289116/1‑93 (MQ=255)

                                  agagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGt                                                        >  1:5591939/1‑93 (MQ=255)

                                  agagCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGt                                                        <  1:1540493/93‑1 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                     <  1:5692801/93‑1 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                     <  1:43827/93‑1 (MQ=255)

                                     gCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTGTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                     <  1:2976915/93‑1 (MQ=255)

                                     acAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCaa                                                     <  2:4237808/92‑1 (MQ=255)

                                      cAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAAt                                                    <  2:1425129/93‑1 (MQ=255)

                                      cAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAAt                                                    <  2:1172942/93‑1 (MQ=255)

                                      cAGAGTTTTCTCCCGGACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAAt                                                    <  2:2005844/93‑1 (MQ=255)

                                       aGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATc                                                   >  1:4803294/1‑93 (MQ=255)

                                        gATTTTTCCCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCGGCGCCCATTTCCAGCCCAACAATACGGTCAATcc                                                  <  2:594771/93‑1 (MQ=255)

                                        gATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATcc                                                  <  2:2154657/93‑1 (MQ=255)

                                        gATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATcc                                                  <  1:1912470/93‑1 (MQ=255)

                                        gATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATcc                                                  <  2:4670501/93‑1 (MQ=255)

                                           ttttCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGAt                                               <  1:5916419/93‑1 (MQ=255)

                                              tCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGc                                            <  2:1117774/93‑1 (MQ=255)

                                              tCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGc                                            <  2:2787857/93‑1 (MQ=255)

                                               cACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCt                                           <  1:835160/93‑1 (MQ=255)

                                                aCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTg                                          >  2:2323252/1‑93 (MQ=255)

                                                aCCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTg                                          >  2:950951/1‑93 (MQ=255)

                                                 cccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGc                                         >  1:5766966/1‑93 (MQ=255)

                                                 cccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGc                                         >  1:4984977/1‑93 (MQ=255)

                                                  ccGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCg                                        >  1:5433313/1‑93 (MQ=255)

                                                   cGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGt                                       >  2:6197239/1‑93 (MQ=255)

                                                    gTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTc                                      <  1:4768855/93‑1 (MQ=255)

                                                    gTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTc                                      >  2:5691931/1‑93 (MQ=255)

                                                    gTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTc                                      <  1:3730455/93‑1 (MQ=255)

                                                     tactacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTcc                                     >  1:4183363/1‑93 (MQ=255)

                                                       ctacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCgg                                   <  1:5479522/93‑1 (MQ=255)

                                                       ctacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCgg                                   <  1:5049539/93‑1 (MQ=255)

                                                         acCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTa                                 <  2:4471530/93‑1 (MQ=255)

                                                         acCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTa                                 >  1:3074035/1‑93 (MQ=255)

                                                          ccAGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTaa                                >  1:990845/1‑93 (MQ=255)

                                                            aGTTCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAAcc                              >  2:3013979/1‑93 (MQ=255)

                                                              gtCGCGCAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCATCCCAACAATACGGTCAATCCGATCGCTGCGTCCCGTAACC‑Ag                           <  2:1872016/92‑1 (MQ=255)

                                                                   gcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGaataa                      >  1:3107045/1‑93 (MQ=255)

                                                                   gcAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGaataa                      <  2:1414345/93‑1 (MQ=255)

                                                                    cAGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGaataat                     >  1:1734515/1‑93 (MQ=255)

                                                                     aGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATc                    <  2:5778012/93‑1 (MQ=255)

                                                                     aGTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATc                    <  1:102191/93‑1 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATcc                   >  2:3821716/1‑93 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATcc                   >  2:4116873/1‑93 (MQ=255)

                                                                      gTTCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATcc                   >  1:2153347/1‑93 (MQ=255)

                                                                       ttCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATccc                  >  2:2475787/1‑93 (MQ=255)

                                                                        tCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATcccc                 >  2:5396927/1‑93 (MQ=255)

                                                                        tCCAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATcccc                 <  1:5879618/93‑1 (MQ=255)

                                                                         ccAGGGGTTTGGTGGCGCAATCGTCTGGGCCCATTTCCAGCCCAACAAAACGGGCAATCCGGCGCCTGGGTCCGGGAACCAAAAAAAATcccc                 >  2:1634333/1‑93 (MQ=255)

                                                                         ccAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCa                <  1:461842/93‑1 (MQ=255)

                                                                          cAGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCAc               >  1:673290/1‑93 (MQ=255)

                                                                           aGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCAcc              <  2:1802640/93‑1 (MQ=255)

                                                                           aGCGGTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCAcc              >  2:507691/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          >  2:5632243/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          >  1:1000042/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          >  1:4296479/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          >  1:4348694/1‑93 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          <  1:4572443/93‑1 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          <  1:5338469/93‑1 (MQ=255)

                                                                               gTTTGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGGAACC‑AGAATAATCCGCACCGTTg          >  1:1493385/1‑93 (MQ=255)

                                                                               gTTTAGTGACGCAATCGTCTCCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTg          <  2:1311779/93‑1 (MQ=255)

                                                                                tttGGTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGATAACC‑AGAATAATCCCCACCGTTGa         >  1:1103026/1‑93 (MQ=255)

                                                                                   ggTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTGAGCg      >  1:5157245/1‑93 (MQ=255)

                                                                                    gTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTGAGCGt     <  2:4373482/93‑1 (MQ=255)

                                                                                    gTGACGCAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTGAGCGt     >  1:2714728/1‑93 (MQ=255)

TGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACC‑AGAATAATCCCCACCGTTGAGCGTTCT  >  NC_000913/1056792‑1056971

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