BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,170	T→C	I502I (ATT→ATC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,170	0	T	C	100.0%	935.0 / NA	282	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (146/136); total (146/136)

Only 100 of 289 total aligned reads displayed.

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

   tatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAac                                                                                      <  1:5424807/93‑1 (MQ=255)

    atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                                     >  2:2240409/1‑93 (MQ=255)

    atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                                     >  1:2021620/1‑93 (MQ=255)

      tCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                                   <  1:3177582/93‑1 (MQ=255)

       cAGTTAGTTGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  <  1:2822615/93‑1 (MQ=255)

       cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  <  1:2496683/93‑1 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                <  2:5950325/93‑1 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                <  1:6326782/93‑1 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                >  2:1319527/1‑93 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               <  2:5451072/93‑1 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              <  2:2909412/93‑1 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              >  1:1945389/1‑93 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              >  1:648883/1‑93 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              >  2:1658733/1‑93 (MQ=255)

            aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACgg                                                                             <  1:757234/93‑1 (MQ=255)

               gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGt                                                                          >  2:1974882/1‑93 (MQ=255)

                  ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                                       <  1:4052736/93‑1 (MQ=255)

                  ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                                       <  2:2122694/93‑1 (MQ=255)

                   tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                                      <  1:4796124/93‑1 (MQ=255)

                   tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                                      >  1:4570656/1‑93 (MQ=255)

                     ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCa                                                                    <  1:1271307/93‑1 (MQ=255)

                      cGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAg                                                                   <  2:6657359/93‑1 (MQ=255)

                       gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                                  >  2:2825410/1‑93 (MQ=255)

                       gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                                  >  1:1627356/1‑93 (MQ=255)

                        ttATTCGACCAGCGGATATTTTAATTTCACTGATACAACATACAGTCGAATGAATGGATACAACATCGAAACACGGGACGGAGCTATTCAGGt                                                                 >  1:942568/1‑93 (MQ=255)

                         tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                                                <  2:5135579/93‑1 (MQ=255)

                         tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                                                <  1:5341696/93‑1 (MQ=255)

                          aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                                               <  2:919303/93‑1 (MQ=255)

                           ttCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCAAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                                              <  2:2561484/93‑1 (MQ=255)

                               aCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGCTATTCAGGTTAAGCCg                                                          >  2:5241162/1‑93 (MQ=255)

                                ccAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTCAGCCGa                                                         >  1:3080165/1‑93 (MQ=255)

                                ccAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGa                                                         <  1:4193810/93‑1 (MQ=255)

                                  aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                                       >  2:5080032/1‑93 (MQ=255)

                                  aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAACGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                                       <  2:1088849/93‑1 (MQ=255)

                                   gCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAt                                                      >  1:2440220/1‑93 (MQ=255)

                                   gCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAt                                                      <  2:4378721/93‑1 (MQ=255)

                                    cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                     <  2:3099824/93‑1 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGCCGGAGGTATTCAGGTTAAGCCGACATTc                                                    >  2:487572/1‑93 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                    >  1:349685/1‑93 (MQ=255)

                                        tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                 >  2:3640513/1‑93 (MQ=255)

                                        tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGCTATTCAGGTTAAGCCGAAATTCAcc                                                 >  2:1352224/1‑93 (MQ=255)

                                         atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAACTTCACCg                                                >  2:6034334/1‑93 (MQ=255)

                                          ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTGACCGa                                               >  1:2647928/1‑93 (MQ=255)

                                          ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                                               >  1:1097031/1‑93 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              >  2:1730744/1‑93 (MQ=255)

                                             tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTa                                            >  1:5816567/1‑93 (MQ=255)

                                             tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTa                                            <  2:1455005/93‑1 (MQ=255)

                                               aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                          <  2:5527534/93‑1 (MQ=255)

                                                tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACAGAGGACGGAGTTATTCAGGTTGAGCCGAAATTCCCCGACTATTa                                         >  1:6621843/1‑93 (MQ=255)

                                                 ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                        <  1:13879/93‑1 (MQ=255)

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAACTTCACCGACTATTACAAc                                     >  2:6187565/1‑93 (MQ=255)

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                                     <  2:4964418/93‑1 (MQ=255)

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                                     >  1:4895667/1‑93 (MQ=255)

                                                     cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAcc                                    >  2:785726/1‑93 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   >  2:5303893/1‑93 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   >  1:6057229/1‑93 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   >  1:555971/1‑93 (MQ=255)

                                                       taTACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  1:3741307/92‑1 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  >  1:2165187/1‑93 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  2:2316437/93‑1 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  >  2:3490971/1‑93 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                <  2:153396/93‑1 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                <  1:6681353/93‑1 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                <  2:6231827/93‑1 (MQ=255)

                                                           caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt                              <  2:662539/93‑1 (MQ=255)

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             <  1:6240380/93‑1 (MQ=255)

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             >  1:2124654/1‑93 (MQ=255)

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             <  2:324547/93‑1 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat                            >  2:4764894/1‑93 (MQ=255)

                                                              caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           <  1:1805235/93‑1 (MQ=255)

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          >  2:3839688/1‑93 (MQ=255)

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          <  2:650052/93‑1 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         <  2:217501/93‑1 (MQ=255)

                                                                 aCAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGTTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACa                        <  2:4699189/93‑1 (MQ=255)

                                                                  cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaa                       <  2:3436312/93‑1 (MQ=255)

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      <  2:1238618/93‑1 (MQ=255)

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      <  1:4439956/93‑1 (MQ=255)

                                                                    gTCGAATGGATGGCTACAAAATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     >  1:6568565/1‑93 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     <  2:3309103/93‑1 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     >  1:1040243/1‑93 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     >  1:3964956/1‑93 (MQ=255)

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  >  1:1739584/1‑93 (MQ=255)

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  >  1:5762209/1‑93 (MQ=255)

                                                                         atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCggg                >  1:6137221/1‑93 (MQ=255)

                                                                           gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaa              <  2:4941812/93‑1 (MQ=255)

                                                                            aatgGCTACAACATCGAAACACAGGACGGTGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaaa             <  1:3445512/93‑1 (MQ=255)

                                                                            aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaaa             >  1:1472108/1‑93 (MQ=255)

                                                                            aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaaa             <  2:3933567/93‑1 (MQ=255)

                                                                             atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt            <  2:5749011/93‑1 (MQ=255)

                                                                              tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           >  1:357476/1‑93 (MQ=255)

                                                                               ggCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTa          >  1:2965601/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       <  2:472445/93‑1 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      >  2:4066043/1‑93 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      >  1:3879449/1‑93 (MQ=255)

                                                                                    caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACt     >  2:2043374/1‑93 (MQ=255)

                                                                                    caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACt     <  2:5615277/93‑1 (MQ=255)

                                                                                    caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACt     <  2:420055/93‑1 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  2:5844259/1‑93 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  1:5205172/1‑93 (MQ=255)

                                                                                       caTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc  <  1:1974882/93‑1 (MQ=255)

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

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