BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	374,060	G→A	E323E (GAG→GAA)	mhpE →	4‑hydroxy‑2‑ketovalerate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	374,060	0	G	A	100.0%	563.9 / NA	156	E323E (GAG→GAA)	mhpE	4‑hydroxy‑2‑ketovalerate aldolase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (80/76); total (80/76)

Only 100 of 167 total aligned reads displayed.

TCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAC  >  NC_000913/373970‑374149

tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                         >  1:2252608/1‑93 (MQ=255)

tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                         <  1:3356851/93‑1 (MQ=255)

tCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAggaaga                                                                                         >  1:2582510/1‑93 (MQ=255)

   cTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGatat                                                                                      <  2:1433866/93‑1 (MQ=255)

   cTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGatat                                                                                      >  2:3478376/1‑93 (MQ=255)

    tgtgAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATg                                                                                     >  1:675317/1‑93 (MQ=255)

    tgtgAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATg                                                                                     >  2:2343812/1‑93 (MQ=255)

       gAATCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATGCTCGTAGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                  <  2:1349896/93‑1 (MQ=255)

       gAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                  <  1:2215918/93‑1 (MQ=255)

       gAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATc                                                                                  <  2:4496849/93‑1 (MQ=255)

        aatcggcggcggcggGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCg                                                                                 <  2:2855384/90‑1 (MQ=255)

        aaaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCg                                                                                 >  1:2104903/1‑93 (MQ=255)

         aaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                >  2:3857651/1‑93 (MQ=255)

         aaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                <  2:1039762/93‑1 (MQ=255)

         aaCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGt                                                                                >  1:701017/1‑93 (MQ=255)

          aCGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGtt                                                                               <  2:341513/93‑1 (MQ=255)

            ggcggcggcGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGa                                                                             <  2:2726831/93‑1 (MQ=255)

            ggcggcggcGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGa                                                                             <  2:1026467/93‑1 (MQ=255)

              gggcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACg                                                                           <  2:520687/92‑1 (MQ=255)

                gcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTg                                                                         >  1:89022/1‑93 (MQ=255)

                gcggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTg                                                                         >  2:3831507/1‑93 (MQ=255)

                  ggcgCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGgc                                                                       <  1:4315548/93‑1 (MQ=255)

                    cgcgTTATGGCTTAAGTGCGGTGGATATTGTCGTTGAGCTGGGCACACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGgcgc                                                                     <  1:4207533/93‑1 (MQ=255)

                     gcgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCt                                                                    <  1:3443246/93‑1 (MQ=255)

                      cgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTg                                                                   <  2:1850203/93‑1 (MQ=255)

                      cgTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTg                                                                   >  2:2694023/1‑93 (MQ=255)

                       gTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTgg                                                                  >  1:3209358/1‑93 (MQ=255)

                       gTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTgg                                                                  >  1:19565/1‑93 (MQ=255)

                       gTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTgg                                                                  >  2:605883/1‑93 (MQ=255)

                        ttATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGa                                                                 >  1:3077473/1‑93 (MQ=255)

                        ttATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGa                                                                 >  1:3077472/1‑93 (MQ=255)

                         tATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGAt                                                                >  2:135247/1‑93 (MQ=255)

                          aTGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATc                                                               >  1:2976450/1‑93 (MQ=255)

                           tGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCt                                                              >  2:1750399/1‑93 (MQ=255)

                            ggCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTg                                                             <  1:297585/93‑1 (MQ=255)

                             gCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgc                                                            <  2:3695186/93‑1 (MQ=255)

                               ttAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTTGCGCTGGATCTgcgc                                                          >  2:1235663/1‑93 (MQ=255)

                               ttAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgcgc                                                          >  2:81580/1‑93 (MQ=255)

                               ttAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTgcgc                                                          >  2:2031124/1‑93 (MQ=255)

                                tAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGca                                                         >  1:4033254/1‑93 (MQ=255)

                                 aaGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaa                                                        <  1:2955137/93‑1 (MQ=255)

                                 aaGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaa                                                        >  1:2987328/1‑93 (MQ=255)

                                  aGTGCGTTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaac                                                       <  2:1701206/93‑1 (MQ=255)

                                     gCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaacaac                                                    <  1:51297/93‑1 (MQ=255)

                                       ggTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGcaacaacaa                                                  >  2:12986/1‑93 (MQ=255)

                                        gTGGATATGCTGGTTGAGCGCGGCAAACGGCGGATGGTCGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACaaa                                                 <  2:1161055/93‑1 (MQ=255)

                                         tGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGGAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaat                                                >  2:4096645/1‑93 (MQ=255)

                                         tGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaat                                                <  1:1363428/93‑1 (MQ=255)

                                         tGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaat                                                <  2:3603425/93‑1 (MQ=255)

                                          ggATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaata                                               >  2:554905/1‑93 (MQ=255)

                                           gATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaataa                                              >  1:4637079/1‑93 (MQ=255)

                                           gATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaataa                                              >  1:107894/1‑93 (MQ=255)

                                           gATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaataa                                              <  1:3153314/93‑1 (MQ=255)

                                            atatCCTCGTTGAGCTGGGCAATCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAaataat                                             <  2:2545062/93‑1 (MQ=255)

                                             tatTCTCGTTGAGCTGGGCAAACGGCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatg                                            <  1:2670482/93‑1 (MQ=255)

                                             gatTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatg                                            <  1:1492259/92‑1 (MQ=255)

                                               ttCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAatgat                                          <  2:2911388/93‑1 (MQ=255)

                                                 ctcGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAAtgatga                                        <  1:2583089/93‑1 (MQ=255)

                                                  tcGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGAc                                       >  2:1470771/1‑93 (MQ=255)

                                                  tcGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGAc                                       >  2:2812881/1‑93 (MQ=255)

                                                  tcGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGAc                                       <  1:2305501/93‑1 (MQ=255)

                                                    gTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTg                                     <  1:3717586/93‑1 (MQ=255)

                                                    gTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTg                                     >  1:2779029/1‑93 (MQ=255)

                                                     ttGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGc                                    >  2:148899/1‑93 (MQ=255)

                                                     ttGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGc                                    >  2:3533297/1‑93 (MQ=255)

                                                      tGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGcc                                   >  1:4666596/1‑93 (MQ=255)

                                                         gCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCgag                                >  2:3091729/1‑93 (MQ=255)

                                                              gCAAACGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGt                           <  1:2759938/93‑1 (MQ=255)

                                                                aaaCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGc                         <  2:1138190/93‑1 (MQ=255)

                                                                 aaCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGCCTGCCGAGAACGCGCa                        >  2:2223492/1‑93 (MQ=255)

                                                                 aaCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCa                        >  2:4405492/1‑93 (MQ=255)

                                                                  aCGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAt                       >  1:642419/1‑93 (MQ=255)

                                                                   cGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtt                      <  1:4589015/93‑1 (MQ=255)

                                                                   cGCCGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtt                      <  1:3514394/93‑1 (MQ=255)

                                                                     ccGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtttt                    >  1:858895/1‑93 (MQ=255)

                                                                     ccGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAtttt                    <  1:3659780/93‑1 (MQ=255)

                                                                      cGGATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCAttttt                   >  2:3642575/1‑92 (MQ=255)

                                                                       ggATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGt                  >  2:1239270/1‑93 (MQ=255)

                                                                        gATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTa                 >  2:45512/1‑93 (MQ=255)

                                                                        gATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTa                 >  1:1346710/1‑93 (MQ=255)

                                                                        gATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTa                 <  2:4316476/93‑1 (MQ=255)

                                                                        gATGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTa                 <  1:11580/93‑1 (MQ=255)

                                                                         aTGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTAt                <  1:4207471/93‑1 (MQ=255)

                                                                          tGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg               >  2:4588697/1‑93 (MQ=255)

                                                                          tGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg               >  1:4643914/1‑93 (MQ=255)

                                                                          tGGTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATg               >  1:2467676/1‑93 (MQ=255)

                                                                           ggTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGc              >  1:2203515/1‑93 (MQ=255)

                                                                           ggTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGc              <  1:1819146/93‑1 (MQ=255)

                                                                           ggTTGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGc              <  2:1764601/93‑1 (MQ=255)

                                                                              tGGCGGCCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCa           >  2:1686953/1‑93 (MQ=255)

                                                                                gcggcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATc         <  1:3816025/93‑1 (MQ=255)

                                                                                 cggcCAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATcc        <  2:1258686/93‑1 (MQ=255)

                                                                                    ccAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAg     >  1:1133244/1‑93 (MQ=255)

                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc    >  2:460736/1‑93 (MQ=255)

                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc    >  2:1944369/1‑93 (MQ=255)

                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGc    >  2:3066152/1‑93 (MQ=255)

                                                                                     cAGGAAGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCACAGc    >  1:3599621/1‑93 (MQ=255)

                                                                                      aggaagATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCa   <  2:3049997/93‑1 (MQ=255)

                                                                                       ggaagaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAc  >  1:2247944/1‑93 (MQ=255)

                                                                                       ggaagaTATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAc  >  2:3371269/1‑93 (MQ=255)

TCACTGTGAAACGGCGGCGGCGCGTTATGGCTTAAGTGCGGTGGATATTCTCGTTGAGCTGGGCAAACGCCGGATGGTTGGCGGCCAGGAGGATATGATCGTTGACGTGGCGCTGGATCTGCGCAACAACAAATAATGATGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCAC  >  NC_000913/373970‑374149

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: `‑`