BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,349,989	T→C	M136V (ATG→GTG)	glpT ←	sn‑glycerol‑3‑phosphate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,349,989	0	T	C	100.0%	530.7 / NA	167	M136V (ATG→GTG)	glpT	sn‑glycerol‑3‑phosphate permease
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (83/84); total (83/84)

Only 100 of 180 total aligned reads displayed.

TGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349899‑2350078

  tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       <  2:3041624/93‑1 (MQ=255)

  tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       <  1:1795729/93‑1 (MQ=255)

  ggcgcACAGTTCCACACTGACACAATGCCGCCACGTTCCTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       <  2:4204635/92‑1 (MQ=255)

   gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCt                                                                                      >  1:4057952/1‑93 (MQ=255)

    cgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCTCCAGTGCCCCATAGTACGACCACACGGCGGCCACCCCACCCCCTg                                                                                     <  2:5002431/93‑1 (MQ=255)

     gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    <  2:2775189/93‑1 (MQ=255)

     gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    >  2:5717856/1‑93 (MQ=255)

      cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  2:769833/1‑93 (MQ=255)

       acaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaa                                                                                  >  2:4480107/1‑93 (MQ=255)

        caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 >  2:2706038/1‑93 (MQ=255)

        caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 >  1:2048357/1‑93 (MQ=255)

        caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 >  1:3712714/1‑93 (MQ=255)

        caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 <  1:5786658/93‑1 (MQ=255)

         aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                >  2:6254414/1‑93 (MQ=255)

         aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                <  2:2475213/93‑1 (MQ=255)

            tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:6163067/93‑1 (MQ=255)

            tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:3348511/93‑1 (MQ=255)

            tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             >  1:2610574/1‑93 (MQ=255)

             ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            <  2:4917081/93‑1 (MQ=255)

                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         >  1:1552814/1‑93 (MQ=255)

                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         >  1:5637052/1‑93 (MQ=255)

                   tGACACAATGCCCCCACGTTCTTTCTGCGACCGCCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      <  1:526812/93‑1 (MQ=255)

                     acacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGa                                                                    >  1:3970393/1‑93 (MQ=255)

                       acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCGACCGCAGAGGaac                                                                  >  1:3605185/1‑93 (MQ=255)

                       acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaac                                                                  >  2:1726008/1‑93 (MQ=255)

                       acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaa                                                                  >  1:171144/1‑92 (MQ=255)

                        caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 <  1:5650965/93‑1 (MQ=255)

                        caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 >  1:645824/1‑93 (MQ=255)

                         aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  1:4669048/93‑1 (MQ=255)

                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               >  1:5413159/1‑93 (MQ=255)

                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               <  2:4165960/93‑1 (MQ=255)

                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               >  1:237183/1‑93 (MQ=255)

                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               <  2:5565631/93‑1 (MQ=255)

                           tGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  1:1016072/93‑1 (MQ=255)

                           tGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  2:6193065/93‑1 (MQ=255)

                                ccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACa                                                         <  2:4800959/93‑1 (MQ=255)

                                 cACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                                                        >  2:5973526/1‑93 (MQ=255)

                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCCCCCCACCCCCTGGAACCAACCGCAGAGGGACAACAGTACaaa                                                       >  1:1234977/1‑93 (MQ=255)

                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:4177482/1‑93 (MQ=255)

                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:2390791/1‑93 (MQ=255)

                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:5360389/1‑93 (MQ=255)

                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:2341567/1‑93 (MQ=255)

                                   cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                                                      >  2:4717347/1‑93 (MQ=255)

                                   cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                                                      >  1:4992564/1‑93 (MQ=255)

                                     ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    >  1:825851/1‑93 (MQ=255)

                                     ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    >  1:6082611/1‑93 (MQ=255)

                                      tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAATACAAACATc                                                   >  1:4883427/1‑93 (MQ=255)

                                       ttttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa                                                  <  2:206078/92‑1 (MQ=255)

                                        tttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCAc                                                 <  2:267613/93‑1 (MQ=255)

                                        tttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACACTACAAACATCAc                                                 >  1:192694/1‑93 (MQ=255)

                                          tctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg                                               <  1:1091829/93‑1 (MQ=255)

                                             gCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCaa                                            <  2:4517156/93‑1 (MQ=255)

                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  2:5580933/93‑1 (MQ=255)

                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           >  1:4707039/1‑93 (MQ=255)

                                               gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          >  1:5937204/1‑93 (MQ=255)

                                                accaccAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGc                                         <  2:1116102/93‑1 (MQ=255)

                                                 ccaccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCt                                        <  2:2516282/93‑1 (MQ=255)

                                                  caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       >  2:5476559/1‑93 (MQ=255)

                                                   accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  2:1696476/1‑93 (MQ=255)

                                                   accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  2:3448930/1‑93 (MQ=255)

                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  1:2848166/93‑1 (MQ=255)

                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  1:3589152/93‑1 (MQ=255)

                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  2:3110372/1‑93 (MQ=255)

                                                     caGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAc                                    >  2:1068058/1‑93 (MQ=255)

                                                      aGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACg                                   >  1:2782797/1‑93 (MQ=255)

                                                            ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACCACAGTACAAACATCACCGCCATGCTCGACGTCGccc                             >  2:1444594/1‑93 (MQ=255)

                                                            ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGccc                             <  2:690953/93‑1 (MQ=255)

                                                            ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGacc                             >  1:5185780/1‑93 (MQ=255)

                                                            ccATAGTACGACAACCCGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCCCCGCAATGCTCGACGCCGccc                             <  2:117156/93‑1 (MQ=255)

                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:2507371/93‑1 (MQ=255)

                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:2163696/93‑1 (MQ=255)

                                                               tAGTACGACCGCACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          <  1:3001533/93‑1 (MQ=255)

                                                               tAGTACGACCACACGGCGGCCGCCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          <  1:1407038/93‑1 (MQ=255)

                                                               tAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          <  1:1444941/93‑1 (MQ=255)

                                                               tAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          <  1:3148818/93‑1 (MQ=255)

                                                               tAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          >  1:4855273/1‑93 (MQ=255)

                                                                aGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATgg                         >  2:5361741/1‑93 (MQ=255)

                                                                  tACGACCACACGGCGGCCACCCCACCGCCGGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       <  1:2002586/93‑1 (MQ=255)

                                                                  tACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       <  2:511197/93‑1 (MQ=255)

                                                                  tACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       >  2:755089/1‑93 (MQ=255)

                                                                   aCGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcac                      >  1:3266428/1‑93 (MQ=255)

                                                                   aCGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcac                      <  2:5285583/93‑1 (MQ=255)

                                                                     gACCACACGGCGGCCACCCCTCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaa                    <  2:2873786/93‑1 (MQ=255)

                                                                      aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   >  2:2101340/1‑93 (MQ=255)

                                                                        cacaCGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGc                 >  2:3748525/1‑93 (MQ=255)

                                                                        cacaCGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGc                 <  2:1354364/93‑1 (MQ=255)

                                                                         acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAgac                >  2:583009/1‑91 (MQ=255)

                                                                         acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                >  2:3641456/1‑93 (MQ=255)

                                                                         acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                >  1:984928/1‑93 (MQ=255)

                                                                           acGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              <  2:5172832/93‑1 (MQ=255)

                                                                           acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAGCAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  1:4432207/1‑93 (MQ=255)

                                                                           acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  2:4368594/1‑93 (MQ=255)

                                                                            cggcggCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCAt             >  1:1057547/1‑93 (MQ=255)

                                                                            cggcggCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCAt             <  1:4465561/93‑1 (MQ=255)

                                                                               cggcCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaaa          <  1:6105039/93‑1 (MQ=255)

                                                                                  cccccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaa       <  1:557086/90‑1 (MQ=255)

                                                                                    cccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     <  2:4546274/92‑1 (MQ=255)

                                                                                    cccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     <  2:1954288/92‑1 (MQ=255)

                                                                                    accccaccccCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     >  1:4581836/1‑93 (MQ=255)

                                                                                       ccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  <  2:330547/93‑1 (MQ=255)

TGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349899‑2350078

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