BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,274	T→C	L537P (CTC→CCC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,274	0	T	C	100.0%	770.3 / NA	236	L537P (CTC→CCC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (132/104); total (132/104)

Only 100 of 240 total aligned reads displayed.

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

 aGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCgg                                                                                        <  2:5426609/93‑1 (MQ=255)

 aGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCgg                                                                                        >  1:3761272/1‑93 (MQ=255)

 aGTTATTCAGGTTAAGCCAAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAGCGTTACTCAGCAACCCgg                                                                                        <  1:3138307/93‑1 (MQ=255)

  gTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCggg                                                                                       >  1:6115798/1‑93 (MQ=255)

     aTTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgcgc                                                                                    >  2:5958504/1‑93 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   <  1:1877564/93‑1 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   >  1:4900154/1‑93 (MQ=255)

         aGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAt                                                                                <  2:4288654/93‑1 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               <  1:1755055/93‑1 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               >  2:339343/1‑93 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               <  2:787713/93‑1 (MQ=255)

          gcttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               <  2:5872165/91‑1 (MQ=255)

            ttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCaa                                                                             <  1:6258501/93‑1 (MQ=255)

             tAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAac                                                                            >  1:3044464/1‑93 (MQ=255)

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAacac                                                                          >  2:2069577/1‑93 (MQ=255)

                gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACt                                                                         <  2:503183/93‑1 (MQ=255)

                gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACt                                                                         <  1:5807886/93‑1 (MQ=255)

                gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACt                                                                         >  1:1536276/1‑93 (MQ=255)

                  cGAAATTCACCGACTATTACACCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGt                                                                       <  1:1307191/93‑1 (MQ=255)

                  cGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGt                                                                       >  1:3258495/1‑93 (MQ=255)

                    aaaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAt                                                                     <  2:162167/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    >  1:6147085/1‑93 (MQ=255)

                      aTTCACCGACTATTACATCCTCGCTTATAACAAACGCGGGAAATTACAACTCGCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   <  2:2831872/93‑1 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  >  2:4611085/1‑93 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  <  2:3623952/93‑1 (MQ=255)

                         cACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAg                                                                >  2:5137254/1‑93 (MQ=255)

                         cACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAg                                                                <  1:2850369/93‑1 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               <  1:960513/93‑1 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               >  2:5729169/1‑93 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               >  1:363460/1‑93 (MQ=255)

                           ccGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTg                                                              <  1:3159233/93‑1 (MQ=255)

                           ccGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTg                                                              >  1:5758856/1‑93 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             >  2:149051/1‑93 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             >  1:1575716/1‑93 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            <  1:2133649/93‑1 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            >  1:6297409/1‑93 (MQ=255)

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           <  2:2945436/93‑1 (MQ=255)

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           <  1:591996/93‑1 (MQ=255)

                                      aaCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaa                                                   >  1:1672259/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGCAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  2:3990708/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  1:2892026/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  <  2:860667/93‑1 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                <  1:4778209/93‑1 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                <  1:891017/93‑1 (MQ=255)

                                           cGCTTATAACAAACGCGGGAAATTACAACTCACCGTCACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTa                                              <  1:3306183/93‑1 (MQ=255)

                                            gCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAt                                             <  1:2616114/93‑1 (MQ=255)

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            <  1:3072066/93‑1 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:2462551/93‑1 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:1233248/93‑1 (MQ=255)

                                              ttAGAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:2569209/93‑1 (MQ=255)

                                               tataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                          <  1:6307744/93‑1 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         >  1:3570985/1‑93 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         <  2:3027296/93‑1 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         >  2:2063341/1‑93 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        >  1:193155/1‑93 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        <  1:1087462/93‑1 (MQ=255)

                                                  aaCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa                                       <  1:4634235/93‑1 (MQ=255)

                                                    cAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAAc                                     >  1:1186150/1‑93 (MQ=255)

                                                     aaaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGATAGCCATCAAACTTATTGGGGAACg                                    <  2:1072827/93‑1 (MQ=255)

                                                      aaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGa                                   >  2:3124689/1‑93 (MQ=255)

                                                         gcgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATTAAACTTATTGGGGAACGAGTa                                >  2:2173681/1‑93 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  2:493612/93‑1 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  1:678693/93‑1 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              >  2:389380/1‑93 (MQ=255)

                                                            ggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACCAGTAATg                             >  2:5351184/1‑93 (MQ=255)

                                                             gAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGt                            >  1:5489723/1‑93 (MQ=255)

                                                              aaaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           <  1:73760/93‑1 (MQ=255)

                                                              aaaTTACAACACACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           <  2:2011336/93‑1 (MQ=255)

                                                               aaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCg                          >  2:5164462/1‑93 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  1:4598881/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         >  2:2273141/1‑93 (MQ=255)

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        <  2:5023269/93‑1 (MQ=255)

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        <  2:5594838/93‑1 (MQ=255)

                                                                  tACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                       <  2:487462/93‑1 (MQ=255)

                                                                     aaCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                    >  2:381306/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                   >  1:3718259/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTAGATGAGCa                   >  1:4370065/1‑93 (MQ=255)

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               >  2:1736084/1‑93 (MQ=255)

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               >  2:1441917/1‑93 (MQ=255)

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               <  1:3977773/93‑1 (MQ=255)

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               <  1:3365669/93‑1 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              >  2:5929657/1‑93 (MQ=255)

                                                                            cGTTACTCAGCATCCCGCGCGCACATCATCTCTGTATTTGAGTGGTAGCAAGCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCa             <  2:3224529/93‑1 (MQ=255)

                                                                            cGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCa             <  1:4233730/93‑1 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            <  1:3546191/93‑1 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            >  1:1628446/1‑93 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          <  1:2244549/93‑1 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          >  2:2434682/1‑93 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          <  1:2743206/93‑1 (MQ=255)

                                                                                aCTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCt         <  2:1840192/93‑1 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  2:5252777/1‑93 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  1:427692/1‑93 (MQ=255)

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       <  2:6158611/93‑1 (MQ=255)

                                                                                    aGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt     <  1:802717/93‑1 (MQ=255)

                                                                                     gCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAtt    >  2:406432/1‑93 (MQ=255)

                                                                                      cAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTa   <  1:5891239/93‑1 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  >  1:3833742/1‑93 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  >  2:4447223/1‑93 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  >  2:4110729/1‑93 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  >  1:4994169/1‑93 (MQ=255)

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

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