BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,221,280	T→C	L336P (CTG→CCG)	b3076 →	evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,221,280	0	T	C	100.0%	648.1 / NA	204	L336P (CTG→CCG)	b3076	evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (110/94); total (110/94)

Only 100 of 214 total aligned reads displayed.

ACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGC  >  NC_000913/3221190‑3221369

aCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAtgccgc                                                                                         >  1:2478565/1‑93 (MQ=255)

 cGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCa                                                                                        <  1:3281899/93‑1 (MQ=255)

  gTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCAc                                                                                       <  1:556476/93‑1 (MQ=255)

   ttCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCa                                                                                        >  2:241075/1‑91 (MQ=255)

   ttCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCa                                                                                        <  1:241075/91‑1 (MQ=255)

      tGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCg                                                                                   <  1:1052991/93‑1 (MQ=255)

      tGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCg                                                                                   >  2:4045022/1‑93 (MQ=255)

        gAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTc                                                                                 >  1:2916561/1‑93 (MQ=255)

        gAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTc                                                                                 <  2:643983/93‑1 (MQ=255)

         aaGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCa                                                                                >  2:4472531/1‑93 (MQ=255)

          aGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCaa                                                                               <  2:3375198/93‑1 (MQ=255)

           gtggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAAc                                                                              <  1:981604/93‑1 (MQ=255)

            tggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAAcc                                                                             >  2:1744600/1‑93 (MQ=255)

              gtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGt                                                                           <  2:4213790/93‑1 (MQ=255)

                 ccACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCAc                                                                        <  1:1722249/93‑1 (MQ=255)

                 ccACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCAc                                                                        <  2:1477833/93‑1 (MQ=255)

                  cacaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACg                                                                       >  2:4559755/1‑93 (MQ=255)

                  cacaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACg                                                                       <  1:3765527/93‑1 (MQ=255)

                   acaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGa                                                                      >  2:1022683/1‑93 (MQ=255)

                      aCGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACaa                                                                   >  1:1575221/1‑93 (MQ=255)

                          gTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGat                                                               >  1:1943769/1‑93 (MQ=255)

                          gTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGat                                                               >  2:2574264/1‑93 (MQ=255)

                             ggCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGatcat                                                            <  2:3779672/93‑1 (MQ=255)

                             ggCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGatcat                                                            >  2:2176088/1‑93 (MQ=255)

                              gCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGatcatc                                                           <  1:700614/93‑1 (MQ=255)

                                ttCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGc                                                         >  1:1203626/1‑93 (MQ=255)

                                ttCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGc                                                         >  2:2755734/1‑93 (MQ=255)

                                 tCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCa                                                        >  1:413343/1‑93 (MQ=255)

                                 tCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCa                                                        >  2:199889/1‑93 (MQ=255)

                                      gATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGc                                                   >  1:3079696/1‑93 (MQ=255)

                                      gATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGc                                                   >  2:770579/1‑93 (MQ=255)

                                      aatatCAAAGCGCGTGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCTCGGCGTAAACCGTCGCGTCACCGATCATCGCAAAGGc                                                   <  2:1317843/92‑1 (MQ=255)

                                        tatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcg                                                 >  1:447483/1‑93 (MQ=255)

                                        tatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcg                                                 <  1:2802301/93‑1 (MQ=255)

                                          tCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcgcg                                               >  1:1252643/1‑93 (MQ=255)

                                          tCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAACgcggcg                                               >  1:2740921/1‑93 (MQ=255)

                                           cagaGTGCGCTACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCgcgc                                              <  1:1135863/90‑1 (MQ=255)

                                           cAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCgcgc                                              >  2:3752245/1‑93 (MQ=255)

                                            aaaGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc                                             >  2:3529138/1‑93 (MQ=255)

                                            aaaGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc                                             >  2:2193372/1‑93 (MQ=255)

                                             aaGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCg                                            >  1:3096738/1‑93 (MQ=255)

                                              aGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGt                                           >  2:2651016/1‑93 (MQ=255)

                                              aGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGt                                           <  2:334765/93‑1 (MQ=255)

                                                 gcgcGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTgg                                        >  1:4595448/1‑93 (MQ=255)

                                                    cgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAAt                                     <  1:1320223/93‑1 (MQ=255)

                                                     gACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATg                                    <  1:3133521/93‑1 (MQ=255)

                                                      aCGGTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATgg                                   >  2:452601/1‑93 (MQ=255)

                                                       cagTCTGTTCTGGATCAATAACCGTTACGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGa                                  <  1:1499582/91‑1 (MQ=255)

                                                        ggTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGAt                                 <  2:117353/93‑1 (MQ=255)

                                                        ggTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGAt                                 <  1:3580952/93‑1 (MQ=255)

                                                        ggTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGAt                                 >  1:4018831/1‑93 (MQ=255)

                                                         gTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATc                                >  2:4363613/1‑93 (MQ=255)

                                                         gTCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATc                                <  1:2197658/93‑1 (MQ=255)

                                                          tCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcg                               >  1:322048/1‑93 (MQ=255)

                                                          tCTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcg                               >  2:1119950/1‑93 (MQ=255)

                                                           cTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgc                              <  2:1368593/93‑1 (MQ=255)

                                                           cTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgc                              <  2:1400703/93‑1 (MQ=255)

                                                           cTGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgc                              <  2:489445/93‑1 (MQ=255)

                                                            tGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgcg                             >  2:2385926/1‑93 (MQ=255)

                                                            tGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgcg                             >  1:1620343/1‑93 (MQ=255)

                                                            tGTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcgcg                             <  1:814991/93‑1 (MQ=255)

                                                             gTTCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGt                            <  2:121368/93‑1 (MQ=255)

                                                               tCTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAATGCCGCGCCGTTGGAATGGATCGTGTCg                          >  2:2921166/1‑93 (MQ=255)

                                                                cTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCga                         <  1:2624228/93‑1 (MQ=255)

                                                                cTGGATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCga                         <  2:1239000/93‑1 (MQ=255)

                                                                   gATCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaa                      <  2:4659509/93‑1 (MQ=255)

                                                                    aTCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaaa                     <  1:900354/93‑1 (MQ=255)

                                                                    aTCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaaa                     <  2:1075509/93‑1 (MQ=255)

                                                                    aTCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaaa                     <  1:3853643/93‑1 (MQ=255)

                                                                     tCAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAg                    <  2:365472/93‑1 (MQ=255)

                                                                      cAATAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGa                   <  1:1895748/93‑1 (MQ=255)

                                                                       aaTAACCGTTATGTGCTGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCGAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAACAt                  >  2:2104634/1‑93 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  >  2:3696430/1‑93 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  >  1:277655/1‑93 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  <  1:3449343/93‑1 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  >  2:855407/1‑93 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  <  1:10266/93‑1 (MQ=255)

                                                                       aaTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAt                  <  2:2495594/93‑1 (MQ=255)

                                                                        aTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtc                 <  2:3598755/93‑1 (MQ=255)

                                                                        aTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtc                 <  1:1434522/93‑1 (MQ=255)

                                                                        aTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtc                 >  2:1684458/1‑93 (MQ=255)

                                                                        aTAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtc                 <  2:2090892/93‑1 (MQ=255)

                                                                         tAACCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtct                >  2:2312836/1‑93 (MQ=255)

                                                                          aaCCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGTAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtctc               <  1:3303830/93‑1 (MQ=255)

                                                                          aaCCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtctc               >  1:88037/1‑93 (MQ=255)

                                                                           aCCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTcc              >  1:3998256/1‑93 (MQ=255)

                                                                           aCCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTcc              >  2:344867/1‑93 (MQ=255)

                                                                           aCCGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTAGAGAAAGATCTcc              >  2:4137379/1‑93 (MQ=255)

                                                                            ccGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCa             <  2:4256040/93‑1 (MQ=255)

                                                                             cGTTATGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAg            <  1:3715891/93‑1 (MQ=255)

                                                                                 aTGTGATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtga        >  2:590709/1‑93 (MQ=255)

                                                                                  tgtgATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgat       <  1:1659590/93‑1 (MQ=255)

                                                                                  tgtgATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgat       >  2:1550242/1‑93 (MQ=255)

                                                                                  tgtgATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgat       >  2:69139/1‑93 (MQ=255)

                                                                                   gtgATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgatg      <  1:222438/93‑1 (MQ=255)

                                                                                   gtgATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgatg      >  1:4322493/1‑93 (MQ=255)

                                                                                     gATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtgatgag    >  2:3857833/1‑92 (MQ=255)

                                                                                     gATGCCGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATATCCAGTTGATGaa    <  2:3418277/93‑1 (MQ=255)

                                                                                       tgccgcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCCGTTGATGaacc  >  1:586269/1‑91 (MQ=255)

                                                                                       tgccgcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAagc  <  2:3923531/93‑1 (MQ=255)

ACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGC  >  NC_000913/3221190‑3221369

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