BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,197,178	T→C	V35V (GTA→GTG)	glnE ←	GlnE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,197,178	0	T	C	92.7%	987.8 / 47.5	356	V35V (GTA→GTG)	glnE	GlnE
Reads supporting (aligned to +/- strand): ref base T (13/13); new base C (182/148); total (195/161)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 384 total aligned reads displayed.

TTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGA  >  NC_000913/3197088‑3197267

 tcgtcgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCCGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACt                                                                                        >  2:1752709/1‑93 (MQ=255)

 tcgtcgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACt                                                                                        <  2:8165528/93‑1 (MQ=255)

 tcgtcgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACt                                                                                        <  2:6621177/93‑1 (MQ=255)

 tcgtagGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACt                                                                                        <  1:7802133/93‑1 (MQ=255)

   gtcgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGa                                                                                      >  2:1064257/1‑93 (MQ=255)

      ggCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACtt                                                                                   <  1:6376118/93‑1 (MQ=255)

       gCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTc                                                                                  >  2:9791754/1‑93 (MQ=255)

        ccTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTcc                                                                                 >  2:2612819/1‑92 (MQ=255)

         cTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGc                                                                                >  1:4564951/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc                                                                               >  2:6438488/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               >  2:5976334/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               >  1:7122981/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               <  2:6893046/93‑1 (MQ=255)

           gcggcgGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCt                                                                              <  2:4981836/93‑1 (MQ=255)

           gcggcgGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCt                                                                              >  1:2070549/1‑93 (MQ=255)

           gcggcgGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCt                                                                              >  1:6804169/1‑93 (MQ=255)

              gcggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCtgtg                                                                           >  2:1498001/1‑93 (MQ=255)

               cggTTGCCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgc                                                                          <  1:9311121/93‑1 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCATCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  1:1838507/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  1:2840808/1‑93 (MQ=255)

                  ttGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCt                                                                       >  1:5538030/1‑93 (MQ=255)

                  gtGGGTTTCCAGTTCCGTCAGCCACTCTGGGTGGGCAATTACGCCGTCCTGGCAAAAATCACTAAAAGTAAGCAACGACTTCGCCCGGGCGCt                                                                       >  2:5135811/2‑93 (MQ=255)

                   tGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTa                                                                      <  1:8511582/93‑1 (MQ=255)

                    ggCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTaa                                                                     <  2:4900051/93‑1 (MQ=255)

                    ggCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTaa                                                                     >  1:7398076/1‑93 (MQ=255)

                      cTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGt                                                                   <  2:3093386/93‑1 (MQ=255)

                      cTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGt                                                                   <  1:6978362/93‑1 (MQ=255)

                       tttCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTg                                                                  >  2:8572930/1‑93 (MQ=255)

                        tgccAGTCCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGa                                                                 <  1:1324695/91‑1 (MQ=255)

                         tCCAGTTCCGTCAGCCACTCTGGATGCGCACTCACGCTGTACTGCACAAATTCTCTAAAAGTAAGCACTAACTTCGCCTGTGCGCTAAGTGAt                                                                <  2:8933645/93‑1 (MQ=255)

                          ccAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGAtt                                                               <  2:2465465/93‑1 (MQ=255)

                             gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                            >  1:4422810/1‑93 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           <  2:8348224/93‑1 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           <  1:1025311/93‑1 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           >  1:9874066/1‑93 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           <  1:9006596/93‑1 (MQ=255)

                               tCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCg                                                          >  2:2374803/1‑93 (MQ=255)

                               tCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCg                                                          <  2:2027451/93‑1 (MQ=255)

                                    cAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTaa                                                     <  2:1319904/93‑1 (MQ=255)

                                      gCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAg                                                   <  2:8508332/93‑1 (MQ=255)

                                      gCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAg                                                   >  1:9913974/1‑93 (MQ=255)

                                        cACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGc                                                 >  1:4008340/1‑93 (MQ=255)

                                        cACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGc                                                 >  1:2037517/1‑93 (MQ=255)

                                         aCTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGct                                                >  2:405870/1‑93 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               >  1:1518083/1‑93 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               <  1:7222730/93‑1 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               <  1:3196674/93‑1 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               <  2:7935472/93‑1 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGCACAAAATCACAAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               >  1:8730942/1‑93 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  2:9613824/93‑1 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  2:7413891/93‑1 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  1:413816/93‑1 (MQ=255)

                                             cggATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTgg                                            <  1:6293823/92‑1 (MQ=255)

                                              ggATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGc                                           >  1:9628074/1‑93 (MQ=255)

                                               gATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCa                                          >  1:6001083/1‑93 (MQ=255)

                                                 tGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                        <  2:5842126/93‑1 (MQ=255)

                                                    gcAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGc                                     >  1:8810075/1‑93 (MQ=255)

                                                      aaTCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTc                                   >  2:4309090/1‑93 (MQ=255)

                                                      aaTCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGCGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTc                                   >  1:5184315/1‑93 (MQ=255)

                                                        tCACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaa                                 >  2:9147803/1‑93 (MQ=255)

                                                         cACGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaac                                >  1:4102290/1‑93 (MQ=255)

                                                          aCGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaaca                               <  1:3736847/93‑1 (MQ=255)

                                                          aCGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaaca                               >  2:4645168/1‑93 (MQ=255)

                                                           cGCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaacaa                              <  1:2716178/93‑1 (MQ=255)

                                                            gCTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCaacaac                             <  2:9496733/93‑1 (MQ=255)

                                                             cTGTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACg                            >  1:4653632/1‑93 (MQ=255)

                                                               gTCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          <  2:5211731/93‑1 (MQ=255)

                                                                tCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                         >  1:3269783/1‑93 (MQ=255)

                                                                tCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                         >  2:9388105/1‑93 (MQ=255)

                                                                tCCTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                         >  2:7550044/1‑93 (MQ=255)

                                                                  cTGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTg                       >  1:3534877/1‑93 (MQ=255)

                                                                   tGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                      >  2:3789582/1‑93 (MQ=255)

                                                                   tGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                      >  2:1646021/1‑93 (MQ=255)

                                                                   tGCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGc                      >  2:8882371/1‑93 (MQ=255)

                                                                    gCACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGcc                     >  2:630480/1‑93 (MQ=255)

                                                                     cacaAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCa                    <  2:6474025/93‑1 (MQ=255)

                                                                      acaAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                   >  1:8069935/1‑93 (MQ=255)

                                                                      acaAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                   <  2:3279885/93‑1 (MQ=255)

                                                                       caAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGt                  <  1:560678/93‑1 (MQ=255)

                                                                        aaaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa                 <  2:2878911/93‑1 (MQ=255)

                                                                        aaaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa                 <  2:2998439/93‑1 (MQ=255)

                                                                         aaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                <  2:4509811/93‑1 (MQ=255)

                                                                          aaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAct               >  1:5917225/1‑93 (MQ=255)

                                                                           aTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCCACAACGGTCTGCCAGTActg              >  2:4198108/1‑93 (MQ=255)

                                                                           aTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg              <  1:2180215/93‑1 (MQ=255)

                                                                           aTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg              >  2:3672710/1‑93 (MQ=255)

                                                                            tCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgc             >  1:9696966/1‑93 (MQ=255)

                                                                            tCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgc             <  1:1258907/93‑1 (MQ=255)

                                                                            tCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGCCTGCCAGTActgc             >  2:7755707/1‑93 (MQ=255)

                                                                              aCTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg           <  1:7935459/93‑1 (MQ=255)

                                                                              aCTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg           >  1:390383/1‑93 (MQ=255)

                                                                               cTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGt          <  2:3912121/93‑1 (MQ=255)

                                                                                tAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAAGAACGGTCTGCCAGTACTGCTGTa         >  2:5562244/1‑93 (MQ=255)

                                                                                tAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTa         >  2:7111088/1‑93 (MQ=255)

                                                                                 aaaaGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTaa        <  1:4925769/93‑1 (MQ=255)

                                                                                   aagtaagcACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACg      >  2:2593919/1‑93 (MQ=255)

                                                                                    agtaagtaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACgg     <  1:10027937/93‑1 (MQ=255)

                                                                                    agtaagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACgg     <  1:2518466/93‑1 (MQ=255)

                                                                                       aagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa  <  2:7190337/93‑1 (MQ=255)

                                                                                       aagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa  >  2:8736827/1‑93 (MQ=255)

TTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGA  >  NC_000913/3197088‑3197267

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